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EDITORIAL 

El cambio climático global es tal vez el más grande desafío que enfrenta la humanidad en 
el siglo XXI. La evidencia científica acumulada confirma la asociación entre la realización 
de actividades humanas tales como el consumo de energía fósil o el cambio de uso de 
suelo y las crecientes emisiones de gases de efecto invernadero (gei). Existe evidencia 
científica que demuestra la estrecha relación entre el aumento continuo de emisiones de 
gei y los impactos climáticos, en particular el aumento paulatino de la temperatura, modi- 
ficaciones en los patrones de precipitación, cambios en la intensidad o en la frecuencia 
de eventos climáticos extremos y el alza del nivel de mar. Valorar el efecto del cambio en 
el clima global y precisar los costos económicos que representa es una tarea indispen- 
sable para poder instrumentar una estrategia global de mitigación. El análisis económico 
del cambio climático es un tema complejo que requerirá el concurso de investigadores de 
diversas disciplinas y demandará la creación de metodologías que permitan entender sus 
impactos y consecuencias. 

En este número de la revista se presentan artículos de investigación que seguramente 
contribuirán de manera sólida a entender este fenómeno tan preocupante. Uno de estos 
trabajos aborda el tema de los seivicios ambientales hidrológicos bajo escenarios de 
cambio climático en una región del centro de la República Mexicana. El análisis presen- 
tado señala que los servicios ambientales que proveen los ecosistemas están en función, 
entre otras, de las condiciones climáticas predominantes en un determinado lugar, así 
como de la estructura y composición de los tipos de vegetación. El objetivo de este impor- 
tante estudio fue simular y cuantificar los impactos que un posible cambio climático puede 
ejercer sobre la regulación hídrica y la capacidad de recarga de acuíferos en un Parque 
Nacional. Se realizó la modelación de la distribución geográfica de las variables tempe- 
ratura y precipitación, lo que permitió evaluar el comportamiento h id rico y la capacidad de 
infiltración a partir del balance de humedad. Los resultados muestran una disminución en 
la recarga al acuífero en importantes porcentajes para las próximas décadas. Estos resul- 
tados deben de considerarse para fines de planeación de medidas de adaptación y miti- 
gación, a fin de que se impida el deterioro para que el parque nacional continúe otorgando 
servicios ambientales ante una amenaza como el cambio climático. 

Otro trabajo trata sobre la reconstrucción de 350 años de precipitación para una región 
en el noroeste de México. Como parte de este estudio se desarrollaron series de tiempo 
dendro-cronológicas de madera temprana, tardía y anillo total de Pseudotsuga menziesii 
con una longitud de tres siglos y medio, con el objetivo de analizar la variabilidad hidro- 
climática histórica de la región. La reconstrucción de precipitación invierno-primavera 
indica gran variabilidad interanual, decenal y multidecenal de los patrones de precipita- 
ción para la región. Las últimas tres décadas del siglo XX y los años corrientes de la 
primera década del siglo XXI indican un periodo de intensa sequía para la región, con 
impactos en lo ecológico y socioeconómico aun no cuantificados. 

Otro estudio analiza la influencia de la edad en el crecimiento de una especie forestal de 
crecimiento rápido conocida como vinal, en una provincia de Argentina. Este trabajo tuvo 
como objetivos analizar las posibilidades de crecimiento de árboles individuales en 
función de la edad, estimar el tiempo de tránsito y el turno biológico. Aunque se generaron 
resultados importantes para conocer las variables de interés, se destaca que se requiere 
continuar estudiando el comportamiento de la especie en otras condiciones de creci- 
miento. 
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Editorial 


También de Sudamérica, se incluye un trabajo que tuvo como objetivo evaluar el poten- 
cial energético de dos especies de la selva del Amazonas, Brasil. En esta región del 
mundo hay un gran consumo de madera por parte de las comunidades locales, tanto para 
producción de carbón como para uso como leña para la producción de ladrillos, activi- 
dades ambas, dicho sea de paso , con fuerte generación de sustancias contaminantes a 
la atmósfera y el suelo. Los resultados de este trabajo indican que las dos especies estu- 
diadas se pueden cultivar y utilizar con ventaja para los fines de producción energética, 
disminuyendo con las plantaciones, la presión sobre los bosques naturales, lo que es un 
resultado positivo, después de todo. 

Otro manuscrito trata sobre la anatomía, índices físicos e hidráulicos de la madera de una 
importante especie que crece en el este de México, entre otros lugares. Se presentan 
resultados sobre Indices de vulnerabilidad, de resistencia al colapso y conductividad 
hidráulica. Aunque la variación en estas características es baja, explican la estabilidad 
dimensional de la albura y el duramen a los cambios de régimen hídrico ocasionados por 
la sequía y la lluvia. Los índices de resistencia al colapso y de vulnerabilidad muestran un 
xilema muy vulnerable a embolismos, pero con mucho éxito para desarrollarse en 
ambientes perturbados. 

El número se completa con un estudio de la distribución espacial del muérdago enano en 
una de las más altas montañas de México, decretada hace tiempo Parque Nacional, reco- 
nociendo su importancia como sitio de recarga de los acuíferos, así como el valor de sus 
recursos naturales, pero que hoy es, sin embargo, escenario de una acelerada depreda- 
ción tanto humana como natural. El muérdago es una de las plantas parásitas más impor- 
tantes debido a que causa serios perjuicios sobre otros vegetales y que se ha convertido 
en un problema reciente debido a los niveles de infestación del parásito. En este trabajo 
se obtuvieron datos que permitieron conocer la distribución espacial de este parásito y 
mapas sobre su densidad en la zona de estudio. El resultado obtenido permitirá realizar 
programas de manejo de este parásito que sea sustentada por información confiable y 
que lleve a pretender un control adecuado de este problema. 

La función de revistas como Madera y Bosques es poner a la disposición de sus lectores, 
la información relevante sobre temas de importancia y trascendencia y así contribuir a 
entender mejor los fenómenos naturales y los inducidos por el hombre. Los artículos de 
este número seguramente ayudarán a cumplir esta función de diseminación del conoci- 
miento sobre los ecosistemas forestales y los bienes y servicios que se derivan de los 
mismos. 


Raymundo Dávalos Sotelo 


Editor 
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IARTÍCULO DE INVESTIGACIÓN! 


Servicios ambientales hidrológicos 
bajo escenarios de cambio climático 
en el Parque Nacional “El Chico”, Hidalgo 

Hydrologic environmental Services under climate change 
scenarios in the “El Chico” National Park, Hidalgo, México 

Alejandro Ismael Monterroso-Rivas, Jesús David Gómez-Díaz 

y Juan Ángel Tinoco-Rueda* 

RESUMEN 

Los servicios ambientales que proveen los ecosistemas están en función, entre otros aspectos, de 
las condiciones climáticas predominantes en un determinado lugar, así como de la estructura y composi- 
ción de los tipos de vegetación. Por lo anterior, el principal objetivo del presente estudio fue simular y cuan- 
tificar los impactos que un posible cambio climático puede ejercer sobre la regulación hídrica y la capa- 
cidad de recarga de acuíferos en el Parque Nacional "El Chico", Hidalgo. Se realizó la modelación de la 
distribución geográfica de las variables temperatura y precipitación tomando como periodo base 
1961-1990 y se aplicaron las razones de cambio obtenidas de los modelos de cambio climático nortea- 
mericano e inglés para el escenario A2 y los años 2020 y 2050. El comportamiento hídrico y la capacidad 
de infiltración en los diferentes escenarios se evaluaron a partir del balance de humedad obtenido con la 
metodología de Thornthwaite III modificado. Los resultados muestran una disminución en la recarga al 
acuífero aplicando el modelo inglés en 26% y 32% para los periodos 2020 y 2050, respectivamente. El 
modelo norteamericano estima un aumento en la infiltración de 6% y 22% para el 2020 y 2050, respecti- 
vamente. El periodo de crecimiento no presenta variaciones significativas en los escenarios estudiados. 

El modelo ingles prevé condiciones más drásticas, por lo cual debe de considerarse para fines de planear 
medidas de adaptación y mitigación, a fin de que se impida el deterioro del parque y que éste continúe 
otorgando servicios ambientales amenazados por el cambio climático. 

PALABRAS CLAVE: 

Cambio climático, escurrimiento, estado de Hidalgo, infiltración, servicios ambientales. 


ABSTRACT 

The environmental Services provide ecosystems depend on predominant dimate conditions of 
local scales, as well as in the vegetation structure and composition. The main goal of the present study 
was to simúlate and quantify the impacts that a possible climate change have affect on the hydrologic 
regulation and on the water infiltration capacity in the National Park “El Chico”, Hidalgo. It was modeled 
the geographic distribution of climatic variables like temperatura and precipitation, taking as baseline 
scenario the period of 1961-1990. It were applied the North American and English climate change 
models for the scenarios A2 and the years 2020 and 2050. The water balance was estimated accor- 
ding to the Thornthwaite modified or third versión in order, to evalúate the hydrologic behavior and the 
infiltration capacity in different scenarios. The main results obtained by followins the english model 
show a decrease in aquifer recharge in 26% and 32% for 2020 and 2050, respectively. The North 
American model considers an increase in the infiltration rate of 6% and 22% for the 2020 and 2050, 
respectively. The growing period did not show significant variations in the studied scenarios. The 


Universidad Autónoma Chapingo, Departamento de Suelos. Km. 38.5 Carretera México-Texcoco, CP 56230, 
Texcoco, Edo. de México. Email: aimrivas@correo. chapingo. mx 
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English model pointed out more drastic 
conditions, because of it that musí be 
considered to plan adoptive and mitiga- 
tion measures to prevent park deteriora- 
tion and ¡t can continué providing envi- 
ronmental Services deterioration that 
climate change put in risk. 

KEY WORDS: 

Climate change, runoff, Hidalgo State, 
infiltration, environmental Services. 


INTRODUCCIÓN 

Muchas son las definiciones que se 
le han dado al concepto de servicios 
ambientales o ecosistémicos, aunque la 
mayoría hacen siempre referencia a 
términos de la composición, estructura y 
función de los ecosistemas. Dentro de un 
ecosistema, su composición se refiere a 
todo aquello que se encuentra dentro del 
área funcional (Meffe, 2002), como por 
ejemplo la composición genética, las 
poblaciones, las especies, las comuni- 
dades y los paisajes. Mientras que la 
estructura de un ecosistema se refiere a 
las especies contenidas dentro de éste, 
su masa y su arreglo (Westman, 1977), 
es decir, cómo están distribuidas en el 
espacio y tiempo. Algunos ejemplos son 
la estructura genética, de poblaciones y 
de los hábitats, así como el comporta- 
miento del paisaje son algunos ejemplos. 

La función de un ecosistema es la 
forma en que los componentes del 
sistema interactúan (Westman, 1977). Se 
refiere al qué es lo que realizan, por 
ejemplo, los procesos genéticos, los 
procesos demográficos, las interacciones 
específicas así como los procesos y 
disturbios dentro del paisaje. Es la parte 
dinámica del ecosistema, por lo que sus 
procesos han sido entendidos como los 
servicios que presta el ecosistema o 
servicios ambientales (Westman, 1977; 
MA, 2003). El hombre obtiene algún tipo 
de beneficio a partir de estas funciones a 
cualquier escala, la mayoría no tiene 


precios de mercado y aún no pueden ser 
reemplazos por algún tipo de tecnología 
(Costanza, 1997). 

Las primeras referencias sobre la 
función y los servicios de los ecosistemas 
datan de mediados de los años sesentas 
y principios de los setentas (De Groot, et 
al., 2002). A partir de entonces han apare- 
cido una gran cantidad de publicaciones 
sobre el beneficio que los ecosistemas 
naturales proveen a la sociedad humana 
(Daily, 1997; Costanza, et al., 1997; de 
Groot, et al., 2002 y Lobo, 2001). En 
México, la Secretaría de Medio Ambiente 
y Recursos Naturales (semarnat, 2003) 
los define como los beneficios intangibles 
que los diferentes ecosistemas ponen a 
disposición de la sociedad, ya sea de 
manera natural o por medio de su manejo 
sustentable. 

Los beneficios que ofrecen los 
bosques a nivel local, regional y en 
algunos casos a nivel mundial, derivan de 
su valor como fuente de abastecimiento 
de agua, centros de diversidad biológica, 
origen de diversos productos madereros y 
no madereros, lugar de recreación y esta- 
bilizadores del suelo frente a los procesos 
erosivos. En general, los bosques se han 
asociado con diversos servicios ambien- 
tales al nivel de cuenca hidrológica, 
destacando entre otros: 1) regulación de 
los flujos de agua, 2) conservación de la 
calidad del agua, 3) control de la erosión 
y sedimentación, 4) reducción de la salini- 
zación del suelo, 5) regulación del nivel 
freático y, 6) conservación de hábitats 
acuáticos (García et al., 2004; Manson, 
2004 y Pattanayak, 2004). Myers (1997) 
señala que los bosques, por su particular 
estructura y abundante follaje, experi- 
mentan un “efecto esponja” al absorber la 
humedad proveniente de las precipita- 
ciones y liberarla gradualmente a las 
corrientes hídricas. El mismo autor señala 
que un bosque cerrado puede interceptar 
35% de la precipitación, mientras que un 
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bosque abierto menos del 20% y una 
plantación forestal sólo 12%. La mayor 
parte del agua que es interceptada por las 
hojas y el follaje forestal, disminuye su 
velocidad y reduce su tamaño, lo que 
hace posible que ésta caiga con menor 
fuerza y se evite así el desprendimiento 
del suelo, permitiendo una mejor percola- 
ción o una escorrentía gradual hacia las 
corrientes de agua. 

Uno de los principales servicios 
ambientales que ha influido en el éxito o 
fracaso de las sociedades es el abasteci- 
miento de agua potable, debido a su 
importancia tanto para el consumo 
humano como para la producción de 
alimentos. Sin embargo, el aumento de 
las áreas urbanas, la deforestación y el 
cambio de uso de suelo han provocado 
una disminución significativa en la calidad 
y cantidad de dichos servicios, en espe- 
cial el abastecimiento de agua (Maderey 
y Carrillo, 2005), por lo que se perfila 
como uno de los retos ambientales 
futuros más importantes para el ser 
humano (Manson, 2004; ipcc, 2007). 

Aunado a los factores antes mencio- 
nados, el cambio climático es otra 
amenaza a la estabilidad en la distribu- 
ción del agua en los ecosistemas (ipcc, 
2007) ya que no se tiene un panorama 
específico de cuáles serán las conse- 
cuencias que un posible cambio en el 
sistema climático pueda ejercer sobre el 
medio ambiente. Sin embargo, se han 
realizado simulaciones que pronostican 
resultados con un alto grado de incerti- 
dumbre, como por ejemplo largos 
periodos de sequía o incremento en la 
cantidad y periodo de lluvias (ipcc, 2007). 

En este sentido, Arnell (1999) simuló 
bajo el modelo HadCM3, los cambios en 
el escurrimiento medio anual para el año 
2050 para todo el mundo. En dicho 
trabajo se presenta que México se 
encuentra en rangos de van de menos 


250 a más de 25 mm/año, teniendo la 
mayor parte del territorio nacional zonas 
con cambios negativos. Rosegrant et al. 
(2002) establece que para el año 2025 la 
extracción mundial de agua para uso 
industrial, doméstico y ganadero aumen- 
tará en más del 50%, lo que limitará de 
manera considerable el riego para la 
producción de alimentos. También en 
México, Manson (2004) revisó los servi- 
cios hidrológicos brindados por estos 
ecosistemas, incluyendo la captación de 
agua y la prevención de ciclos de inunda- 
ción y sequía, así como la conservación 
de los suelos, la regulación del clima 
regional y la reducción del azolve de los 
cauces de los ríos, proponiendo final- 
mente el establecimiento de esquemas 
de pago por servicios ambientales como 
un mecanismo prometedor para 
aumentar la cobertura boscosa y favo- 
recer el manejo sustentable de los 
recursos naturales en el país. 

Es por lo anterior que surge la nece- 
sidad de evaluar el posible comporta- 
miento de los servicios ambientales hidro- 
lógicos que proveen los ecosistemas 
forestales frente a escenarios de cambio 
climático, ya que de su análisis se contará 
con más y mejores elementos para la 
toma de decisiones desde el presente, 
pero enfocado al mantenimiento de dicho 
servicio en un futuro no muy lejano. Para 
ello, en la presente investigación se plan- 
tearon los siguientes objetivos: 

1) Elaborar cartografía actualizada 
de las variables climáticas de temperatura 
y precipitación en el Parque Nacional El 
Chico, Hidalgo. 

2) Obtener y aplicar las razones de 
cambio sugeridas por los modelos de 
cambio climático norteamericano (GFDL- 
R30_A2) e inglés (HadCM3_A2) para el 
escenario socioeconómico A2 y los 
periodos 2020 y 2050. 
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3) Evaluar el posible comporta- 
miento del servicio ambiental hidrológico 
mediante balances de humedad en el 
suelo, aplicando la metodología de 
Thornthwaite III modificado y conside- 
rando condiciones actuales y futuras de 
cambio climático. 

METODOLOGÍA 
Área de estudio 

El Parque Nacional “El Chico”, fue 
decretado en 1982 como área natural 
protegida con una superficie de un poco 
más de 2 700 hectáreas. Se ubica al 
centro del estado de Hidalgo (Figura 1) e 
incluye parte de los municipios de Mineral 
del Chico, Mineral del Monte y Pachuca 
de Soto, la capital del estado. Presenta 
un rango altitudinal que va de los 2 000 a 


poco más de 3 000 metros sobre el nivel 
del mar, derivando en algunas zonas con 
fuertes pendientes y paisajes exube- 
rantes. De acuerdo con el programa de 
manejo (conanp, 2006), el parque “tiene 
especial relevancia dado que el macizo 
montañoso constituye una mínima frac- 
ción del parteaguas que separa a los 
sistemas hidrológicos correspondientes a 
la gran cuenca del río Pánuco y a la 
cuenca del Valle de México”. En este 
sentido, el parteaguas corre en sentido 
este-oeste y supera altitudes de 3 050 
msnm. La vertiente orientada hacia el 
norte abarca la mayor superficie del 
parque, 2 273 hectáreas que significan 
83%, mientras que la expuesta al sur 
comprende 466 hectáreas, alrededor del 
17% respecto al área. La vegetación 
predominante (coede, 2007) es de 
bosque de oyamel (1 856 ha), bosque de 


9B'45XTW SBMZCTW 



Figura 1. Localización del área de estudio. 
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pino-encino (403 ha), bosque de oyamel- 
encino (95 ha), bosque de tlaxcal o 
táscate (63 ha), pastizal (46 ha), bosque 
de cedro (27 ha) y bosque de encino (25 
ha). El mismo autor señala que la flora 
está constituida por 545 especies, 264 
géneros y 73 familias, de las cuales 4 
especies son sujetas a protección espe- 
cial, una especie amenazada y otra en 
peligro de extinción. 

Delimitación de áreas de influencia 
climática 

Las áreas de influencia climática son 
unidades geográficas que poseen valores 
de temperatura media anual y precipita- 
ción media anual únicos, se construyen a 
partir de isoyetas e isotermas medias 
anuales, para lo cual, primero se selec- 
cionaron las estaciones meteorológicas 
ubicadas en la zona de estudio y en áreas 
aledañas para ser utilizadas en la elabo- 


ración de la cartografía (Figura 2). La 
búsqueda de información climática y su 
condensación se llevó a cabo con ayuda 
de las normales climatológicas publi- 
cadas por el Servicio Meteorológico 
Nacional (smn, 2000). Se obtuvieron los 
datos promedio mensual y anual de preci- 
pitación y temperatura para el periodo 
1961-1990, conformando éste el año 
base. Para el trazo de las isoyetas se 
utilizó la metodología propuesta por 
Gómez et al. (2008) que utiliza como 
información base espaciomapas y/o 
imágenes de satélite, modelos digitales 
de elevación, curvas de nivel y los datos 
de precipitación de las estaciones meteo- 
rológicas de la zona de estudio. Utili- 
zando el software de Sistemas de Infor- 
mación Geográfica ArcView v.3.2 se 
sobrepusieron estos productos cartográ- 
ficos y se trazaron las isoyetas tomando 
en cuenta la condición de la vegetación 
como un indicador de la cantidad de 
precipitación, la dinámica de circulación 



98°50'W 


98°40'W 


98 * 30 ^ 


Figura 2. Localización de las estaciones meteorológicas utilizadas. 
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de los vientos con los que se asocia la 
entrada de humedad a la zona y los dife- 
rentes fenómenos meteorológicos que 
afectan al Parque a lo largo del año. Lo 
anterior para conformar el mapa de preci- 
pitación media anual y mensual. Para 
obtener las isotermas se crearon doce 
modelos mensuales y uno anual de regre- 
sión lineal simple para observar el 
comportamiento de la temperatura con 
referencia a la altitud (Gómez et al., 2007 
y 2008), utilizando un grupo de esta- 
ciones meteorológicas ubicadas en la 
misma región climática de la zona de 
estudio. Aplicando los modelos obtenidos 
se trazaron las isolíneas de temperatura 
por cada rango altitudinal, conformando 
así los doce mapas de temperatura 
mensual y anual. 

Finalmente, se obtuvieron las áreas 
de influencia climática (aic) al sobreponer, 
mediante el programa de Sistemas de 
Información Geográfica ArcView v.3.2, las 
isoyetas con las isotermas anuales. A 
cada una de las aic resultantes se le 
asignó su correspondiente base de datos 
de temperatura y precipitación mensual. 

Aplicación de modelos de cambio 
climático 

A la base de datos de las áreas de 
influencia climática obtenidas en el paso 
anterior, se les aplicaron las razones de 
cambio a la temperatura (°C) y precipita- 
ción (%) obtenidos por los modelos 
generales de circulación (mgc) nortea- 
mericano (GFDL-R30) e inglés 
(HadCM3) para los periodos 2020 y 
2050, de acuerdo al escenario socioeco- 
nómico A2,en el que se considera que el 
desarrollo de la economía global estará 
basada más en criterios de obtención de 
beneficios económicos que en la preser- 
vación del medio ambiente y que 
seguirá la disparidad entre los países 
desarrollados y los en vías de 


desarrollo. Las tasas de cambio de los 
mgc fueron obtenidas del Canadian 
Institute for Climate Studies (cíes, 
2007). 

Se decidió aplicar los modelos 
GFDL-R30 y HadCM3 debido a que esta- 
blecen escenarios diferentes. En la 
variable temperatura, ambos modelos 
predicen incrementos, sin embargo, en la 
variable precipitación el modelo gfdl 
establece aumentos en la cantidad de 
lluvia, y por el contrario, el modelo 
HadCM establece mermas en esta misma 
variable (Gómez et al. 2007, ipcc, 2007). 
Por su parte, el escenario A2 se eligió 
para esta investigación, el cual se basa 
en el supuesto de que el desarrollo global 
seguirá una tendencia tal y como se 
presenta actualmente, diferenciándose 
del escenario B, el cual establece que el 
desarrollo global se presentará bajo un 
enfoque conservacionista del medio 
ambiente (ipcc, 2007). 

Cálculo del balance de humedad en el 
suelo 

El balance se realizó utilizando la 
metodología de Thornthwaite (Dunne y 
Leopold, 1978) modificada versión III, la 
cual fue adaptada en la década de los 
90's por países de la Unión Europea para 
contar con mayor información del impacto 
del cambio climático en sus regiones 
(Kenny y Harrison, 1992), así como incor- 
porada por la Convención de las 
Naciones Unidas de Lucha Contra la 
Desertificación, en países afectados por 
sequía grave y desertificación, para medir 
en todo el mundo las zonas áridas y 
semiáridas con un mismo parámetro 
(Monterroso y Gómez, 2003; Gómez et 
al., 2007 y wrs, 2004). 

Como lo señalan Gómez et al. 
(2007) el sistema incluye la aplicación del 
concepto del balance hídrico, que 
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consiste en la relación entre el agua apor- 
tada por la precipitación y la requerida por 
la vegetación en sus funciones de evapo- 
ración y transpiración (Dunne y Leopold, 
1978). El sistema de Thornthwaite consi- 
dera, para la clasificación climática, la 
eficiencia de los parámetros meteoroló- 
gicos de precipitación y temperatura. La 
evaluación del primero se determina a 
través de un balance de humedad con la 
finalidad de determinar qué tan seco o 
qué tan húmedo es el clima en relación 
con el crecimiento de las plantas. En la 
evaluación del segundo parámetro se 
estima qué tan caliente o qué tan frío es 
el clima para el crecimiento de las 
plantas. 

El sistema considera, además, la 
capacidad de almacenamiento de 
humedad del suelo, que es la cantidad 
máxima de agua disponible en las capas 
del suelo ocupadas por las raíces y que 
puede ser tomada por las mismas cuando 
la planta lo requiera para realizar sus 
funciones. En su estimación se considera 
la cantidad de humedad en el suelo a 
partir de la situación del punto de 
marchitez permanente y hasta capacidad 
de campo. 

El balance de agua por el método de 
Thornthwaite modificado no considera a 
la intercepción (int), entendida como el 
volumen de agua de lluvia que no alcanza 
a tocar el suelo por ser detenida por el 
follaje, ramas y troncos. Para evaluar este 
parámetro se siguió lo propuesto por 
Jiménez (2005) y lo observado por Flores 
et al., (2005), quienes señalan que en 
promedio los bosques tropicales inter- 
ceptan de 10 a 30% del total de la lluvia, 
mientras que para los bosques templados 
el rango se estima de un 10 a 20%, en 
promedio. Bajo este entendido se aplicó 
un valor intermedio de quince por ciento 
de intercepción de agua de lluvia para 
toda la cuenca. 


Periodos de crecimiento 

Además del balance hídrico 
mensual, se evaluó el posible impacto 
que se tendría bajo escenarios de cambio 
climático en el comportamiento del 
periodo de crecimiento (pecre), definido 
por la fao (1978) como el número de días 
durante el año en los que existe disponi- 
bilidad de agua y una temperatura favo- 
rable para el desarrollo de los cultivos. Lo 
anterior permite evaluar el periodo 
continuo del año durante el cual la 
cantidad de precipitación es mayor que la 
mitad de la evapotranspiración potencial, 
más un número determinado de días 
requeridos para evaporar una cantidad 
asumida de agua del suelo almacenada 
al final de las lluvias, que no puede ser 
mayor a 100 milímetros. 

RESULTADOS 

Variables climáticas detalladas 

En la tabla 1 se presenta el listado de 
las estaciones meteorológicas (smn, 2000) 
que se encuentran en el área de influencia 
del Parque. A partir de los valores 
mensuales de temperatura, y de acuerdo 
con Gómez et al. (2008) se calcularon 
modelos de regresión lineal simple para 
cada mes (Tabla 2), estableciendo la rela- 
ción de la variación de la temperatura a 
partir del cambio en la altitud media sobre 
el nivel del mar. Sobre el mapa de curvas 
de nivel se aplicaron los modelos esti- 
mados para obtener de manera gráfica los 
mapas de temperaturas medias 
mensuales y el anual (Figura 3a). Se 
observa que las temperaturas promedio 
anuales dentro del parque oscilan entre los 
11 y los 16°C, siendo menores en las 
partes altas del parque. El valor máximo 
de temperaturas, entre 15 y 16°C, ocupa 
una pequeña porción del área natural y se 
encuentra en las partes bajas hacia los 
límites de la zona norte de éste. 
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Tabla 1. Estaciones meteorológicas en el área de influencia del Parque Nacional El 

Chico, Hidalgo. 


CLAVE 

(SMN) 

NOMBRE 

LATITUD 

(N) 

LONGITUD 

(W) 

ALTURA 

(msnm) 

13017 

Mineral del Chico 

20,21667 

-98,73333 

2329 

13056 

Pachuca 

20,13333 

-98,75000 

2483 

13060 

Actopan 

20,26667 

-98,95000 

1991 

13061 

Alcholoya 

20,21667 

-98,45000 

2100 

13079 

Presa El Girón 

20,06440 

-98,65240 

2439 

13096 

Atotonilco 

20,28333 

-98,66667 

2080 

13098 

Huasca 

20,21667 

-98,56667 

2054 

13100 

Omitían 

20,16667 

-98,63333 

2540 

13115 

Real del Monte 

20,14588 

-98,66737 

2638 

13116 

El Sembo 

20,15000 

-98,56667 

2435 

13117 

Tlaxiaca 

20,11119 

-98,89247 

2336 

13121 

Santa María Amajac 

20,33333 

-98,70000 

1975 


Tabla 2. Modelos de regresión simple obtenidos para temperatura. 


MES 

MODELO 

Enero 

T = -0,0048h + 23,096 

Febrero 

T = -0,0051 h + 24,949 

Marzo 

T = -0,0054h + 27,692 

Abril 

T = -0,0061 h + 31,194 

Mayo 

T = -0,0072h + 34,693 

Junio 

T = -0,0074h + 34,921 

Julio 

T = -0,0068h + 32,503 

Agosto 

T = -0,0067h + 32,061 

Septiembre 

T = -0,0065h + 31,379 

Octubre 

T = -0,0062h + 29,508 

Noviembre 

T = -0,0055h + 26,434 

Diciembre 

T = -0,0049h + 23,780 

Anual 

T = -0,0060h + 28,719 
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En lo que se refiere a precipitación 
(Figura 3b), el valor medio anual de las 
áreas de influencia climática fue 881,1 mm, 
con régimen de lluvias de verano y 
presencia de canícula ligera en el mes de 
agosto. El porcentaje de precipitación 
invernal respecto a la total anual es inferior a 
5%. Un aspecto importante es que los 
vientos dominantes en la mayor parte de 
año provienen del noreste, aunque en la 
época seca son frecuentes también los que 
soplan del noroeste. Durante casi todo el 
año la zona está sujeta a la acción de 
vientos moderadamente fuertes por la tarde. 

La región de estudio corresponde a 
una zona de barlovento, donde los 
vientos provenientes del noreste se 
asocian con el ascenso forzado del aire 
sobre las laderas orientales del macizo 
montañoso y con la formación de nubes 
al enfriarse el aire por debajo del punto de 
condensación del vapor de agua. 
Asimismo, se presenta la formación de 
nubes, sobre todo cuando el gradiente de 
presión es considerable entre la zona del 


Golfo de México y el altiplano central que 
establece el flujo de los vientos alisios 
(Tinoco, 2005). 

Las zonas con mayor precipitación 
se presentan hacia el sur del parque, 
sobre el macizo montañoso que lo forma, 
observándose el máximo entre los 1 100 
y 1 200 milímetros anuales, en promedio. 
Por su parte, el rango de menor precipita- 
ción, aquel entre los 700 y 800 milímetros 
anuales, se ubica hacia las partes más 
bajas de la montaña y hacia la porción 
norte del parque. 

Aplicación de los modelos de cambio 
climático 

Los resultados de aplicar las tasas de 
cambio sugeridas por los modelos gene- 
rales de circulación (mgc) a la temperatura 
media anual del escenario base (1961- 
1990) se presentan en la tabla 3. Se 
observó un incremento, en promedio anual 
de 0,9°C según el modelo GFDL-R30 al 



Figura 3. Se muestra en a) la temperatura media anual (°C) y en b) la precipitación 

promedio anual (mm) obtenida para el parque. 



Tabla 3. Temperatura (°C) para el año base y bajo los cuatro escenarios de cambio climático evaluados. 
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año 2020 y de 1,6°C según el mismo 
modelo para el año 2050. La marcha anual 
de la temperatura bajo el mismo modelo se 
comporta de manera similar que en el 
escenario base, ya que la temperatura 
promedio más alta se presenta en el mes 
de mayo con 16,3°C para el año 2020. Sin 
embargo, para el año 2050 el mes de abril 
será el más caluroso con 17,1°C de 
acuerdo con el modelo ya señalado. Por 
su parte, bajo el modelo HadCM3 se 
presentará un aumento, en promedio, en 
la temperatura media anual de 1 °C para el 
año 2020 y de 2,4°C para el 2050. El mes 
de mayo se mantendrá como el más calu- 
roso con 16,6°C y 17,6°C para los años 
2020 y 2050. 

En cuanto a la precipitación media 
anual (Tabla 4), se observa un aumento 
en los dos periodos bajo el modelo norte- 
americano, 14% para el año 2020 y 
17,7% para el año 2050. Bajo condi- 
ciones actuales, el mes con la mayor 
cantidad de lluvia es septiembre (170 
mm) pero bajo el modelo GFDL-R30, el 
mes de julio será en donde se presente la 
mayor precipitación pluvial con 192,1 mm 
para el año 2020 y 210,8 mm hacia el año 
2050. El modelo inglés, por el contrario, 
muestra una disminución en los valores 
de precipitación media anual en compara- 
ción con las condiciones actuales, ya que 
para el año 2020 supone que disminuye 
en 8,5% y para el 2050 en 5,6%. Bajo 
este escenario, el mes de junio será el 
que registre la mayor cantidad de lámina 
precipitada con 142,5 mm en el 2020 y 
155,3 mm para el 2050. 

Balance de humedad en el parque 

Los datos obtenidos anteriormente, 
junto con los valores de latitud y los tipos 
de vegetación, sirvieron para obtener las 
variables del balance de humedad en el 
suelo: evapotranspiración (et), escurri- 
miento (esc), infiltración (inf) e intercep- 


ción (int). Los resultados del balance 
hídrico general se presentan en la tabla 5. 

El comportamiento de la evapotrans- 
piración anual es a la alza bajo los dos 
modelos para ambos periodos de tiempo 
evaluados, ya que para el modelo GFDL- 
R30 se presenta un aumento de 3,6% en 
el año 2020 y de 6,2% para el año 2050. El 
modelo HadCM3 establece un incremento 
mayor. Para el periodo 2020 6,2%. Lo 
anterior se puede explicar debido al 
aumento de la temperatura sugerido por 
ambos modelos de cambio climático. 

El escurrimiento por su parte 
presenta, al igual que la precipitación, 
cambios dependientes del modelo apli- 
cado: con el norteamericano registra un 
aumento de 30,8% para el año 2020 y de 
37,9% para el 2050. Por su parte, el 
modelo inglés registra una disminución de 
25,6% en el volumen escurrido para el 
año 2020 y de 34,5% para el año 2050. 

En lo que se refiere a la cantidad de 
infiltración anual, bajo el escenario base 
representa 6% del total de la lluvia obser- 
vada para el periodo, es decir, un volumen 
de 1,48x1 0 6 m 3 de agua (Figura 4). Los 
modelos de cambio climático sugieren que 
para el año 2020 estos valores disminuirán 
en 70% , según el modelo inglés y 15% 
según el modelo norteamericano, obser- 
vando volúmenes de 0,4x1 0 6 m 3 y 
1,26x1 0 6 m 3 de agua, respectivamente. 
Para los escenarios al año 2050 la ten- 
dencia señalada por el modelo HadCM3 
prevalece, al disminuir en 99%, lo que signi- 
fica una pérdida de un volumen de 1 ,4x1 0 6 
m 3 de agua, para quedar en tan solo un 
volumen de 0,008x1 0 6 m 3 El modelo 
GFDL-R30 para el año 2050 es el único que 
supone un incremento en la cantidad de 
volumen infiltrado de aproximadamente 
30% más de agua. Entendiéndose lo ante- 
rior porque el mismo modelo sugiere mayor 
cantidad de lluvia precipitada para dicho 
periodo. 
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Tabla 5. Lámina de lluvia promedio anual y de las variables del balance de humedad 

para los cinco escenarios evaluados (mm/año). 


ESCENARIO 

AÑO 

P 

INT 

ET 

ESC 

INF 

BASE 

1961-1990 

881,1 

132,1 

432,2 

264,6 

52,2 

HadCM3_A2 

2020 

806,5 

120,9 

474,3 

195,1 

16,5 

GFDL-R30_A2 


1004,2 

150,6 

463,6 

346,1 

43,9 

HadCM3_A2 

2050 

832,0 

124,8 

533,7 

173,4 

0,1 

GFDL-R30 A2 


1037,4 

155,5 

446,7 

365,0 

70,2 


Donde P=Precipitación, INT=lntercepción, ET=Evapotranspiración, ESC=Escurrimiento e 
INF=lnfiltración 


30 



5 


o 



Figura 4. Volumen de agua (10 6 m 3 ) del balance hídrico para el escenario base y cuatro 
escenarios de cambio climático evaluados en el Parque Nacional El Chico, Hidalgo. 
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Analizando la distribución mensual 
de las variables promedio de la zona 
(Figura 5a), para el año base se observa 
que en la precipitación existe una dismi- 
nución de aproximadamente 20 milíme- 
tros en el mes de agosto respecto de los 
meses anteriores, junio y julio Para el 
mes siguiente, es decir, septiembre, se 
presenta el máximo anual de precipita- 
ción con más de 160 milímetros de lluvia 
en el mes y aproximadamente 30 milíme- 
tros más que en agosto. Este fenómeno 
ha sido ampliamente descrito como un 
periodo de sequía intraestival o canícula. 
Para la región de estudio, se puede 
resumir el comportamiento de la precipita- 
ción señalando que el periodo de lluvias 
es de régimen de verano y el mes más 
húmedo presenta más de 10 veces la 
precipitación del mes más seco y ocurre 
en la mitad que es la más caliente del año 
y para este caso particular después de la 
canícula. Al aplicar los escenarios de 
cambio climático se observó un posible 
cambio en este comportamiento descrito. 

Los cuatro escenarios de cambio 
climático aplicados (Figura 5b, 5c, 5d y 
5e) sugieren que la canícula se seguirá 
presentando en la región de estudio en el 
mes de agosto. Sin embargo, los modelos 
apuntan a que el mes más húmedo ya no 
se presentará después de la canícula, 
sino antes de ésta, en el mes de junio (3 
escenarios) o julio (1 escenario). Para el 
primer horizonte de tiempo al año 2020, 
los dos modelos sugieren un comporta- 
miento muy similar de la canícula al 
observado en el año base. 

Para los modelos aplicados al año 
2050, el comportamiento de la canícula 
se anuncia más drástico. Para el modelo 
inglés representa una disminución de 
aproximadamente 60 mm en compara- 
ción con los meses que la anteceden 
(junio y julio), observándose un ligero 
incremento para el mes de septiembre de 
apenas 20 mm (Figura 5d). Éste repre- 


sentará el escenario más radical para el 
balance hídrico, ya que el máximo de 
precipitación se presentará en los meses 
de junio y julio. El mes de agosto será de 
canícula, pero ésta no alcanzará los 
niveles de capacidad de infiltración obser- 
vados en el año base. Por su parte, y para 
el mismo año, el modelo norteamericano 
(Figura 5e) sugiere la presencia de la 
canícula en el mes de agosto, con el mes 
más lluvioso en el mes de julio y con una 
posible elevación considerable de las 
lluvias en el mes de septiembre. Sin 
embargo, la canícula prevista presentará 
una disminución de aproximadamente 90 
mm en comparación con el mes de julio y 
de 80 mm en comparación con el mes de 
septiembre, impactando también en el 
balance hídrico y la capacidad de infiltrar 
agua a los acuíferos, ya que las disminu- 
ciones observadas en el año base fueron 
de 20 y 30 mm, respectivamente. 

En lo que se refiere a la evapotrans- 
piración, el máximo observado en el 
escenario base se presenta en el mes de 
mayo con poco más de 77 mm, mientras 
que el mínimo se encuentra en febrero 
con aproximadamente 35 mm. Este 
comportamiento permanece para todos 
los escenarios de cambio climático apli- 
cados, debiéndose a que después de la 
época de lluvias queda agua almacenada 
en el suelo durante la estación de 
invierno, además de la disminución en 
temperatura presentada en el mismo 
periodo, repercutiendo en valores infe- 
riores de evapotranspiración. Conforme 
se acerca la época lluviosa y aumentan 
las temperaturas, la evapotranspiración 
también aumenta hasta el máximo en el 
mes de mayo. 

El escurrimiento, por su parte, se 
presenta con valores máximos en el mes 
de septiembre con aproximadamente 80 
mm, coincidiendo con el máximo de preci- 
pitación en el escenario base. Poco 
después de que inicia la época de lluvias 
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y al saturarse el suelo se inicia el proceso 
de escurrimiento, comprendiendo poco 
más de seis meses, desde el mes de 
junio y hasta diciembre, con valores muy 
pobres en el mes de enero. En los meses 
de febrero a junio no se aprecia escurri- 
miento. Al aplicar los escenarios de 
cambio climático la tendencia es muy 
parecida al escenario base, observán- 
dose que el mes de junio podría presentar 
un poco más de escurrimientos que los 
observados para el periodo base. El mes 
de septiembre continuará siendo el de 
mayor escurrimiento, dado que se ha 
cubierto la capacidad de almacenamiento 
de agua del suelo con la precipitación de 
junio, julio y agosto, aunque para los dos 
modelos del escenario 2020 se espera 
disminuyan a menos de 60 mm y según el 
modelo inglés para el 2050, a menos de 
40 mm. El modelo norteamericano para el 
año 2050 sugiere que los valores sobre- 
pasarán los 100 mm escurridos, debido a 
que se esperarán mayores lluvias 
después de la canícula en el mes de 
agosto. 

En cuanto a la infiltración de agua, 
que es el objetivo central del presente 
estudio, el escenario base indica que se 
tienen ocho meses del año (octubre a 
mayo) con infiltración nula. Los restantes 
cuatro meses, dentro del periodo de 
lluvias, son donde se presenta la posibi- 
lidad de infiltrar agua. El mes de junio es 
donde se presenta la mayor cantidad de 
agua infiltrada con poco más de 50 mm, 
julio con aproximadamente 30 mm y los 
meses de agosto y septiembre con poco 
más de 10 mm en conjunto. Los valores 
en volumen de agua que esto representa 
ya fueron señalados en párrafos ante- 
riores. 

Para los escenarios futuros con 
cambio climático al año 2020 el compor- 
tamiento es similar, aunque el mes de 
agosto, donde se presenta la canícula, 


disminuirá considerablemente la capa- 
cidad de infiltrar agua. El mes de junio 
seguirá siendo el de mayor infiltración, 
pero los valores se reducirán a poco 
menos de 40 mm, lo que representará 
una disminución del 35% de agua infil- 
trada. 

Los escenarios de cambio para el 
año 2050 señalan comportamientos 
distintos entre sí y comparados con el 
escenario base. El modelo inglés señala 
que únicamente en dos meses del año se 
tendrá capacidad de infiltrar agua, el mes 
de junio con cerca de 50 mm y el mes de 
julio con poco más de 15 mm. Para el 
primero representa una disminución del 
20% comparado con el mismo mes en el 
escenario base y para el segundo una 
disminución del 50%. De acuerdo con el 
modelo norteamericano, para el mismo 
año el periodo con capacidad de infiltrar 
agua aumentará, empezando en el mes 
de mayo, disminuyendo considerable- 
mente en agosto y terminando en 
septiembre, posiblemente aumentando el 
volumen de agua infiltrada. 


Periodos de crecimiento 

Para el escenario base se encontró 
que el periodo de crecimiento inicia el 8 
de abril y termina el 1 de diciembre, 
teniendo una duración de siete meses y 
23 días. Los escenarios de cambio climá- 
tico no sugieren modificaciones impor- 
tantes en el pecre. Al contrario, todos 
simulan un ligero aumento de este 
periodo. Para el año 2020, los modelos 
sugieren un aumento de cuatro y 
dieciocho días, de acuerdo con el modelo 
inglés y norteamericano, respectiva- 
mente. Por su parte, para el año 2050 los 
modelos apuntan a aumentar dieciséis y 
trece días, en comparación con el esce- 
nario base para el modelo inglés y el 
norteamericano, respectivamente. 
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Figura 5. Balance hídrico mensual para a) Año Base, b) Modelo HadCM al año 2020, 
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Además, como es de esperarse, al 
aumentar el pecre, los escenarios de 
cambio climático señalan que podrían 
existir cambios en la fecha de inicio y fin 
del periodo. Para el año 2020 el modelo 
inglés sugiere que el pecre podría retra- 
sarse hasta cuatro días, al encontrarse 
que podría iniciar el periodo de creci- 
miento el 12 de abril. Los restantes tres 
escenarios de cambio climático sugieren 
que podría adelantarse el inicio del 
periodo de crecimiento en 7 días (modelo 
inglés al 2050) o en trece días (ambos 
escenarios del modelo norteamericano). 

DISCUSIÓN Y CONCLUSIONES 
Caracterización climática detallada 

El parque presenta un clima, de 
acuerdo con García (2004), del tipo 
Cb(w 2 )(i')gw”, tomando como referencia 
los datos promedios del parque, que 
corresponde a un clima templado, con 
verano fresco y largo, el más húmedo de 
los subhúmedos, con régimen de lluvias 
en verano y un porcentaje de precipita- 
ción invernal entre 5 y 10,2 mm, con poca 
oscilación anual de la temperatura, 
marcha anual de la temperatura del tipo 
Ganges y presencia de canícula. El 
programa de manejo (Meló y López, 
2003, citado en conanp, 2006) también 
señala valores similares de temperatura, 
observándose que para la presente inves- 
tigación se obtuvieron valores más deta- 
llados para esta variable anual, variando 
de los 11 a los 16°C. Lo anterior se puede 
deber a que para dicha descripción se 
utilizó cartografía de escala más 
pequeña. 

Los modelos de regresión lineal 
obtenidos para temperatura presentan 
variaciones en el rango de incremento de 
la temperatura conforme se aumenta la 
altura. Al ser una zona de barlovento, se 
observan variaciones altotérmicas 


menores que en una zona de sotavento 
en el mismo estado, esto se asocia al 
mayor contenido de vapor de agua del 
aire en esta zona, que al ascender sobre 
la pendiente se condensa y se libera calor 
latente al pasar el agua de vapor a 
líquido, mismo que calienta el aire y 
reduce el rango de enfriamiento (Arhens, 
2003). Al aplicar lo modelos de regresión 
obtenidos sobre las curvas de nivel para 
la zona de estudio, se obtuvieron los 
mapas detallados de temperaturas 
mensuales y anual, observándose así el 
mejor detalle en la caracterización de 
esta variable. 

Por otra parte, el trazo de ¡soyetas 
mensuales muestra que la distribución de 
la precipitación en la zona es primordial- 
mente en verano, como resultado del 
efecto de los vientos alisios provenientes 
del Golfo de México. Así como con los 
huracanes que se forman en el Caribe y 
siguen su curso por el Golfo entre los 
meses de julio a octubre. Durante la 
época invernal las lluvias son influen- 
ciadas por los frentes fríos llamados 
“Nortes”, mismos que están asociados al 
componente invernal del Monzón de 
Norteamérica, los cuales también aportan 
lluvia, aunque en menor cantidad que los 
sistemas de verano. 

Modelos de cambio climático 

Los modelos de cambio climático 
concuerdan, en cuanto a que para la zona 
de estudio se incrementará la tempera- 
tura inequívocamente. Dependiendo del 
modelo aplicado se obtuvieron incre- 
mentos que pueden ir desde 1°C y hasta 
los 2,4°C, en promedio para la zona. Este 
incremento de temperatura está amplia- 
mente explicado por la acumulación de 
gases de efecto invernadero, que al ir 
aumentando su concentración en la 
atmósfera causan una mayor retrasmi- 
sión a la superficie de la radiación de 
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onda larga irradiada por la superficie de la 
tierra, lo que calienta más el aire y 
provoca un incremento en la temperatura. 
Para el escenario más tardío se prevé 
mayor acumulación de estos gases y la 
magnificación del fenómeno de calenta- 
miento (ipcc, 2001). 

La razón por la que los modelos de 
cambio climático difieren se asocia con 
las series de tiempo con las que se cali- 
braron dichos modelos. El modelo norte- 
americano considera una cantidad signi- 
ficativa de estaciones de zonas con altas 
concentraciones de aerosoles, que son 
comunes en la costa este de dicho país, 
lo que reduce el calentamiento. Por su 
parte, el modelo inglés utiliza datos más 
globales donde estos elementos tienen 
menor peso (McGregor y Nieuwolt, 
1998). 

Las diferencias en precipitación esti- 
madas en los escenarios de los modelos, 
en donde el modelo inglés reporta una 
disminución de aproximadamente 80 y 50 
mm (-10 y -5%) para los años 2020 y 
2050, respectivamente, se explica por un 
desplazamiento pequeño hacia el sur de 
la zona de corrientes descendentes y el 
menor avance de los centros de alta 
presión hacia el océano en la temporada 
caliente del año, lo que disminuirá la 
precipitación en las zonas aledañas al 
movimiento de estos sistemas. En 
cambio, el modelo norteamericano 
pronostica una mayor evaporación y con 
esto un mayor contenido de vapor de 
agua en la atmosfera, sin que esta zona 
se vea afectada considerablemente por el 
movimiento de los sistemas de vientos 
anteriores, lo que trae consigo un incre- 
mento en la precipitación (McGregor y 
Nieuwolt, 1998) de aproximadamente 200 
mm para ambos años evaluados con el 
modelo (15% más). 


Balance hídrico actual 

Del total de la precipitación obser- 
vada, en promedio para el periodo 1961- 
1990, 49% es evapotranspirada, 30% se 
escurre, 15% es interceptada y 6% se 
infiltra. Cabe señalar que estos valores 
son en promedio para el parque, que 
aunque se trabajó a más detalle geográ- 
fico, el balance permite establecer 
patrones generales de comportamiento 
para la zona. Sin lugar a dudas hacia las 
partes más altas del parque estos valores 
serán diferentes a las zonas planas y de 
menor altitud. La cobertura de vegetación 
predominante son los bosques de 
oyamel, pino-encino y oyamel-encino, 
quienes en conjunto ocupan alrededor de 
2 300 ha que significa poco más del 85% 
del total del parque. Debido a lo anterior, 
se puede entender el porcentaje que 
representa la evapotranspiración con 
relación del total precipitado, ya que el 
tipo de vegetación es altamente deman- 
dante de agua para cubrir sus funciones. 

El porcentaje de intercepción del 
agua de lluvia por la vegetación fue un 
parámetro evaluado de acuerdo a lo 
sugerido por la literatura, donde se consi- 
deró como un valor intermedio de lo 
observado por otros autores, por lo que 
se reconoce que para obtener resultados 
más reales debería de apoyarse con 
estudios locales para la variable. Final- 
mente, el valor de infiltración deberá ser 
considerado como teórico, ya que 
algunas características de los suelos en 
la región de estudio no fueron conside- 
radas en el balance. Sin embargo, el valor 
obtenido no difiere considerablemente a 
lo señalado por cna para los acuíferos 
estudiados. 

Con lo anterior, se considera que el 
balance hídrico aquí presentado refleja el 
comportamiento general del agua en la 
región de estudio, por lo que al aplicar 
escenarios de cambio climático se podrán 
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considerar como escenarios futuros que 
son viables de experimentarse en el 
parque. 

Balance hídrico bajo escenarios de 
cambio climático 

El comportamiento de las variables 
climáticas, bajo escenarios de cambio 
climático, fluctúa dependiendo del modelo 
aplicado. En ambos modelos la tempera- 
tura se incrementará y por lo tanto la 
evapotranspiración también lo hará. Lo 
anterior se puede traducir en déficit 
hídrico en cuanto a humedad retenida e 
infiltrada se refiere. 

En el caso del modelo inglés, que 
sugiere que la precipitación pluvial dismi- 
nuirá, se observa una disminución en la 
cantidad de agua que escurre de 26% y 
35% y de la cantidad de agua infiltrada de 
26% y 32% para los escenarios al 2020 y 
2050, respectivamente. Por su parte, en 
el modelo norteamericano, aunque 
sugiere aumentos en la evapotranspira- 
ción, se esperan incrementos tanto del 
escurrimiento (de 31% y 38%) como en la 
infiltración (de 6% y 23%), para los esce- 
narios de 2020 y 2050, respectivamente, 
esto como consecuencia del incremento 
en las lluvias pronosticadas. 

Este comportamiento es de espe- 
rarse debido a que el modelo norteameri- 
cano indica que el mes más húmedo será 
el mes de julio, en comparación con el 
promedio de los años 1961-1990, donde 
se ha observado que es el mes de 
septiembre. Esta situación corresponderá 
a que de presentarse más precipitación 
en los meses próximos al final de la esta- 
ción seca (como lo sugiere el modelo 
inglés) el suelo tendrá mayor capacidad 
de retener humedad. En cambio, de 
presentarse enmedio de la temporada 
húmeda (como lo sugiere el modelo 
norteamericano) el suelo ya se encon- 


trará húmedo, por lo que un porcentaje 
mayor de la precipitación se escurrirá al 
rebasarse la capacidad de infiltración de 
los suelos. 

Como se presenta en el programa 
de conservación y manejo (conanp, 
2006), el Parque Nacional El Chico repre- 
senta una importante zona de recarga de 
agua al acuífero. De acuerdo con cna 
(2007), el parque es una zona de recarga 
de los acuíferos Actopan-Santiago y 
Cuautitlán-Pachuca. Para el primero, se 
señala que la recarga natural es de 171,9 
hm3/año, mientras que la extracción es 
del orden de 40,5 hm3/año, observando 
una condición de subexplotación. Por su 
parte, para el acuífero Cuautitlán- 
Pachuca, las cifras señalan una recarga 
de 204 hm3/año y una extracción de 
483,3 hm3/año, lo que conlleva a una 
condición de sobreexplotación (cna, 
2007; dof, 31 enero de 2003). 

En caso de presentarse condiciones 
como las establecidas bajo el modelo 
norteamericano, el abastecimiento del 
acuífero no se vería tan afectado, ya que 
la infiltración y el escurrimiento posible- 
mente aumentarán. Sin embargo, como 
ya se presentó, existe la posibilidad de 
que se registren eventos extremos de 
lluvia, lo que pondría en riesgo a las 
poblaciones que se encuentren cerca de 
los cauces de los ríos o pendiente abajo 
de las presas o lagos, ya que habría la 
incidencia de inundaciones o desborda- 
mientos. Además, el comportamiento de 
este modelo esconde un marcado 
aumento en la canícula, que aunque se 
prevé mayor cantidad de lluvia, la dismi- 
nución sugerida en el mes de agosto 
podrá ser considerable e impactará en el 
comportamiento de las especies animales 
y vegetales. 

Por su parte, bajo las condiciones 
del modelo inglés, que subraya una 
posible disminución en la cantidad de 
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agua infiltrada y escurrida, se presenta- 
rían condiciones de mayor estrés por 
déficit de agua en las comunidades vege- 
tales del parque. Al disminuir la recarga 
de los acuíferos, como consecuencia de 
lo antes señalado, y si la demanda de 
agua sigue en aumento, éstos empeza- 
rían a sufrir un abatimiento, lo que impli- 
caría desabasto de agua potable en las 
comunidades que dependen de pozos y 
manantiales. Por consiguiente, también 
se tendrían impactos negativos en la 
producción agrícola y pecuaria. Este 
modelo simula escenarios más extremos, 
lo que para fines de planeación de 
medidas de adaptación y mitigación debe 
de ser considerado sobre el modelo 
norteamericano para tener las mejores 
propuestas y estrategias a implementar y 
que reduzcan efectivamente los impactos 
negativos del cambio climático en el 
balance de humedad de los acuíferos a 
los que la superficie del parque abastece. 

Una de las ventajas de elaborar 
balances hídricos mensuales es que se 
puede observar el comportamiento de las 
variables estudiadas bajo las condiciones 
actuales y con más detalle para los esce- 
narios de cambio climático. Por ejemplo, 
los cuatro escenarios de cambio climático 
aplicados sugieren que la canícula se 
seguirá presentando en la región de 
estudio en el mes de agosto. Sin 
embargo, los modelos apuntan a que el 
mes más húmedo ya no se presentará 
después de la canícula, sino antes de 
ésta, en el mes de junio o julio. Este 
comportamiento no se ve reflejado en el 
balance hídrico anual. 

La conducta de la canícula prevista 
por los modelos sugiere incongruencias, 
ya que de acuerdo con la circulación 
general de los vientos se torna difícil 
entender porqué se podría presentar en 
el mes de junio el máximo de precipita- 
ción. Sin embargo, de acuerdo al modelo 
norteamericano, para el año 2050, el 


máximo de precipitación podría presen- 
tarse en el mes de julio, lo que pareciera 
ser más lógico con dicho comportamiento 
de los vientos. Este comportamiento 
deberá ser mejor estudiado en futuros 
estudios. 


Servicio ambiental hidrológico 

La realización de un balance hídrico 
mensual bajo condiciones actuales y con 
escenarios de cambio climático permitió 
conocer el comportamiento hidrológico 
actual y futuro posible en el área de 
estudio. El enfoque que se siguió fue el 
de evaluar la capacidad que tienen los 
ecosistemas en el parque de regular y 
mantener dicho balance hidrológico, para 
contar con mayores elementos que 
sustenten el análisis del servicio 
ambiental estudiado. De tal manera que 
se reconoce el valor que tienen las varia- 
bles ambientales físicas (por ejemplo 
suelo y pendientes) como las biológicas 
(particularmente la vegetación) en la 
capacidad de los ecosistemas de proveer 
servicios ambientales. 

El Parque Nacional El Chico es 
importante debido a la cantidad de servi- 
cios ambientales que provee a las comu- 
nidades que habitan a su alrededor, como 
es la capacidad de retener e infiltrar agua. 
Además, no solo se debe considerar la 
recarga de agua hacia los acuíferos, sino 
también otros servicios como la produc- 
ción de oxígeno, captura de carbono y 
recreación, entre otros. 

Sin embargo, es prioritario consi- 
derar el abastecimiento de agua hacia el 
subsuelo debido a la gran cantidad de 
actividades que dependen de ella. La 
descripción técnica y numérica del 
balance hídrico mensual permitió reco- 
nocer el papel que juegan otras variables 
ambientales de los ecosistemas (suelos y 
vegetación, por ejemplo), ya que, si bien 
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no fueron evaluadas directamente, si 
permiten reconocer el servicio ambiental 
que suministran a la sociedad y de 
acuerdo al presente estudio, a la regula- 
ción hídrica. 

Es por lo anterior que se recomienda 
proteger y aumentar las zonas mejor 
conservadas del parque, como aquéllas 
con vegetación natural de bosque, para 
evitar que, como lo señalan los modelos, 
disminuya la cantidad de agua disponible 
para el uso humano. De tal manera que el 
servicio ambiental hidrológico, de captura 
y provisión de agua, no se vea alterado 
en sus diferentes componentes. 

Finalmente consideramos que los 
futuros estudios en la región deberán 
encaminarse a la evaluación de otros 
servicios ambientales, ya que el presente 
documento trata sólo de uno de ellos, 
pero para comprender el valor del 
conjunto de servicios es menester la 
evaluación total de los servicios ambien- 
tales proporcionados por los ecosistemas 
en el parque. 
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Reconstrucción de 350 a os de precipitación 
para el suroeste de Chihuahua, México 

Reconstructng 350 years of precipitation in southeast 

Chihuahua, México 

Julián Cerano Paredesi, José Villanueva Díazi, 

PeterZ. Fulé 2 , Jesús Guadalupe Arreóla Avilas, 

Ignacio Sánchez Cohén 1 y Ricardo David Valdez Cepeda 3 

RESUMEN 

En el suroeste del estado de Chihuahua, México, en la reserva “Cerro El Mohinora” se desarro- 
llaron series de tiempo dendrocronológicas de madera temprana, tardía y anillo total de Pseudotsuga 
menziesii con una longitud de tres siglos y medio (1657-2005), con el objetivo de analizar la variabi- 
lidad hidroclimática histórica de la región. La reconstrucción de precipitación invierno-primavera indica 
gran variabilidad interanual, decenal y multidecenal de los patrones de precipitación para la región. 
Sequías severas se reconstruyeron para los periodos 1695-1715, 1753-1760, 1785-1792, 1798-1806, 
1819-1830, 1841-1870, 1890-1897, 1906-1912, 1924-1941, 1971-1977 y 1994-2005, aunque las 
sequías más prolongadas ocurrieron en los periodos 1695-1715, 1841-1870 y 1924-1941. Las últimas 
tres décadas del siglo XX y los años corrientes de la primera década del siglo XXI (1971-2005) indican 
un periodo de intensa sequía para la región, con impactos en lo ecológico y socioeconómico aun no 
cuantificados. 

PALABRAS CLAVE: 

Madera tardía, madera temprana, Pseudotsuga menziesii , reconstrucción de precipitación. 


ABSTRACT 

Earlywood, latewood, and total ring width series from Pseudotsuga menziesii were developed in a 
mixed conifer forest at Cerro El Mohinora”, Chihuahua. The tree-ring series extended for the period 
1657-2005 (349 years) and were used to analyze hydroclimatic variability for southwestern Chihuahua. 
A seasonal winter-spring precipitation reconstruction was developed with the objective to determine high 
and low frequency precipitation variability and to analyze the impact of dominant atmospheric circulatory 
patterns. The reconstruction indicated high interannual, decadal and multiannual variability. Severe 
droughts were reconstructed for the periods 1695-1715, 1753-1760, 1785-1792, 1798-1806, 1819-1830, 
1841-1870, 1890-1897, 1906-1912, 1924-1941, 1971-1977, and 1994-2005. The most extended 
droughts took place in the periods 1695-1715, 1841-1870, and 1924-1941 . The last three decades of the 
20 th century and the current years of the 21 st century have been particularly dry for northwestern 
Chihuahua but the ecological and socieoconomical impacts have not been determined. 
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Latewood, earlywood, Pseudotsuga menziesii , precipitation reconstruction. 
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INTRODUCCIÓN 

Hoy en día, el entendimiento de la 
variabilidad del clima es uno de los temas 
más relevantes de las ciencias ambien- 
tales. El clima varía de un año a otro, 
entre décadas e incluso entre milenios. 
La actividad industrial y los gases de 
efecto invernadero liberados a la atmós- 
fera a consecuencia de la actividad antro- 
pogénica, han alterado el clima de la 
tierra (Houghton et al., 1990). En este 
sentido, se han producido cambios en los 
patrones de precipitación asociados al 
calentamiento global (Bradley etal., 1987, 
Díaz et al., 1989), alterando el comporta- 
miento natural del ciclo hidrológico (ipcc, 
1996). Aunque no existen evidencias 
claras a escalas locales sobre los efectos 
del calentamiento global de la tierra, 
existe consenso científico en cuanto a 
que hay una gran incertidumbre con 
respecto a la magnitud y tasas de estos 
cambios a escalas regionales (eei, 1997). 

Los registros del clima procedentes 
de satélites y registros instrumentales 
cubren por lo general menos de 80 años. 
Estos registros son demasiado cortos 
para examinar el rango completo de la 
variabilidad climática. Existen diferentes 
fuentes naturales indirectas que registran 
dicha variación en el tiempo, tales como: 
anillos de árboles, capas de hielo, polen 
fósil, sedimentos marinos, corales y docu- 
mentos históricos. Al analizar los registros 
tomados de éstas y de otras fuentes paleo- 
climáticas, se puede extender el conoci- 
miento de la variación del clima mucho 
más allá de los registros climáticos instru- 
mentales (Bradley, 1999). 

Los anillos de árboles constituyen la 
fuente indirecta más apropiada de resolu- 
ción anual para entender el comporta- 
miento histórico del clima, su crecimiento 
anual está influenciado por las condi- 
ciones climáticas dominantes previas o 
durante su formación. Son una herra- 


mienta para el análisis de las condiciones 
climáticas en periodos en los que no 
existen registros instrumentales (Fritts, 
1976). 

Los bosques templados de la Sierra 
Madre Occidental poseen coniferas con 
excelente potencial dendrocronológico, 
en el que el desarrollo tanto de madera 
temprana como tardía está influenciado 
por las condiciones dominantes de preci- 
pitación de invierno y de verano, que a su 
vez son modulados por patrones atmosfé- 
ricos de circulación global (Stahle et al., 
1998). 

Los estudios de variabilidad hidrocli- 
mática han sido poco desarrollados en la 
región de Chihuahua (Díaz et al., 2002; 
Villanueva-Díaz et al., 2006). La reserva 
Cerro El Mohinora, al suroeste del estado 
de Chihuahua, presenta ecosistemas de 
bosques mixtos con especies de exce- 
lente potencial dendroclimático como es 
el caso de Pseudotsuga menziesii, repre- 
sentando estos sitios un gran potencial 
para reconstruir la precipitación, informa- 
ción que ayudará a un mejor entendi- 
miento de la variabilidad hidroclimática a 
nivel regional y a entender su impacto en 
los ecosistemas. 

OBJETIVOS 

1) Generar una serie dendrocronoló- 
gica con base en los anillos de creci- 
miento de Pseudotsuga menziesii. 

2) Reconstruir la precipitación de los 
últimos tres siglos y medio y analizar su 
variabilidad hidroclimática en el tiempo. 

METODOLOGÍA 
Localización del área 

El sitio Cerro El Mohinora se localiza 
al suroeste del estado de Chihuahua, 



Madera y Bosques 15(2), 2009:27-44 

cerca de los límites con los estados de 
Durango y Sinaloa (Figura 1). El área 
comprende un gradiente de elevación de 
2 500 a 3 000 m, donde se desarrolla un 
bosque de pino y Abies-Pseudotsuga- 
Picea en la parte alta (25° 57' 16" N, 
107° 01' 38" W). En la parte media se 
localiza un bosque mixto de Abies-Pseu- 
dotsuga-Picea (25° 58' 27" N, 107° 02' 
26" W) y en la parte baja bosque de pino 
(25° 58' 23" N, 107° 01' 35" W). 

El área presenta temperaturas 
promedio máximas de 25 a 27°C para los 
meses de mayo a septiembre; tempera- 
turas promedio mínimas de -1 a -3°C para 
los meses de diciembre a marzo. Registra 
una precipitación promedio para los 
meses más lluviosos del año de 110 mm 
en junio, 274 mm en julio, 249 mm en 
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agosto, 164 mm para septiembre y un 
promedio anual de 1 130 mm. Los prome- 
dios de las variables se generaron de un 
periodo de 25 años (eric ii, 2000). 

Para reconstruir la variabilidad 
climática de la región suroeste del estado 
de Chihuahua, se tomaron núcleos de 
crecimiento y secciones transversales de 
Pseudotsuga menziesii. Se empleó un 
muestreo selectivo; en cada árbol selec- 
cionado, a la altura del pecho y perpendi- 
cular a la pendiente, con taladro Pressler 
de 5 mm de diámetro interno y exten- 
siones variables, se extrajeron dos virutas 
(núcleos de crecimiento) por árbol (Figura 
2); asimismo, con motosierra se cortaron 
secciones transversales de tocones, 
árboles muertos o madera subfósil 
(secciones de madera sobre la superficie 
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Figura 1. Localización geográfica del sitio de muestreo en la reserva “Cerro El 
Mohinora” al suroeste del estado de Chihuahua. 
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Figura 2. Toma de muestras de árboles vivos (virutas o núcleos de crecimiento) 
empleando el Taladro de Pressler, a) muestra la posición de la barrena para la toma de 
la muestra y b) la forma en que se extrae el núcleo de crecimiento. 


del suelo o semienterrada). Los núcleos 
de crecimiento se identificaron, montaron 
y pegaron en secciones acanaladas de 
madera para facilitar su manejo. 

Las muestras se pulieron con 
diversos grados de lija (del número 30 la 
más gruesa, a la más fina del número 
1200) para resaltar sus estructuras de 
crecimiento. En el Laboratorio de Dendro- 
cronología del inifap ceñid-raspa se 
contaron los anillos con el apoyo de un 
estereomicroscopio de alta resolución, los 
anillos se contaron y fecharon al año 
exacto de su formación con base en 
técnicas dendrocronológicas estándar 
(Stokes y Smiley, 1968). 

El principio del fechado cruzado es 
el principio básico de la dendrocrono- 
logía. Es la herramienta principal por 
medio de la cual es determinado el creci- 
miento exacto de cada uno de los anillos 
anuales. Un simple conteo de anillos 


posiblemente generaría errores debido a 
que existen anillos ausentes (perdidos) y 
falsos anillos (dobles). El fechado 
cruzado es imprescindible cuando las 
mediciones del ancho de los anillos son 
comparadas con alguna variable climática 
anual o estacional, tal como precipitación, 
temperatura, presión atmosférica, etc. Sin 
un fechado preciso de los anillos de creci- 
miento, es imposible una calibración 
adecuada debido a que la cronología 
estará desfasada en el tiempo. Esta 
técnica permite observar en una muestra, 
donde existen anillos perdidos o donde 
podríamos tener formados dos o más 
anillos en un año. El fechado al año de 
formación de cada uno de los creci- 
mientos anuales es el objetivo principal 
de esta técnica (Figura 3). 

Una vez fechadas las muestras, se 
midieron individualmente las bandas de 
madera temprana, tardía y anillo total con 
un sistema de medición velmex, equipo 
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Figura 3. Fechado Cruzado, Principio de la Dendrocronología, que permite determinar 
el año exacto de formación de cada anillo de crecimiento y que con base en la similitud 
de los patrones de crecimiento es factible extender en el tiempo las series 
dendrocronológicas, al lograr datar árboles vivos, troncos muertos, enterrados o 
semienterrados, que a su vez, determinarán la extensión en el tiempo para propósitos 

de reconstrucción climática (Fritts, 1976). 


compuesto por una platina deslizable y un 
codificador linear que tiene una precisión 
de medición de 0.001 mm, sistema inte- 
grado a un procesador pc y un estereomi- 
croscopio con una retícula montada en el 
ocular para ubicar los límites de cada 
anillo de crecimiento. 

Después de generada la base de 
datos, se verificó la calidad del fechado, 
la respuesta climática y la exactitud de la 
medición de cada anillo con el programa 
cofecha (Holmes, 1983). De las series de 
crecimiento, las tendencias biológicas 
como competencia y liberación y las 
geométricas (el área de fuste se incre- 
menta con la edad y el crecimiento anual 
tiende a disminuir al distribuirse en una 
mayor superficie) no relacionadas con 
clima, se removieron con el programa 


arstan, al insertar una curva exponencial 
negativa o línea recta a la serie de medi- 
ción y luego al dividir cada valor anual de 
medición entre el valor obtenido de la 
curva. Esto creó una serie de índices 
normalizados (cronologías) con media de 
1,0 y varianza homogénea (Cook, 1987). 

Para conocer el periodo de lluvia 
anual que influye de manera significativa 
en el crecimiento del arbolado, se utili- 
zaron datos de la estación climática 
Guadalupe y Calvo, Chihuahua, (Instituto 
Mexicano de Tecnología del Agua, 2000) 
que es la estación más cercana al área de 
estudio. La respuesta climática o Análisis 
de Función de Respuesta se investigó al 
relacionar los índices de anillo con los da- 
tos observados de precipitación emplean- 
do el programa precon (Fritts, 1999). 
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Definido el periodo de respuesta 
climática, con los datos climáticos disponi- 
bles y los índices de anillo de madera 
temprana para el mismo periodo, se 
generó un modelo de predicción. El 
modelo se aplicó para generar la precipita- 
ción de este mismo periodo. Para lo ante- 
rior en la mitad de los datos se realizó una 
calibración entre la precipitación obser- 
vada y la generada, y en la otra mitad se 
realizó una verificación (Fritts, 1991). Final- 
mente, se aplicó la ecuación de transfe- 
rencia del periodo total de datos de preci- 
pitación disponible para desarrollar la 
reconstrucción de precipitación en la 
longitud total de la cronología. A la serie de 
alta frecuencia (resolución anual) se le 
ajustó una curva decenal flexible (baja 
frecuencia) para resaltar eventos secos o 
húmedos (Cook y Peters, 1981). 

La reconstrucción de precipitación 
generada se comparó con otras recons- 
trucciones realizadas para el norte de 
México. Los periodos de sequía 
presentes en la reconstrucción se compa- 
raron con periodos similares en otras 
reconstrucciones de precipitación reali- 
zadas para sitios ubicados en las sierras 
Madre Oriental y Occidental. En la 
medida de lo posible, se analizaron 
archivos históricos que nos permitieron 
verificar la calidad de la reconstrucción 
realizada; asimismo, se analizó la 
influencia de fenómenos atmosféricos 
circulatorios como El Niño, Oscilación del 
Sur (enso, por sus siglas en inglés) y el 
Monzón de Norte América (mna) en la 
variabilidad climática de la región. 

RESULTADOS 

Reconstrucción de la variabilidad climá- 
tica para el suroeste de Chihuahua 

Con el análisis de 53 núcleos de 
crecimiento se desarrollaron cronologías 
de madera temprana, tardía y anillo total 


con una extensión de 349 años, que 
corresponden al periodo de 1657 a 2005. 
Según el programa cofecha la crono- 
logía de madera temprana tuvo una 
intercorrelación entre series alta 
(r = 0,61, p<0.05), lo cual es indicativo 
de una respuesta climática común 
(Holmes, 1983). 

Función de respuesta 

Para definir la respuesta climática 
de la cronología se utilizó el programa 
precon (Fritts, 1999). El programa 
usado produjo un análisis de función 
de respuesta con base en los índices 
de crecimiento y los registros de preci- 
pitación disponibles. En este caso, 
para el periodo 1966-1985, se deter- 
minó a qué periodo de precipitación 
está respondiendo el crecimiento de 
los árboles del sitio. La función de 
respuesta incluye 14 meses de precipi- 
tación de julio del año previo a agosto 
del año en curso, que comprende la 
estación de crecimiento. Tres elemen- 
tos adicionales están asociados con el 
índice de anillo, el crecimiento de 1, 2 
y 3 años previos. 

Esto permitirá determinar si el creci- 
miento actual está influido por el creci- 
miento de años previos y no sólo por la 
variable climática. Un coeficiente signifi- 
cativo indica una relación directa de la 
variable clima o crecimiento previo con el 
ancho del anillo y un coeficiente negativo 
indica una relación inversa. Las líneas 
verticales indican un nivel de significancia 
del 95% (p<0,05) (Figura 4). 

Los meses individuales de precipita- 
ción acumulada que se correlacionaron 
significativamente (p<0,05) con la crono- 
logía de madera temprana son: diciembre 
del año anterior, marzo, abril, mayo, julio 
y agosto del año actual de crecimiento 
(Figura 4). 
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Figura 4. Función de respuesta entre el índices de madera temprana de Pseudotsuga 
menziesii y los registros de precipitación para el período 1966-1985 de la estación 

climática Guadalupe y Calvo, Chihuahua. 


No obstante que la correlación entre 
el índice de madera temprana y precipita- 
ción para los meses de enero, febrero y 
junio presentan una relación positiva pero 
no significativa (p>0,05), se encontró una 
respuesta significativa (p<0,05) al periodo 
estacional enero-junio. Para los tres años 
previos asociados al análisis de 
respuesta se encontró que las condi- 
ciones climáticas previas no influyen de 
manera significativa (p>0,05) en el creci- 
miento anual. 

Reconstrucción de la precipitación 

Al relacionar la cronología de 
madera temprana con datos instrumen- 
tales del periodo de precipitación esta- 
cional enero-junio, periodo 1966-1985, se 
encontró una respuesta significativa entre 
ambas variables con una correlación de 
0.74 (p<0.001), lo que sugiere un poten- 
cial adecuado para la reconstrucción de 
la precipitación (figuras 5 y 6). 


Calibración y verificación 

Para que un modelo de reconstruc- 
ción sea estadísticamente válido, es 
preciso desarrollar dos procesos; uno de 
calibración, realizado en la mitad de los 
datos climáticos disponibles y uno de 
verificación, en la mitad restante. La cali- 
bración es la asociación entre dos varia- 
bles. En este caso el índice de la crono- 
logía de ancho de anillo es una variable y 
los datos promedio mensuales de precipi- 
tación representan los valores de la 
variable climática con los cuales el índice 
se calibrará. El éxito de la calibración se 
mide por el porcentaje de la variación que 
la variable independiente (índice de 
madera temprana) explica de la variable 
dependiente reconstruida (precipitación). 
Por otra parte, la verificación tiene como 
objetivo validar la reconstrucción, mien- 
tras que la calibración únicamente 
examina la relación entre ambas varia- 
bles. 
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Indice de Madera Temprana Precipitación 



Figura 5. Asociación entre el índice de madera temprana y precipitación estacional 
observada para el período enero-junio, que comprende el intervalo de 1966-1985. 
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Figura 6. Gráfico de dispersión del modelo de regresión lineal que muestra la 
variabilidad de la cronología de madera temprana con la precipitación total de enero- 

junio, período 1966-1985. 
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El modelo utilizado para la recons- 
trucción incluyó el periodo total de datos 
disponibles (1966-1985) y los índices de 
madera temprana para el mismo periodo. 
Una vez aplicado el modelo, éste se 
calibró con datos del subperiodo 1974- 
1985 y se verificó en el subperiodo 1966- 
1973. 

Los resultados del modelo indican una 
correlación de r = 0,72 (r 2 = 0,52, p<0,001) 
entre la precipitación observada y la precipi- 
tación reconstruida, y explica 52% de la 
variabilidad en precipitación para el periodo 
de calibración. La verificación presentó una 
correlación de r = 0,88 y explica 77% (r 2 = 
0,77, p<0,001) de la variabilidad. El periodo 
total de precipitación observada y recons- 
truida muestra una r = 0,76 (R 2 = 0.,58, 
p<0,001) (Figura 7). 

Debido a que los subperiodos utili- 
zados tanto para realizar la calibración 
como la verificación del modelo entre la 
precipitación observada y reconstruida 
estuvieron correlacionados significativa- 
mente, se utilizó todo el periodo de datos 
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climáticos disponibles (1966-1985) para 
generar la ecuación de reconstrucción 
que mostró una r = 0,74 (r 2 = 0,55, 
p<0,001) (tablas 1 y 2). 

La subrutina verifys de la Base de 
Programas Dendrocronológicos de la 
Universidad deArizona (dpl, por sus siglas 
en inglés), indicó que tanto la calibración 
como la verificación pasaron la prueba de 
significancia (p<0,05) para la correlación, 
reducción de error, valor de “t” y primera 
diferencia significativa. El modelo de regre- 
sión obtenido para el periodo 1966-1985 
se consideró estadísticamente válido para 
reconstruir el periodo total de la serie 
dendrocronológica. 

El modelo lineal bivariado utilizado 
fue el siguiente: 

Yt = -374,647 + 670,357 * Xt 

Donde: Yt = Valor de precipitación 
reconstruido enero-junio para un año 
específico (mm). Xt = índice de madera 
temprana. 


Reconstruida Obserada IVtedia 



Año 


Figura 7. Precipitación observada (línea sólida) y reconstruida (línea punteada) para el 
período estacional enero-junio, mostrando una r = 0.76 (R 2 = 0. 58, p<0.001) para el 

período total (1966-1985). 
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Tabla 1. Modelo lineal generado para la reconstrucción estacional de precipitación 

invierno-primavera. 


Periodo 

R2 

Coeficientes 

Error 

Prueba de 

Probabilidad 




estándar 

ti ^ tt 


1966 - 1985 

0.55 

a -374.647 

137.74 

-2.71 

0.014 



b 670.357 

134.44 

4.98 

0.000 


Tabla 2. Análisis de varianza del modelo de reconstrucción 


Periodo 

Suma de 
cuadrados 

df 

Cuadrados 

medios 

f 

P 

Regresión 

247357.4 

1 

18 

247357.4 

24.86059 

.000096 

Residual 

Total 

179096.1 

426453.5 

9949.8 




Análisis de la variabilidad climática 

La reconstrucción de precipitación 
invierno-primavera para la reserva Cerro 
El Mohinora, indica gran variabilidad inter- 
anual (alta frecuencia), decenal y multide- 
cenal (baja frecuencia) de los patrones de 
precipitación para la región (Figura 8). 

Sequías severas se reconstruyeron 
para los periodos 1695-1715, 1753-1760, 
1785-1792, 1798-1806, 1819-1830, 1841- 
1870, 1890-1897, 1906-1912, 1924-1941, 
1971-1977 y 1994- 2005. Las sequías 
más prolongadas ocurrieron en los 
periodos 1695-1715, 1841-1870, 1924- 
1941. Las últimas tres décadas del siglo 
XX y los años actuales de la primera 
década del siglo XXI (1971-2005) indican 
un periodo de intensa sequía para la 
región, que continúa hasta el presente, 
afectando el suroeste de Chihuahua. De 
igual manera, se reconstruyeron eventos 
con precipitaciones por arriba de la media 
regional, sobresaliendo los periodos de 


1673-1695, 1735-1753, 1759-1785, 1809- 
1818, 1829-1839, 1898-1906, 1911-1922, 
1941-1969 y finalmente de la década de 
los setentas, a la fecha sobresalen dos 
periodos con precipitaciones por encima 
de la media 1976-1981 y 1990-1993 
(Figura 8). 

DISCUSIÓN 
Respuesta climática 

La reconstrucción de precipitación 
se fundamentó en la madera temprana 
como variable independiente, cuya forma- 
ción está ligada significativamente 
(p<0.05) con la precipitación invernal y 
ésta a su vez por la fase cálida de enso 
(Ropeiewski y Harper, 1989; Stahle et al., 
1998 y Magaña et al., 1999). En México 
esta relación se registra de manera clara 
en la porción de madera temprana de 
Pseudotsuga menziesii en rodales de la 
Sierra Madre Occidental (Stahle et al., 
1998 y Cleaveland et al., 2003). 
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Figura 8. Precipitación estacional reconstruida para los últimos 350 años en la reserva 
ecológica Cerro El Mohinora, Chihuahua. La línea de color gris en el fondo indica la 
variabilidad anual de la precipitación, la curva flexible de color negro es una línea 
suavizada de 10 años para observar de mejor manera los eventos secos (precipitación 
por debajo de la media) y húmedos (precipitación por arriba de la media); la línea recta 
muestra la precipitación promedio (294 mm). Las sequías reconstruidas más intensas 
que fueron verificadas con archivos históricos y otros trabajos de reconstrucción de 
precipitación para el norte y centro de México, se muestran resaltadas con relleno de 

color negro. 


Una explicación a dicha respuesta 
es que la lluvia invernal generalmente es 
de baja intensidad y al no sobrepasar la 
capacidad de infiltración queda almace- 
nada en el perfil del suelo para ser utili- 
zada al inicio de la estación de creci- 
miento (Fritts, 2001). La falta de una 
asociación significativa entre el creci- 
miento y la precipitación de verano, se 
puede atribuir al hecho de que gran parte 
de la precipitación de verano ocurre 
cuando el árbol prácticamente ha cesado 
su división celular. Además, muchas de 
las lluvias superan la capacidad de infil- 
tración del suelo y escurren como flujo 
superficial (García, 1978). 

Resultados similares se han encon- 
trado para reconstrucciones de precipita- 


ción en el norte de México y suroeste de 
EUA (Stahle y Cleaveland, 1993; Stahle 
et al., 1999; Díaz et al., 2002; Cleaveland 
et al., 2003; Cerano, 2004 y Villanueva et 
al., 2007). 


Verificación de la reconstrucción 

En los 350 años de precipitación 
reconstruida, uno de los periodos más 
secos de la reconstrucción se presentó 
de 1695 a 1715, lapso en el que ningún 
año supero la precipitación promedio. En 
esta sequía se presentaron epidemias 
intensas de sarampión y viruela que 
provocaron gran mortalidad entre niños y 
mujeres de Chihuahua (Dunne, 1958). 
Información histórica indica que la región 
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Año 


Figura 9. Influencia del ENSO y el MNA en la variabilidad climática del suroeste de 
Chihuahua, para los siglos XX y XXI. Sequías intensas a consecuencia del efecto de 
ENSO en su fase calida y períodos con precipitaciones por arriba de la media 
originadas por efecto del MNA en su fase positiva y ENSO en su fase fría. 


oeste del estado sufrió la muerte de 1817 
adultos y más de 3000 niños en diversas 
epidemias acontecidas entre 1690 y 1720 
(Gerhard, 1996). Una sequía similar es 
reportada por Villanueva et al., (2006) 
para la región del noroeste de Chihuahua 
y este de Sonora. Esta sequía también 
estuvo presente en los estados de 
Durango y Chihuahua (Stahle et al., 1 999; 
Díaz et al., 2002 y Cleaveland et al., 2003). 

Los años y periodos secos de 1752, 
1758-1759 y 1785-1786, provocaron falta 
de grano, enfermedades y pestes en el 
estado de Chihuahua (Enfield y Tejedo, 
2006). Los años con sequía de 1785 y 


1786 conocidos como “El año del 
hambre”, provocaron una hambruna 
generalizada en el Valle de México, 
debido a la carestía de alimento (Therrell 
et al., 2006), afectando estas sequías 
hasta el norte de México (Florescano y 
Swam, 1995). 

Durante la sequía de 1798-1810 dn 
el norte y centro de México el descon- 
tento social aumentaba debido a los altos 
precios del maíz; este periodo crítico de 
sequía es considerado como un impor- 
tante factor más que contribuyó a la deto- 
nación de la Guerra de Independencia 
(Florescano y Swam, 1995). Este periodo 
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también se ha detectado en diversas 
reconstrucciones paleoclimáticas para el 
norte, noreste, centro de México y Texas 
(Stahle y Cleaveland, 1988; Stahle et al., 
1999; Therrel et al., 2002; Díaz et al., 
2002; Cleaveland et al., 2003; Pohl et al., 
2003 y Villanueva et al., 2006). 

Otro de los periodos secos en el 
siglo XIX es el que ocurrió de 1841 a 
1870, y que causó severos estragos a la 
región norte de México. En este periodo 
ningún año superó la precipitación 
promedio. Esta sequía se extendió para 
otras regiones del país como lo constatan 
reconstrucciones de precipitación esta- 
cional de invierno-primavera para 
Chihuahua y Durango (Stahle etal., 1998; 
Díaz etal., 2002; Cleaveland etal., 2003), 
Coahuila (Cerano, 2004) y Nuevo León 
(Villanueva et al., 2007). Archivos histó- 
ricos señalan que en los años de la 
década de 1850, los estados que se 
vieron más perjudicados debido a las 
intensas sequías fueron Nuevo León, 
Coahuila y San Luis Potosí. La situación 
se complicaba según avanzaba el siglo, 
para los estados de Coahuila, Chihuahua, 
Durango, Nuevo León, Sonora y Tamau- 
lipas ante la carestía de alimentos 
(García, 1997). 

La reconstrucción también registra el 
periodo de 1906 a 1912 con una disminu- 
ción en precipitación. Los años de 1908 y 
1910 se consideran años con sequía 
severa, en especial 1910 que afectó el 
norte de México y originó descontento en 
la población por falta de alimento, lo cual 
es considerado como un factor más que 
propició la detonación de la Revolución 
Mexicana (Florescano, 1980). 

A mediados del siglo XX y específi- 
camente en la década de 1950 se 
presentó un periodo seco que afectó el 
norte y centro de México. Este periodo 
seco se ha detectado en reconstruc- 
ciones paleoclimáticas para todo el país 
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(Díaz etal., 2002; Cleaveland et al., 2003; 
Cerano, 2004 y Villanueva et al., 2006, 
2007). Para el caso específico de este 
estudio la reconstrucción desarrollada no 
indica la presencia de sequías severas 
para este periodo; no obstante, el periodo 
1924-1941 registró condiciones de preci- 
pitación por debajo de la media, constitu- 
yendo una intensa sequía con mayor 
impacto social y económico en esta 
región que la que se presentó más tarde 
en la década de 1950 y que fue más 
generalizada. 

Finalmente, el periodo 1971-2005 
muestra una sequía prolongada donde 
solo los años de 1978-1982 y 1992-1993 
excedieron la precipitación promedio. La 
década de 1970 se ha caracterizado 
como de extrema sequía para el norte de 
México, en especial los años de 1971, 
1972, 1974, 1975, 1976 y 1977 (Flores- 
cano, 1980). 

Estos eventos secos también se 
encuentran presentes en reconstruc- 
ciones paleoclimáticas para Durango y 
noreste de México (González, 2003; 
Cerano, 2004 y Villanueva et al., 2007). 
La sequía de 1994-2005 representa una 
de las sequías más intensas de los 
últimos 300 años para el suroeste del 
estado de Chihuahua y hasta la fecha 
continúa afectando la región. 

En los tres siglos y medio de precipi- 
tación reconstruida se observan fuertes 
sequías periódicas cada 100 años y 
sequías menos intensas pero más prolon- 
gadas a mitad de siglo. En otras recons- 
trucciones en el norte de México que se 
extienden para un periodo más prolon- 
gado, se observa un comportamiento 
similar con sequías intensas a finales de 
1500 e inicios de 1600 (Villanueva et al., 
2005). La sincronía de las sequías 
reconstruidas para los últimos 600 años 
en el norte de México permite inferir que 
si en el futuro este comportamiento 
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continúa, es probable la ocurrencia de 
intensas sequías en la década del 2050 y 
finales del 2100; la sequía de finales del 
siglo XX aún está presente, aunque para 
algunas regiones del norte de México 
como el suroeste de Chihuahua, se ha 
mostrado con mayor intensidad. 

Influencia de patrones circulatorios El 
Niño, Oscilación del Sur y el Monzón 
de Norteamérica 

Las condiciones climáticas de la 
región de estudio están fuertemente 
moduladas por fenómenos como el enso 
y el mna. El mna produce más del 60% de 
la precipitación anual para el suroeste del 
estado de Chihuahua, en el periodo de 
verano (eric ii, 2000). El enso propicia un 
aumento de las lluvias durante el invierno 
en su fase cálida y una disminución de la 
precipitación durante el verano, en su 
fase fría (Niña) ocasiona lluvias normales 
o por arriba de la media en el verano 
(Magaña et al., 2003). 

Algunos estudios recientes indican 
que el mna en la región de Sonora y 
Sinaloa es débil durante veranos Niño. 
El inicio del mna en el noroeste mexicano 
puede presentar adelantos o retrasos 
sustanciales (30 días), con respecto a la 
fecha de inicio de la temporada de 
lluvias (17 de junio). Un inicio tardío del 
mna generalmente significa una tempo- 
rada de lluvias deficientes, lo cual 
parece ocurrir en años Niño (Higgins et 
al., 1999). 

Englehart y Douglas (2002) 
indican que la región oeste de la Sierra 
Madre Occidental presenta correla- 
ciones bajas con enso, durante la fase 
positiva de la Oscilación Decadal del 
Pacifico (pdo, por sus siglas en inglés), 
fenómeno que es un estimativo de la 
variabilidad del mna. En contraste, se 
presentan correlaciones significativas 


entre el enso y la precipitación durante 
la fase negativa del pdo. 

Eventos intensos de enso provocan 
que la señal del mna sea negativa, origi- 
nando una caída en la precipitación y por 
ende sequías de gran magnitud en el 
periodo primavera-verano. En contraparte 
una señal negativa del enso propicia 
condiciones de lluvia normales en el 
verano, esto atribuido a que no se inhibe 
o se retrasa la señal del mna. 

La relación de ambos fenómenos es 
clara en la variabilidad de la precipitación 
reconstruida para el siglo XX y XXI. La 
sequía reconstruida en 1908 da inicio 
debido a los eventos intensos del Niño de 
1902, 1904 y 1905; posteriormente se 
reconstruye un periodo libre de sequías 
de 1911-1921, atribuible a ningún registro 
de ocurrencia del fenómeno enso, 
presentándose el mna en fase positiva. 
Otra fuerte sequía es la ocurrida en el 
periodo de 1935 y que fue propiciada por 
eventos El Niño de 1929, 1930 y 1932; 
posterior a estos eventos comienza un 
incremento en la precipitación alcan- 
zando sus puntos máximos en los años 
de 1942 y 1944, que precisamente son 
años del Niño en su fase fría, originando 
un incremento en la precipitación de 
verano (Figura 9). 

La sequía reconstruida para la 
década de los setenta da inicio a conse- 
cuencia del evento de El Niño de 1969, 
que ocasionó intensas sequías en el 
verano de 1970. En el invierno de 1970 se 
registra el mismo fenómeno, que aunque 
se presentó en su fase fría, no se regis- 
tran precipitaciones por arriba de la 
media, dando con esto paso al año más 
seco de la década de los setentas, el 
verano de 1971. Finalmente, la constante 
aparición de El Niño en la década de 
1980 y 1990 ha llevado a declarar zonas 
de desastre a gran parte de los estados 
del norte del país por la falta de lluvias. 
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Los eventos intensos El Niño 1982-1983, 
1986-1987, 1994-1995 y 1997-1998 han 
ocasionado las intensas sequías recons- 
truidas para las últimas dos décadas del 
siglo XX y los años actuales de la primera 
década del siglo XXI (Figura 9). 

CONCLUSIONES 

La similitud o sincronización de 
eventos secos reconstruidos para Cerro 
Mohinora y otros sitios del norte y centro de 
México, permiten determinar la cobertura 
de los fenómenos climáticos. Las fuertes 
sequías que afectaron a esta región no son 
eventos locales, sino que afectan amplias 
regiones a consecuencia del impacto de 
patrones de circulación global, lo cual 
determina la variabilidad de la precipitación 
en grandes áreas de México. 

En los 350 años de precipitación 
reconstruida se detectan fuertes sequías, 
tanto al final como al principio de cada 
siglo, es decir, 1698 a 1718, 1798 a 1807, 
1890 a 1908 y 1994 a 2005, así como 
sequías de menor intensidad pero de 
mayor duración cada mediados de siglo, 
como en 1650, 1750, 1850 y la más 
reciente de 1950. Actualmente, para el 
suroeste de Chihuahua se vive una 
sequía que inició en la década de 1990 y 
continúa hasta el presente. Si esta 
tendencia continúa se esperaría una 
sequía en el 2050. No obstante, de 1900 
a la fecha las sequías se han repetido 
cada 20 años; situación que se ha inten- 
sificado después de 1970 con sequías 
cada década, probablemente afectada 
por una mayor variabilidad de enso. 

Es fundamental explicar los eventos 
de baja frecuencia que se han presentado 
de manera periódica cada 50 o 100 años 
en el norte de México, y contar con los 
elementos técnicos para modelar posi- 
bles escenarios futuros del clima, lo que 
es de gran importancia para la planeación 
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del adecuado aprovechamiento de los 
recursos hídricos con que cuenta esta 
región. 

Los fenómenos circulatorios enso y 
mna modulan la variabilidad climática de 
la región del suroeste de Chihuahua. Sin 
embargo, es importante poner especial 
énfasis en eventos enso, los que actual- 
mente tienen cierta predictibilidad. El 
conocimiento con antelación de la inci- 
dencia de este fenómeno permitirá inferir 
la variabilidad hidroclimática para el suro- 
este del estado de Chihuahua, y tomar las 
medidas preventivas para tal efecto, parti- 
cularmente en términos de planeación de 
los recursos hídricos. 
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Influencia de la edad en el crecimiento 
de vinal ( Prosopis ruscifolia Burkart.), en la 
Provincia de Santiago del Estero, Argentina 

Age influence in growth of vinal (Prosopis ruscifolia 
Burkart.), in the Province of Santiago del Estero, Argentina 

Ana María Giménez, N. Ríos, P. Hernández y J. G. Moglia* 

RESUMEN 


El vinal (Prosopis ruscifolia) es un árbol de uso múltiple, constituye bosques halófilos, colonizador 
por excelencia de ambientes inestables sometidos a inundaciones. El objetivo del trabajo fue analizar las 
posibilidades de crecimiento de árboles individuales en función de la edad, estimar el tiempo de tránsito 
y el turno biológico. El estudio fue realizado con 15 individuos de Prosopis ruscifolia (Mimosaceae), de 
Loreto, Departamento San Martin, Santiago del Estero. Argentina. La elección de los árboles se efectuó 
siguiendo los criterios de: árbol dominante, sano, de DAP superior a 15 cm y fuste sin defectos. Para el 
estudio de crecimiento se trabajó con las secciones transversales obtenidas a las alturas de 0,30, 1,30 
m y hasta la primera ramificación. La marcación y medición de anillos se efectuó con el Equipo Compu- 
tarizado ANIOL y el programa CATRAS, en las orientaciones Norte, Sur, Este y Oeste. Prosopis rusci- 
folia es una especie de crecimiento rápido en la Región Chaqueña Seca, con anillos de espesor 
promedio de 5,66 mm (2,20-9,8) con un incremento anual diametral (IA) de 0,93 cm e IMA de 1 ,075 cm. 
El diámetro mínimo de corta es de 25 cm, lo que equivale a un árbol de 23 años. El IMA e IA calculado, 
en volumen de fuste alcanza el turno biológico de corta a la edad de 45 años. Se requiere continuar estu- 
diando el comportamiento de la especie en otras condiciones de crecimiento. 

PALABRAS CLAVE: 

Crecimiento, edad, Prosopis ruscifolia. 


ABSTRACT 

Vinal ( Prosopis ruscifolia) is a multi-purpose tree, halophyte, colonizer par excellence of unstable 
environments proved to fiooding. The research objective of this work are: Analyseng the growth prospects 
of individual trees in relation to age, and estímate biological tum short. Fifteen trees of Prosopis ruscifolia , 
Mimosaceae, were selected from the natural forest of Loreto, San Martin Department, Santiago del Estero 
province, Argentina. Trees selection were based on the following criteria dominant canopy, healthy, with 
good form, more than 15 cm of DBH and free of defects. Transverse sections were obtained at different 
heights from 0,30, 1 ,30 m, and until the first branching. Growths rings were determined in four radii at north, 
south, east and west orientations. Rings width was measured with the Computerized Equipment ANIOL 
and the program CATRAS. Prosopis ruscifolia is a rapid growth tree in the Semi arid Chaco región, with 
ring average thickness of 5,66 mm (2,20-9,8). Mínimum cut diameter is 25 cm, which is equivalent to 23 
years. IMA and Al calculated by volume of stem reache the biological cut shift at the age of 45 years. The 
behaviour of this species in other growing condítions requires further study. 

KEY WORDS: 

Growth, age, Prosopis ruscifolia. 
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INTRODUCCIÓN 

Los bosques de Prosopis de las 
zonas áridas y semiáridas de Argentina 
han sido fuente de subsistencia para 
numerosas comunidades humanas 
durante varios siglos, y hasta el presente 
siguen siendo explotado por los habi- 
tantes de la zona (Villagra et al., 2002). 
Desde el punto de vista económico, los 
rodales de Prosopis fueron sometidos a 
una explotación con concepto minero, sin 
ajustar la velocidad de extracción a la 
velocidad de renovación de los recursos 
extraídos. A esto se le suma la escasa 
capacidad de recuperación natural que 
tienen estos sistemas dadas las caracte- 
rísticas ambientales extremas a las que 
han sido sometidos y a las condiciones 
actuales de alta degradación. El conoci- 
miento de la estructura y dinámica de 
estos bosques permitirá plantear estrate- 
gias de manejo y recuperación de los 
mismos. 

Mediante el Decreto Núm. 85.584 
del 1 de marzo de 1941, se declaró al 
vinal como plaga. Se debió a las caracte- 
rísticas de especie colonizadora e inva- 
sora de zonas deforestadas y usadas en 
actividades agropecuarias. Tal decreto, 
comprendido dentro de la Ley Núm. 4863 
en su Artículo 2 o encomienda al Poder 
Ejecutivo Nacional efectuar la nomencla- 
tura de las plagas incluidas en las dispo- 
siciones de dicha norma, asimismo, a 
modificarla cuando lo estime conveniente. 
Con fecha del 14 de junio de 2007, el 
presidente de la nación Argentina 
mediante el Decreto 746/2007, deroga el 
Núm. 85.584. 

El vinal ( Prosopis ruscifolia) habita 
en dos distritos del Chaco Argentino: el 
oriental (húmedo) y en el occidental 
(semiárido). Es una especie abundante, 
de hasta 12 m de altura, de uso múltiple, 
que constituye bosques halófitos con un 
subbosque de quenopodiaceae (Cabrera, 


1971). En la Provincia de Santiago del 
Estero se distribuye naturalmente en 
áreas relacionadas a los derrames de los 
ríos Dulce y Salado. Se destaca por su 
resistencia máxima a la salinidad junto a 
P vinalillo, P. pugionata, P strombulífera, 
P reptans y P sericantha. 

El vinal se instala sólo en bordes de 
ambientes pantanosos de agua dulce y 
borde de salitrales, su óptimo ecológico 
es el deslinde entre pantano y tierra firme. 
En ocasiones forma densos matorrales 
cuyo factor determinante son los suelos 
formados por capas aportadas por 
derrames fluviales, salinos, de bajo conte- 
nido en materia orgánica y textura limo- 
arcillosa. 

Morello et al. (1971) describen la 
presencia de vinal asociada directamente 
a la precipitación anual. Está ausente a 
valores mayores a 1000 mm; entre los 
rangos de 900-1100 mm, es fuerte su 
presencia; entre 750-900 mm, se 
encuentra en zonas de desagües y 
suelos halohidromórficos; entre 500-750 
mm su presencia se manifiesta sólo en 
las márgenes de ríos permanentes y con 
menos de 500 mm es manifiesta la 
ausencia de vinal. Actualmente cubre 
8,106 ha, habiendo invadido intensas 
áreas por efecto de derrames de los ríos 
Teuco, Bermejo, Salado y Dulce. 

La presencia del vinal se asocia a 
una acumulación hídrica en el suelo cuyo 
perfil presenta síntomas de anaerobiosis 
temporaria (Adámoli et al., 2001). Es 
evidente que las especies que soportan 
las inundaciones han evolucionando 
respondiendo al ritmo de procesos 
geomorfológicos. 

Su principal fuente de expansión es 
el ganado, el desmonte y la presencia de 
áreas agrícolas abandonadas que facili- 
taron su establecimiento. El agua juega 
un papel importante, ya que arrastra y 
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deposita las semillas. La erosión hídrica 
produce colmatación de sitios inundables, 
permitiendo la instalación del vinal (Karlin 
y Palacios, 1988). Es la leñosa invasora 
más importante del país, continuando su 
proceso de expansión. 

Crea condiciones microambientales 
para la instalación de especies forestales 
de mayor valor como Aspidosperma 
quebracho-blanco y Schinopsis lorentzii, 
generalmente crece acompañado por 
Prosopis affinis, P. vinalillo, P kuntzei y 
P. hassleri. 

Extensas zonas del Chaco argentino 
y boliviano están sometidas a desbordes, 
inundaciones, originadas por crecidas y 
desviación de los ríos que lo cruzan, 
abriéndose nuevos cauces o rellenán- 
dose otros. Se adaptó el vinal a ese dina- 
mismo, siendo el árbol colonizador por 
excelencia de los ambientes inestables 
originados por acción pluvial, lo que exige 
investigación sobre su biología. 

El Grupo de Estudios sobre Ecología 
Regional de la Facultad de Ciencias 
Exactas y Naturales de la uba, geser, 
implemento una nueva visión sobre los 
bosques de vinal. Adamoli et al. (2001) 
diseñaron en Formosa un modelo produc- 
tivo para el manejo de vinal en áreas 
improductivas correspondientes a zonas 
agrícolas-ganaderas abandonadas. 

En los últimos años ha surgido un 
gran interés por esta especie que tiene 
potencialidades de uso múltiples. La prin- 
cipal fuente de valor actual es la madera, 
que se caracteriza por su color castaño 
amarillento, dura y pesada; con peso 
específico de 0.83 kg/dm3 y excelente 
estabilidad dimensional. Es usada local- 
mente como leña y carbón, por su 
elevado poder calorífico de 4400 
kgcal/kg. De características tecnológicas 
similares a los otros Prosopis, es apta 
para carpintería de obra, tejuelas, pisos, 
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encofrados y carpintería fina, mesas, 
sillas, aberturas y ebanistería (Piarfon, 
2006) 

Constituye una alternativa de uso 
como combustible con rendimientos de 
100-150 ton/ha (Gorleri, 1987). Según la 
Base de Datos Productos Forestales No 
Maderables, son siete los usos princi- 
pales: alimenticio, tintóreo, forrajero, melí- 
fero, productos bioquímicos, medicinal y 
artesanal. Los frutos y hojas son exce- 
lente forraje. La vaina del vinal contiene 
30% de glucosa, 5,7% de proteína, 18% 
de almidón y otros hidratos de carbono 
(Leonardis, 1949), elaborándose 
alimentos tradicionales como patay, aloja 
y añapa. Como uso medicinal destaca las 
hojas que se emplean en el control de 
glucemia, diabetes, desinfectante 
externo, antiséptico, desinflamante y 
conjuntivitis. 

Esta especie, durante años ha sido 
ignorada y hoy se miran con atención las 
áreas marginales donde habita, empren- 
diéndose planes de manejo y uso múltiple 
como los iniciados en Santos Lugares por 
el proyecto Piarfon. 

El manejo sustentable en un ecosis- 
tema forestal exige el conocimiento de la 
renta que produce el capital monte, repre- 
sentado por el crecimiento del bosque. La 
información referida al crecimiento es un 
subsidio básico para la elaboración y apli- 
cación de un plan de manejo, necesario 
además para determinar ciclos de corta y 
para regular la producción del bosque. 
Particularmente el estudio del crecimiento 
diamétrico permite analizar cómo ocurre 
el movimiento de los árboles a través de 
las sucesivas clases de diámetro y 
estimar el número de años requeridos 
para que los individuos que se encuen- 
tran en una clase de tamaño pasen a la 
siguiente. Tales estimaciones sirven para 
calcular el tiempo de tránsito y el módulo 
de rotación o ciclo de corta, determina- 
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ción fundamental para la obtención de un 
rendimiento sostenido y para la organiza- 
ción de la masa en el tiempo y en el 
espacio (Araujo, 1993). 

El crecimiento se puede estimar a 
partir de datos de inventario o por análisis 
epidométrico. La carencia de parcelas 
permanentes para el estudio de creci- 
miento, determinó el uso de la metodo- 
logía del análisis epidométrico a fin de 
sentar bases para el conocimiento de las 
potencialidades económicas y de manejo 
de la especie. 

En el marco del proyecto financiado 
por el cicyt-unse: “Crecimiento y calidad 
de madera de nativas del Chaco Semiá- 
rido”, surge el interés por estudiar pará- 
metros de crecimiento de vinal en 
Santiago del Estero, Argentina. 

Los objetivos del presente trabajo 
son analizar las posibilidades de creci- 
miento de árboles individuales de 
Prosopis ruscifolia, en función de la edad, 
estimar el tiempo de tránsito y el turno 
biológico. 

MATERIAL Y MÉTODO 

El estudio fue realizado con 15 indi- 
viduos adultos de la especie Prosopis 
ruscifolia, perteneciente a la familia Mimo- 
saceae, que provienen del bosque nativo 
de Loreto, Departamento San Martín 
Provincia de Santiago del Estero, Argen- 
tina. 

El rodal objeto de estudio es un 
vinalar invasor que se acerca al área de 
inundación del río Saladillo, ocupando 
zonas bajas de relleno y creando ecosis- 
temas marginales. La vegetación es un 
bosque bajo-arbustal, donde el vinal se 
expresa básicamente en su forma arbus- 
tiva, siendo escasos los individuos arbó- 
reos. La asociación de leñosas arbóreas 


está conformada por las siguientes espe- 
cies: Prosopis ruscifolia, Prosopis vina- 
lillo, Tabebuia nodosa, Celtis tala, Vallesia 
glabra, Prosopis alba y Prosopis nigra. El 
estrato arbóreo tiene una altura promedio 
de 11 metros. 

El estrato arbustivo está formado por 
cactáceas de los géneros: Cereus y 
Opuntia y quenopodiáceas ( Heteros - 
tachys ritteriana, Allenrolfea patagónica, 
A. vaginata, Suaeda divaricata, Atriplex 
argentina y A. lampa). 

Los árboles fueron escogidos según 
criterio de: individuos con dap superior a 
15 cm; posición social dominante, fuste 
sin defectos aparentes; de copa simétrica 
y sin síntomas externos de enfermedad. 
Cada árbol constituyó el centro de una 
parcela circular con un radio de 17,84 
metros (1 000 m 2 de superficie) midién- 
dose el dap, la altura de fuste (definida 
por la primera ramificación principal) y 
altura total. 

El árbol seleccionado fue abatido, se 
marcó la orientación norte en el fuste, se 
extrajeron rodajas de 5 cm de espesor, a 
alturas de 0,30 m, 1,30 m y a partir de 
aquí, a cada metro hasta llegar a la altura 
del fuste. Las muestras fueron prepa- 
radas con cepilladora, lijadora de banda y 
lijadora orbital con juego de lijas de 
granulometría de 600 a 100, hasta conse- 
guir un perfecto pulido. 

La marcación y medición de anillos 
se efectuó con el Equipo Computarizado 
aniol y el programa catras (Aniol, 1991) 
con una precisión de centésima de milí- 
metro, sobre los radios: Norte, Sur, Este y 
Oeste. 

Se realizó el análisis epidométrico 
del fuste de los árboles individuales, que 
es la forma más directa y precisa para 
estudiar la evolución de un árbol, deter- 
minar su crecimiento e interrelaciones. Se 
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calculó dap, sn, altura del fuste y volumen 
de fuste mediante Smallam. 

Para el estudio de las tendencias de 
crecimiento se usó el método basado en la 
edad biológica de los árboles (Giménez, 
1998), trabajándose con individuos de 
edades diferentes. El crecimiento (incre- 
mento gradual de un organismo en un 
periodo de tiempo) depende del parámetro 
y del periodo de tiempo considerado. En el 
presente trabajo se expresa el crecimiento 
en dap, sn y volumen del fuste. Se parte de 
los espesores en relación a la edad deter- 
minada según el marcado y conteo de 
anillos. Se calculó para cada variable el 
incremento medio anual (ima) e incremento 
anual (ia), dado la posibilidad de datos de 
espesores anuales. 

Para el análisis de los datos, el 
espesor medio de los anillos de cada 
árbol se alineó con los de otros individuos 
según la edad biológica y no cronológica, 
fue promediado para producir la curva de 
crecimiento del árbol base. Se usó la 
técnica de suavizado de las series indivi- 
duales de anchos de anillos, a fin de 
eliminar la variación climática de las 
series dendrocronológicas individuales, 
previo a la modelación del crecimiento. 
Se aplicó un filtro de suavizado exponen- 
cial (ssps). 
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La estandarización de la curva de 
crecimiento asume que la forma de la 
estructura a cualquier edad biológica es 
independiente al periodo de tiempo 
durante el cual se produce (Juárez et al., 
2003). Con el fin de visualizar las rela- 
ciones existentes entre las variables de 
referencia y el ajuste a una función mate- 
mática se analizaron los modelos de 
regresión para las variables dendromé- 
tricas en referencia a la edad. Los 
modelos de ecuaciones de regresión utili- 
zados para el ajuste de curvas fueron: 
lineal, exponencial, logarítmico, potencial 
y polinómico. Se seleccionó el modelo 
que mejor R 2 presentó. 

RESULTADOS 

El análisis epidométrico de los 
árboles en estudio permitió elaborar las 
tendencias de crecimiento basadas en la 
edad biológica expresadas en la tabla 1 . 

El patrón de marcación de anillos de 
crecimiento fue descrito por Giménez et 
al. (2005 a y b). En la figura 1 se expresa 
la distribución de espesores de anillo y la 
edad. El espesor promedio es 5.32 mm 
(2,20-9,8), S: 1,69 y CV%: 31,8, corres- 
ponde a anillos de categoría 3 (medianos) 
(Giménez et al., 2007). 


Tabla 1 . Relaciones de crecimiento y sus respectivos modelos de regresión. 


RELACIÓN DE VARIABLES 

ECUACIÓN 

R2 

Espesor de anillos (y) 

y = -1 E-05x4 + 0,0012x3 - 0,0379x2 + 0,3579x + 5,0104 

R2 = 0,4711 

y edad (x) 



DAP y edad (x) 

y = -0,0022x2 + 1 ,11 93x + 0,4318 

R2 = 0.991 

Sección N y edad (x) 

y = 0,7507x2 + 4,21 07x - 7,0216 

R2 = 0,9998 

Altura fuste/ edad (x) 

y = 0, 9491 Ln(x) -0,7593 

R2 = 0,9234 

Volumen /edad (x) 

y = 1,8015x2 + 22, 861x- 212,42 

R2 = ,9917 

IMA en vol / edad (x) 

y = -0,9424x2 + 255, 27x + 213,9 

R2 = 0,9954 

ICA en vol /edad (x) 

y = -0,9849x2 + 434, 13x + 961,2 

R2 = 0,9438 
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El anova de la relación referida al 
modelo es altamente significativo entre 
individuos (6,4315 E-10) para un nivel 
de probabilidad del 0,95. El espesor de 
anillos tiene un máximo entre 5-6 años, 
para luego decrecer. Esta variable es 
determinante para el cálculo del creci- 
miento y resume la capacidad biológica 
de la especie influida por sus caracte- 
rísticas genéticas y su ¡nterrelación con 
el ambiente, lo que de alguna manera 


explica la gran variabilidad individual 
de los espesores de anillos. Razón por 
la cual se aplica el suavizado exponen- 
cial. 

En función del espesor de anillos se 
calcula el dap promedio, que ajusta a una 
ecuación polinómica de segundo grado 
con R 2 : 0.99 (Figura 2). El diámetro 
mínimo de corta establecido por ley, de 
25 cm, se alcanza a los 23 años. 



Figura 2. DAP en función de la edad. 
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Se calculó el incremento anual (ia) y el 
incremento medio anual (ima) en referencia 
a la variable dap. Ambas curvas se inter- 
ceptan a los 9 años, lo que indica un máximo 
de incremento medio a esa temprana edad. 
A los 25 cm al dap le corresponde un ia de 
0,93 cm e ima de 1,075 cm (Figura 3). 

El ia, al principio, aumenta en forma 
rápida para decrecer luego del punto de 
culminación lentamente. Esta tendencia 


básica puede verse enmascarada por 
influencias distorsionantes, no obstante, 
para todo el lapso de vida la curva es 
inconfundible. 

Para calcular el tiempo de tránsito 
de una clase diamétrica (5 cm) a otra se 
partió de la curva de espesor de anillos 
promedio. Se estableció que se tarda 5 
años para pasar de una clase a la 
siguiente. El incremento periódico por 


Incremento Periódico Diamétrico 



o 0,95 


— IPAD 


i 1 1 1 r 

11 16 21 26 31 


Clase diamétrica 


Figura 4. Incremento periódico diamétrico. 
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Sección normal /Edad 



Figura 5. Sección normal y edad. 


clase diamétrica disminuye con la edad 
con valores entre 1,15 a 0,98 cm, en inter- 
valo de 5 a 46 cm al dap. 

La relación entre la sección normal y 
la edad se ajusta a una ecuación polinó- 
mica de segundo grado, obteniéndose 
para un fuste de 48 años una sección de 
0.19 m 2 (Figura 5). El incremento anual (ia) 
y el incremento medio anual (ima) no se 
interceptan a la edad máxima estudiada. 


La altura de fuste se incrementa 
hasta los 20 años, para después mante- 
nerse constante (Figura 6). El ia e ima se 
interceptan a la edad de 13 años con un 
valor de incremento de 0,12 m. El fuste en 
vinal es corto con valores promedio de 
2,50 m, lo que limita el aprovechamiento 
de la madera en piezas largas. 

La mejor forma de expresar el creci- 
miento de árboles maduros es en términos 



Figura 6. Evolución de la altura del fuste. 
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Volumen de fuste 



Vol fuste 


Figura 7. Relación entre el volumen de fuste y la edad. 


de volumen. El volumen de fuste es un 
parámetro válido para expresar integral- 
mente el crecimiento del árbol. La rela- 
ción volumen de fuste y la edad, ajusta a 
una ecuación polinómica, con volumen de 
0,495 m 2 a los 48 años (Figura 7). 

El crecimiento en volumen de fuste 
se expresa en la figura 8 donde se mani- 
fiesta la evolución de los incrementos ia e 


ima. La culminación del crecimiento medio 
se logra más tarde que el crecimiento 
corriente. El máximo de ia jamás es alcan- 
zado por el ima. Ambas curvas se inter- 
ceptan a la edad de 45 años, por lo cual se 
estima que el vinal encuentra su edad de 
culminación biológica. Esta tendencia es 
muy promisoria ya que en la región la edad 
de culminación de otras especies se 
produce a edades superiores a 100 años. 


IMA e IA en Volumen de fuste 



IMA 

A 


Figura 8. Evolución del lAe IMA en volumen. 
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DISCUSIÓN 

Adámoli et al. (2001) indican que el 
vinal se distingue entre los Prosopis por 
su alta velocidad de instalación en 
nuevas tierras y por tener un crecimiento 
importante (considerando el tipo de 
madera que produce), semejante al del 
algarrobo blanco ( Prosopis alba). Las 
comunidades vegetales que forma son 
casi puras, de muy alta densidad, con 
escasa diversidad y casi nulo desarrollo 
del estrato herbáceo. Esto evidencia la 
gran posibilidad de la especie en zonas 
marginales donde crece. 

El crecimiento de vinal es superior al 
de las otras especies de la zona. El 
espesor medio de los anillos (5,32 mm) 
resulta ser el más promisorio de las 
leñosas arbóreas de la región chaqueña 
semiáñda de Argentina. 

En el género Prosopis se destacan: 
Prosopis alba con un espesor de 4,05 mm 
(2,11-5,18) (Giménez et al., 1998); 
Prosopis kuntzei de 3,26 mm (1,9-4, 5) 
(Giménez et al., 1997); Prosopis nigra de 
3,289 mm (1,14-8,59) (Giménez et al. 
2000 ). 

Si se compara el espesor de anillos 
con las dos especies principales de la 
región: Schinopsis quebracho-colorado 
2,19 mm (Giménez Ríos, 1999) y Aspi- 
dosperma quebracho-blanco con 2,2 mm 
(Moglia, 2000), el vinal tiene amplia 
ventaja. 

Juárez et al. (2007) estudiaron el 
crecimiento en diámetro del vinal en la 
localidad Santos Lugares, Santiago del 
Estero, bajo el enfoque de modelos linea- 
les mixtos, empleando para el crecimiento 
promedio poblacional un modelo polino- 
mial de segundo orden e incluyendo 
componentes aleatorios en cada uno de 
sus parámetros para contemplar la varia- 
bilidad entre árboles. Se encontró que el 


vinal crece 0,49 cm por año y que la edad 
de culminación biológica para el diámetro 
es alrededor de los 11 años. Este valor es 
ligeramente inferior al obtenido en el 
presente trabajo, posiblemente por 
influencia del sitio. 

Un aspecto interesante es el tiempo 
de tránsito por clase diamétrica que es de 
5 cm en 5 años. Este valor contrasta con 
los calculados para las otras especies 
(6-12 años). 

La vida óptima de una masa forestal, 
o turno óptimo, puede definirse de 
muchas maneras. El turno biológico en 
función del volumen de fuste es de 45 
años en vinal, valor bajo que despierta 
gran expectativa para la región, si se le 
compara con otras especies del género. 
Prosopis kuntzei (itín) generó gran interés 
como alternativa económica por las 
características tecnológicas de su 
madera, con un turno biológico calculado 
en 59 años para la variable volumen de 
fuste (Ríos et al., 2001). Sobre la capa- 
cidad de producción del chaco semiárido, 
algunos autores determinaron turnos de 
rotación muy largos, en torno a los 80 
años; no obstante, Brassiolo (2004) 
propone turnos de 50 años según el nivel 
de degradación de la masa. 

En otras regiones boscosas argen- 
tinas (provincias fitogeográficas del 
monte, espinal y prepuna) hay referencia 
dendrométricas de otras especies de 
Prosopis que presentan bajas tasas de 
crecimiento, especialmente durante los 
primeros años de vida. 

En Prosopis flexuosa (Villagra et al., 
2002) el mayor aporte en madera se da 
entre los 40 y 45 años y el turno biológico 
de corta sería de aproximadamente 130 
años. Este resultado implica que es nece- 
sario replantear la mayoría de los 
proyectos de producción de madera 
rolliza. La demarcación clara de anillos 
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anuales en P. flexuosa y P. ferox brinda la 
oportunidad de estudiar la estructura y 
dinámica de bosques de zonas áridas y 
semiáridas a través de métodos dendro- 
cronológicos (Villagra et al., 2002). 

Para Prosopis caldenia el turno 
biológico de corta, en función del área 
basal, osciló entre los 66 y 100 años de 
edad. Estos valores muestran tasas de 
crecimiento elevadas, teniendo en cuenta 
que el caldén es una especie de madera 
dura que crece en ambientes semiáridos 
con escasa precipitación anual (Boggino 
y Villalba, 2004). 

Villalba et al. (2006) analizan el turno 
biológico de corta de una especie 
emblema de la región fitogeográfica de 
las Yungas, el Cedrela lilloi, a partir del 
crecimiento promedio corriente (ica) y 
medio anual (ima), calculados en función 
de los incrementos anuales en área basal 
a partir de las 10 series de anchos de 
anillos. La edad del turno biológico de 
corta se alcanza aproximadamente a los 
150 años de edad, donde la curva de ima 
corta la curva de ica. 

Los estudios mencionados son los 
primeros realizados en el país sobre el 
tema y representan un invalorable aporte 
al conocimiento. Posicionan al vinal en 
una situación de privilegio en la produc- 
ción de madera, como la especie arbórea 
del chaco semiárido que presenta un 
turno biológico de corta a edad temprana. 
Esto implica la posibilidad de ¡mplementar 
planes de manejo forestal con rotaciones 
más cortas que lo estipulado en la región, 
además de contribuir científicamente en 
el conocimiento biológico de la especie 
tan frecuente en la zona, pero de la cual 
se carece de información básica para su 
mejor aprovechamiento. Esta caracterís- 
tica, aunada con las bondades de la 
madera, sus usos múltiples y su 
presencia en áreas salinas, inundables, 
marginales, hacen de esta especie un 
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valuarte importante para el manejo 
sustentable de la región chaqueña argen- 
tina. 

CONCLUSIONES 

Prosopis ruscifolia es una especie 
de crecimiento rápido en la región 
chaqueña seca, con anillos de espesor 
promedio de 5,66 mm (2,20-9.8) El 
diámetro mínimo de corta es de 25 cm, lo 
que equivale a un árbol de 23 años. El ima 
e ia calculado en volumen de fuste 
alcanza el turno biológico de corta a la 
edad de 45 años, para un tiempo de trán- 
sito entre clases diamétricas en 5 años. 
Se requiere continuar estudiando el 
comportamiento de la especie en otras 
condiciones de crecimiento. 
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RESUMO 


No municipio de Iranduba, Estado do Amazonas - Brasil há grande consumo de madeira pelas 
comunidades locáis para fabricagáo de carváo, e pelas olarias na queima para produgáo de tijolos. 
Diante deste panorama, a EMBRAPA vem desenvolvendo um projeto visando implementar plantíos 
ñas comunidades deste municipio Neste contexto, o presente trabalho tem como objetivo avaliar as 
espécies Acacia aurículiformis (exótica) e Ormosia paraensis (nativa) a partir das variáveis dendromé- 
tricas, determinagáo da densidade básica, poder calorífico e análise ¡mediata. O experimento foi 
implantado na Estagáo Experimental da EMBRAPA no municipio de Iranduba. O delineamento foi em 
blocos ao acaso com trés repeticóes perfazendo um total de 6 parcelas compostas de 25 plantas 
cada. Foram retiradas cinco árvores dentre as nove centráis da área útil. Destas foram retirados os 
discos com 5cm espessura e desdobrados em corpos de prova para determinar a densidade básica 
e demais testes. Foi efetuada a análise de variáncia e teste de Tukey a 5% de probabilidade para 
todos os testes realizados. A A. aurículiformis (23,61cm-15,90m) superou a O. paraensis (7,26 cm- 
7,06 m) em diámetro a altura, respectivamente. As espécies mostraram densidade média, onde A. 
aurículiformis (0,63 g/cm 3 ) superou a O. paraensis (0,55 g/cm 3 ). O potencial calorífico encontrado 
para A aurículiformis (4383,65 kcal/kg) foi superior estatisticamente ao encontrado na O. paraensis 
(4381,24 kcal/kg). Na análise ¡mediata a A. aurículiformis superou a O. paraensis, mas foi inferior em 
carbono fixo. Os resultados conduzem ao emprego das duas espécies para plantío por apresentarem- 
se promissoras para o cultivo e uso final como fonte energética proposto no trabalho. 

PALAVRAS-CHAVE: 

Análise ¡mediata, biomassa, densidade básica, plantío, poder calorífico. 


ABSTRACT 

At the municipality of Iranduba, ¡n the State of Amazonas, there is a high consumption of wood 
by the local communities, for both charcoal production and to burn for producing bricks. As a result, 
EMBRAPA has developed a project to grow trees at Iranduba’s communities. So, the goal of this study 
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was to evalúate the species Acacia auri- 
culiformis (exotic) and Ormosia 
paraensis (native) by dendrometric 
variables, basic density, heating power 
and immediate analysis, through an 
experiment at Estagáo Experimental da 
EMBRAPA, in Iranduba. The experi- 
mental design was by random blocks 
with three repetitions, in a total of six 
pareéis with 25 plants each. Five trees 
were taken out of the nine middle trees 
of the useful area. Dises 5 cm thick were 
taken from those trees and turned into 
proof bodies In order to determine the 
basic density and other tests. Analysis of 
variance and Tuckey tests, at 5% proba- 
bility, were performed. A. aurículiformis 
(23,61 cm -15,90 m) was larger than O. 
paraensis (7,26 cm-7,06m) in diameter 
and height, respectively. Both species 
presented médium density, but A. auricu- 
liformis (0,63 g/cm 3 ) had greater valúes 
than O. paraensis (0,55 g/cm 3 ). The 
heating potential of A. auriculiformis 
(4383,65 kcal/kg) was statistically higher 
than that of O. paraensis (4381,24 
kcal/kg). For the immediate analysis A. 
auriculiformis showed greater valúes 
than O. paraensis, but presented lower 
valúes for fixed carbón. The results of 
this study showed that both species can 
be grown and used as proposed. 

KEY-WORDS: 

Immediate analysis, Biomass, Basic 
density, Planting, Fleating power. 


INTRODUQÁO 

A madeira oriunda de florestas plan- 
tadas, por suas características próprias e 
pelas condigoes ecológicas altamente 
favoráveis, tem sido apontada como urna 
opgáo lógica a ser utilizada devido a sua 
potencialidade como fonte energética. 
Porém, de toda lenha consumida, a 
grande maioria é proveniente de florestas 
nativas, estimando um corte anual dessas 
florestas de 500 mil hectares, sendo 
necessário reflorestarem grandes áreas 
para garantir o abastecimento futuro de 
madeira como energético. 

De acordo com Shumacher (2003), 
no Brasil a implantagáo de macigos 
florestais na grande maioria formada por 


espécies exóticas, sao conseqüéncias da 
evolugáo de toda urna estrutura industrial. 
Esta tem como objetivo atender a 
demanda das regióes mais desenvolvidas 
do país com matéria-prima necessária 
para produgáo de papel, celulose, 
chapas, aglomerados, carváo vegetal, 
movéis e outros. 

Para Salomáo (1993), o conceito de 
planta exótica sob a limitagáo geopolítica 
está sendo cada vez mais substituido 
pelo de adaptabilidade da planta á regiáo 
onde se pretende cultivá-la, sendo o 
clima e o solo os aspectos mais determi- 
nantes. 

A madeira destinada para fins ener- 
géticos deve basear-se, entre outros, no 
conhecimento do seu poder calorífico e 
no seu potencial para produgáo de 
biomassa, portanto, é necessário haver 
urna avaliagáo de seus constituintes 
químicos e urna análise de seu potencial 
energético (Vale et al., 2000). Um dos 
parámetros para avaliagáo da qualidade 
da madeira é representado pela sua 
densidade básica, que está ligada direta- 
mente com as propriedades físicas, 
mecánicas e anatómicas, podendo assim 
caracterizar o uso final da madeira 
(Rezende et al., 1998). 

O poder calorífico é um parámetro 
que mede a eficiencia energética e serve 
como importante indicador para conhecer 
a capacidade calorífica de urna determi- 
nada espécie (Oliveira, 1982). O teor de 
umidade é um fator que influencia o valor 
do poder calorífico, quanto menor, maior 
será a produgáo de calor por unidade de 
massa. Farinhaque (1981) indica que, 
para queima a madeira náo pode apre- 
sentar teor de umidade superior a 25%. 

A análise ¡mediata é um teste que 
fornece a percentagem de umidade, 
material volátil, carbono fixo e cínzas. Em 
outras palavras, ela fornece a percen- 
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tagem do material que se queima no 
estado gasoso (material volátil) e no 
estado sólido (carbono fixo), bem como 
dá urna indicagáo do material residual 
(cinzas). Aróla (1976) afirma que quando 
a lenha é queimada, geralmente 75 a 
80% sao materiais voláteis, 20 a 24% é 
carbono fixo e 1 a 3% cinzas. 

Ñas zonas rurais a floresta é funda- 
mental como produtora de lenha para o 
consumo diário. Mais de 80% da madeira 
aproveitada é queimada para cozinhar 
alimentos, aquecer casas e abastecer as 
industrias rurais. 

No municipio de Iranduba no Estado 
do Amazonas há um grande consumo de 
lenha por duas partes distintas. A primeira 
refere-se ás pequeñas comunidades 
dessa regiáo. Estas familias sobrevivem 
desta matéria-prima para fabricagáo de 
carváo vegetal que é vendido ñas locali- 
dades próximas e para a cidade de 
Manaus. Porém, estas comunidades 
enfrentam problemas, como a busca da 
matéria-prima que se torna cada vez mais 
escassa, devido á extragáo predatoria e 
falta da reposigáo da cobertura florestal, 
pois a medida que vai se esgotando, mais 
longe váo ficando as florestas para coleta 
da madeira para fabricagáo de carváo. A 
segunda parte é composta pelas otarias. 
A existencia destas se deve a grande 
concentrado de argila ocorrente no solo 
da regiáo, favorecendo a produgáo de 
tijolos e telhas. 

Visando amenizar os impactos 
causados, a Empresa Brasileira de 
Pesquisa Agropecuária-EMBRAPA vem 
desenvolvendo por meio do experimento 
“Selegáo de espécies florestais nativas e 
exóticas com potencial para produgáo de 
energía na regiáo de Iranduba/AM" que 
objetiva a implementagáo de plantíos que 
possam atender a necessidade do 
consumo de lenha de pequeñas comuni- 
dades no municipio de Iranduba podendo 
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ser usufruído por um período de médio 
em curto prazo. Nesse experimento foram 
avaliados os aspectos económicos e 
técnicos utilizando a biomassa de planta- 
góes de madeiras em áreas degradadas 
trabalhando-se com espécies exóticas e 
nativas de crescimento inicial rápido, 
sendo este um fator importante, urna vez 
que para produgáo de lenha o diámetro é 
parámetro de avaliagáo para qualificagáo 
do produto, pois para este fim pode-se 
considerar o diámetro mínimo de 5 cm. 

Das espécies escolhidas para o 
experimento citam-se a A. auriculiformis 
A. Cunn. Ex Benth, da familia Mimosa- 
ceae. É urna espécie arbórea de cresci- 
mento inicial rápido, originária do sul de 
Papua, Nova Guiñé e norte da Austrália. 
É encontrada em diversos tipos de solo, 
tolerando ampia faixa de pH, além de 
possuir capacidade de sobrevivencia em 
solos pobres e em regióes onde a 
estagáo seca é curta. 

A espécie Ormosia paraensis Ducke, 
pertencente á familia Fabaceae, conhecida 
vulgarmente como tento-agaí, que segundo 
Loureiro et al., (1997), habita regióes de 
térra firme e ocorre na Amazonia Legal, 
Estado do Mato Grosso e parte da Vene- 
zuela (Cordilheira Paracaima). 

Nesse contexto o presente trabalho 
teve como objetivo avaliar as espécies 
lenhosas e mensurar o seu potencial 
energético a fim de programar plantíos 
com as espécies selecionadas em comu- 
nidades do municipio de Iranduba. 

METODOLOGIA 

Localizagáo e descrigáo do experi- 
mento em campo 

Os trabalhos de campo foram 
conduzidos na “Estagáo Experimental do 
Caldeiráo” de propriedade da embrapa 
Amazonia Ocidental, localizada no muni- 
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cípio de Iranduba, no estado do 
Amazonas, abrangendo urna área abso- 
luta de 2.354km 2 . 

O municipio está compreendido 
entre as latitudes 3°00’ e 4°00’S e longi- 
tude 60°00’ e 61°00’W Grw, situando-se 
no Médio Amazonas. É constituido em 
grande parte de terrenos quartenários 
recentes, cujos solos sao colmatados 
anualmente pelas endientes dos grandes 
ríos (Brasil, 1979). 

O clima da área é do tipo Afi, perten- 
cente ao grupo de clima tropical chuvoso 
de acordo com a classificagáo climatoló- 
gica de Kóppen. A temperatura para o 
mes mais frió nunca é inferior a 18°C. A 
distribuigáo das chuvas durante o ano 
indica isotermía, com pluviosidade em 
torno de 2.000mm, nao há propriamente 
veráo e invernó distintos (Ribeiro, 1976). 

O experimento foi implantado no 
mes de janeiro de 1995, onde foram 
testadas as espécies (tratamentos) com 
delineamento em blocos ao acaso, com 
trés repetigóes perfazendo 6 parcelas 
compostas de 25 plantas cada, no espa- 
gamento de 4m x 4m. Os blocos I e II 
foram instalados em terrenos degra- 
dados, compactados, anteriormente utili- 
zados com culturas de ciclo curto. O 
bloco III foi instalado ¡mediatamente após 
o corte de urna área de floresta secun- 
dária de aproximadamente 20 anos. 

Amostragem e coleta de dados 

Durante todo o período do experi- 
mento foram realizadas as medigóes dos 
elementos dendrométricos, altura e 
diámetro (dap) executados com o auxilio 
de urna vara métrica e urna fita diamé- 
trica, respectivamente, para acompanhar 
o desenvolvimento das espécies em 
todas as árvores das parcelas. 


Em 2002 foi realizada a derrubada 
das árvores para retirada de amostras 
para os testes experimentáis. De cada 
parcela, eliminou-se a bordadura, e sele- 
cionou-se da área útil, aleatoriamente 
cinco árvores dentre ás nove centráis. 
Após a derrubada, as árvores seccio- 
nadas foram marcadas e délas retirados 
discos com 5 cm de espessura, locali- 
zados em diferentes alturas do fuste 
(base, DAP, 25%, 50%, 75% e 100%), 
tendo como parámetro o diámetro mínimo 
de 5cm na ponta mais fina. De cada 
espécie obtiveram-se trinta discos, totali- 
zando sessenta. Estes foram desdo- 
brados e confeccionados os corpos de 
prova para o teste da densidade básica. A 
Acacia auriculiformis obteve 715 amos- 
tras confeccionadas e a Ormosia 
Paraensis 420. Essa desigualdade em 
número de amostras deu-se devido ao 
diámetro diferenciado do fuste de cada 
urna. 

Descrigáo dos experimentos em labo- 
ratorio 

a) Determinagáo da Densidade básica 
A densidade básica das amostras 
foi determinada pelo método de 
imersáo. Neste, 1,135 corpos de provas 
foram mantidos submersos em água por 
8 dias aproximadamente, até estarem 
acima do ponto de saturagáo das fibras 
para determinagáo do volume verde. 
Para este, utilizou-se urna balanga 
analítica, onde as amostras foram 
submersas em um béquer com 
água, posteriormente foi medido o seu 
volume verde em fungao do peso do 
líquido deslocado pela sua imersáo. As 
amostras foram mensuradas no sentido 
base-topo e medula-casca. O peso seco 
foi obtido após o condicionamento das 
amostras na estufa com circulagáo de ar 
forgada a urna temperatura de 103°C, 
até atingirem peso constante. 
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b) Determinagáo do Poder Calorífico em 
Bomba Calorimétrica 

Os corpos de prova foram transfor- 
mados em cavacos, picados e moídos em 
um moinho, posteriormente classificados 
em peneiras. As fragoes classificadas 
abaixo de 60 mesch foram utilizadas para 
determinagáo do poder calorífico supe- 
riores, segundo a Associagáo Brasileira 
de Normas Técnicas-ABNT-NBR 8633/84. 

Inicialmente pesou-se 1 g da 
amostra seca em balanga analítica. 
Posteriormente, essa amostra foi enviada 
para a bomba calorimétrica que continha 
5 mi de água. Esta foi fechada e carre- 
gada lentamente com oxigénio a urna 
pressáo de 3,0 Mpa, sem deslocar o ar 
inicial. A bomba foi ativada a urna tempe- 
ratura inicial de 27,32°C, sendo acompan- 
hada constantemente por meio de um 
termómetro acoplado a bomba, onde se 
obteve a temperatura final de 28,96°C. 
Esse procedí mentó foi realizado em tres 
repetigóes para cada espécie, visto que 
se tratava da mistura de corpos de prova 
oriundos de árvores diferentes. 

c) Anáiise ¡mediata 

Para determinagáo da anáiise 
¡mediata utiiizaram-se amostras rmoídas 
que foram classificadas em peneiras de 
40 e 60 mesch. A fragáo intermediária foi 
utilizada para realizagáo destes testes. 
Esse procedimento obedeceu as normas 
da abnt nbr 8112/86 que determina o 
método de ensaio do teor de umidade, 
cinzas, materiais voláteis e carbono fixo. 

Com excegáo do carbono fixo, os 
demais testes foram realizados em dupli- 
catas. 

Teor de umidade 

Foram separadas em cápsulas de ago 
previamente taradas, 500 g das amostras 
moídas e pesadas em balanga analítica 
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para obtengáo da massa inicial. Posterior- 
mente foram colocadas em estufa previa- 
mente aquecida, onde permaneceram por 
um período de 48 horas a urna temperatura 
de 1 05 °C. Após esse periodo, as amostras 
foram retiradas e novamente pesadas para 
obtengáo da massa final. 

Teor de cinzas 

Em cadinhos de porcelana previa- 
mente tarados em balanga analítica, 
depositou-se um grama de cada amostra 
seca. Posteriormente, as amostras 
contidas nos cadinhos sem tampa foram 
colocadas na mufla pré-aquecida, por um 
período de 1h a temperatura de 700 °C. 
Após este tempo foram retirados e depo- 
sitados em dessecador para esfriar, e 
novamente pesados para obtengáo da 
massa final. 

Materiais voláteis 

Em urna balanga analítica, pesou-se 
um grama da amostra no cadinho previa- 
mente tarado. Após a pesagem, os 
cadinhos foram tampados e colocados 
por 3 minutos sobre a porta da mufla 
previamente aquecida. Em seguida colo- 
cados dentro da mufla por 7 minutos a 
urna temperatura de 900 °C. Também reti- 
rados e colocados em dessecador para 
obtengáo da massa final. 

Carbono fixo 

O teor de carbono fixo, por se tratar 
de urna medida indireta foi calculado de 
acordo com a equagáo abaixo: 

CF = 100 -(CZ + MV) 

Onde: CF = teor de carbono fixo, em 
CZ = teor de cinza, em %. 

MV = teor de matérias voláteis, em % 
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Análise dos Dados 

A partir dos dados obtidos utili- 
zaram-se os recursos do programa esta- 
tístico sas, para realizar as análises de 
variáncia (anova). Os resultados das 
médias foram contrastados pelo teste de 
Tukey (t), considerando a significancia de 
5% de probabilidade. 

RESULTADOS E DISCUSSÁO 
Variáveis dendrométricas 

O primeiro fator nos resultados que 
coopera para implementagáo dessas 
espécies é o rápido crescimento. Esse 
comportamento foi observado durante os 
sete anos ñas parcelas. Nesse período 
houve acompanhamento do incremento 
da altura e diámetro. Os dados coletados 
foram organizados em planilha e calcu- 
lados para obtengáo de sua média, como 
verificado na tabela 1 . 

Observa-se na tabela 1 que há dife- 
renga significativa ñas características 
dendrométricas das duas espécies estu- 
dadas. Nota-se que a espécie exótica A. 
auriculiformis destacou-se em diámetro e 
altura, enquanto que a espécie nativa O. 
paraensis apresentou dados estatísticos 
inferiores. Ressalta-se aínda a impor- 
táncia da adaptabilidade e desenvolvi- 
miento da espécie exótica no plantío. 

Depreende-se que mesmo com dife- 
renga significativas de dap e a altura, 


ambas as espécies sao promissoras para 
o uso da lenha e/ou fabricagáo do carváo 
vegetal. Urna vez que o diámetro mínimo 
aceito para queima é de 5 cm, nesse 
sentido a O. paraensis enquadra-se nos 
requisitos para tal finalidade. 

Outro fator que contribuiu para o 
desenvolvimento das espécies utilizadas 
neste estudo, foi o espagamento. Com a 
dimensáo 4x4 adotada nao ocorreria 
maior competigáo entre as árvores das 
parcelas, também observado por Brasil & 
Ferreira (1971), quando afirmam que o 
espagamento tem influencia marcante na 
produgáo e qualidade da madeira, princi- 
palmente no manejo de povoamentos 
cujo objetivo é a produgáo de madeira 
industrial de rápido crescimento e densi- 
dade adequada. Rensi Coelho et al. 
(1970), em estudo de diferentes espaga- 
mentos para Eucalyptus propinqua e 
Eucalyptus granáis, concluiram que a 
menor competigáo entre as árvores 
resulta em menor porcentagem de domi- 
nadas e conseqüente maior uniformidade 
do povoamento, com reflexos favoráveis 
ao volume aproveitável final. 

Densidade básica 

Os dados obtidos da densidade 
básica ao longo do fuste, ñas alturas 
propostas, apresentaram diferengas 
pouco significativas, conforme variagáo 
observada na figura 1 , pois o comporta- 
mento das espécies foi aproximadamente 
homogéneo. 


Tabela 1. Valores médios dendrométricos das espécies estudadas aos sete anos. 


ESPÉCIE 

DAP (CM) 

ALTURA (M) 

Acacia auriculiformis 

23,61 

15,90 

Ormosia paraensis 

7,26 

7,06 
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Figura 1. Variagáo dos valores médios ao longo do fuste das espécies. 


Observa-se que em ambas as espé- 
cies a densidade decaí da base para o 
DAP, e a partir da posigáo de 25% ela 
sofre oscilagóes até 100%, sendo que a 
A. auriculiformis segue urna linha cons- 
tante, ao contrário da O. paraensis que 
apresenta variagoes mais evidentes. 
Contudo, essa diferenga de comporta- 
miento entre as espécies é comprovada 
em urna variagáo de 0,48 a 0,69 (g/cm3), 
classificando estas espécies como 
madeiras de densidade média. Nesse 
sentido, o resultado infere o uso da totali- 
dade do fuste da árvore para queima, 
urna vez que o menor valor obtido foi de 


0,48g/crri3 na altura de 100%. Ressalta- 
se que nao haverá necessidade de sele- 
cionar partes do fuste da árvore para 
consumo na hora da queima, evitando 
assim o desperdicio de residuos flores- 
tais, como encontrado muitas vezes em 
trilhas ou beiras de estradas, restos de 
fuste e galhos por apresentarem diámetro 
pequeño. 

Observa-se na tabela 2 que ocorreu 
urna variagáo significativa na média entre 
as duas espécies. Em comparagáo com o 
gráfico apresentado, essa diferenga nao 
se apresentava táo discrepante, como 


Tabela 2. Valores médios da densidade básica das espécies. 


ESPÉCIE 

DENSIDADE BÁSICA 


(g/cm 3 ) 

Acacia auriculiformis 

0,63 

Ormosia paraensis 

0,55 
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constatado no quadro de média geral. 
Contudo, essa variagáo nao é nociva, 
pois as médias encontradas enquadram- 
se nos valores de densidade média para 
madeira, inferindo a qualidade destas 
para utilizagáo final. 

A espécie A. auriculiformis mostrou- 
se com maior densidade (0,63 g/cm 2 ) 
quando comparada a O. paraensis (0,55 
g/cm?). Resultado semelhante quanto á 
média da densidade básica de espécies 
exóticas na regiáo brasileira foi encon- 
trado por Araújo (2000), quando estudou 
a introdugáo da espécie exótica Nim 
indiano (0,57 g/cm 2 ) na regiáo semi-árida 
do Nordeste. Vale (2001) encontrou 
médias similares estudando a densidade 
básica de algumas espécies do cerrado 
para correlagáo entre a densidade e a 
qualidade do carváo, citam-se a Murta 
(0,52 g/cm 2 ), Pequi (0,62 g/cm 2 ) e Jaca- 
randá-do-cerrado (0,77 g/cm 2 ). 

Com base nos valores encontrados 
pode-se afirmar que a A. auriculiformis e 
O. paraensis sao espécies com densi- 


dade equivalente para a geragáo de 
energía na forma de calor, aínda que de 
origens diferentes, responderam bem 
com a adaptagáo ao solo e crescimento. 

Poder calorífico 

Na figura 2 observa-se a média da 
variagáo do poder calorífico superior 
disponibilizado para as espécies anali- 
sadas, sendo possível caracterizar urna 
pequeña diferenga entre elas. 

Os dados médios encontrados para 
A. auriculiformis (4383,65 kcal/kg) e O. 
paraensis (4381,24 kcal/kg) tornam-se 
satisfatório, devido ao período de cultivo 
e condigóes locáis. Mesmo sendo espé- 
cies de origem diferentes alcangaram um 
valor médio de eficiencia energética que 
satisfaz duas faixas de estudos, urna 
realizada por Howard (1973) afirmando 
que o poder calorífico superior para 
folhosas varia na faixa de 4.600 a 4.800 
kcal/kg, enquanto que, para Brito (1993), 
este valor para folhosas trapicáis está 
entre 3.500 a 5.000 kcal/kg. 
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Figura 2. Valores médios do poder calorífico para as espécies. 
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Em estudos realizados, Araújo 
(2000) encontrou um potencial de 
4088,50 kcal/kg para espécie exótica Nim 
indiano. Jara (1989) encontrou 4.550 
kcal/kg para Acacia decurrens (Acacia 
negra) aos cinco anos de idade e 4.790 
kgcal/kg para o Eucalyptus grandis aos 
dez anos de idade. 

Observa-se que os valores médios 
encontrados aos sete anos sao distintos 
aos das espécies com cinco e dez anos. 
Nesse contexto pode-se afirmar que 
mesmo como experimento para implan- 
tagáo as espécies corresponderam bem 
ao potencial energético, e poderáo ser 
aprimoradas com mais técnicas e trata- 
mientos silviculturais. 

Análise Imediata 

Observa-se na tabela 3 que os 
testes realizados na análise imediata 
apontaram urna diferenga significativa 
entre as duas espécies analisadas. 

A espécie A. auriculiformis apre- 
sentou baixa porcentagem de teor de 
umidade, cinzas e carbono fixo, e indicou 
alto valor para materiais voláteis em 
comparagáo a O. paraensis. A perda de 
umidade para igual período exposto em 
temperatura utilizada foi maior para, A. 
auriculiformis. Essa característica é 
importante, urna vez que a quantidade de 
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água na madeira influencia no calor libe- 
rado, o que em processo industrial 
ocasiona redugáo na produtividade, por 
isso a escolha da madeira que possua 
perda de umidade rápida é essencial no 
processo produtivo. Essa característica é 
encontrada ñas espécies estudadas. 
Apesar da O. paraensis ter obtido um 
valor estatístico um pouco elevado ao da 
A. auriculiformis, ela também corres- 
ponde bem ao resultado, já que seu valor 
no teor de umidade foi 8,98%, este se 
encontra na faixa para liberagáo de calor, 
proposto em estudos realizados por 
Farinhaque (1981), quando afirma que o 
teor de umidade para queima nao pode 
ultra passar 25%, 

Quanto ao Teor de Cinzas A. auricu- 
liformis apresentou urna porcentagem de 
0,59% correspondendo a bom aproveita- 
mento da matéria-prima, enquanto que 
para Ormosia paraensis foi 1,35%. Os 
dados encontrados neste estudo estáo 
condizentes com os realizados sobre a 
porcentagem de cinzas, como por 
exemplo, Araújo (2000) que trabalhando a 
espécie Nim indiano detectou 2,11% de 
cinzas. Marabotto et al., (1991) também 
encontraram resultados similares em um 
estudo com 28 espécies da Amazonia 
Peruana-Brasileira que sofreram varia- 
góes como a Moena amarilla (0,11%), 
Osteophloeum plastyspermum (0,28%), 
Loro shungo (0,58%), Caucho masha 


Tabela 3. Valores médios da análise ¡mediata para as espécies estudadas. 




VARIÁVEIS 


ESPÉCIES TEOR DE UMIDADE 

(%) 

CINZAS 

(%) 

MATERIAL VOLÁTIL 
(%) 

CARBONO FIXO 

(%) 

Acacia auriculiformis 

7,59 

0,59 


15,41 

Ormosia paraensis 

8,98 

1,35 

79,17 

19,48 
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(0,76%), Virola albidiflora (1,04%), 
Anacaspi (1,14%), Bombax munguba 
(2,32%) e Maquira coreacea (4,91%). 
Para o mesmo autor, o baixo teor de 
cinzas representa o mínimo de impurezas 
que pode ser encontrado em madeiras 
para combustíveis. 

Para o Teor de Materiais Voláteis, a 
A. auriculiformis (84,00%) apresentou 
média superior a O. paraensis (79,17%) 
possuindo assim maior facilidade de se 
queimar ou liberar calor (gases). Correla- 
cionado com a quantidade de carbono 
fixo, a A. auriculiformis (15,41) obteve 
maior porcentagem que a O. paraensis 
(19,48). Neste sentido pode-se afirmar 
que A. auriculiformis libera maior quanti- 
dade de material gasoso, deixando um 
bom percentual de material sólido 
(carbono fixo) para queima. 

CONCLUSÁO 

Com base nos resultados obtidos, 
conclui-se que: 

Os valores da densidade básica 
encontrados ñas duas espécies corres- 
pondem á classificagáo de espécies 
madeireiras de densidade média; 

O potencial energético da Acacia 
auriculiformis (4383,65 Kcal/kg) foi supe- 
rior a Ormosia paraensis (4381,24 
Kcal/kg) com pequeña diferenga estatís- 
tica; 

Os valores da análise ¡mediata da 
Ormosia paraensis na variável teor de 
umidade, cinzas e carbono fixo foram 
superiores ao da Acacia auriculiformis ; 

O teor de material volátil da Acacia 
auriculiformis foi superior a Ormosia 
paraensis; 


Para um programa de plantío comer- 
cial e/ou agroflorestal com o objetivo de 
produzir biomassa para a geragáo de 
energía, as duas espécies apresentaram- 
se bastante promissoras e grande poten- 
cial para cultivo. 

As especies A. auriculiformis e 
O. paraensis ,de acordo com este estudo 
podem ser usadas em plantíos de rege- 
nerado de áreas degradadas e fins 
madeireiros, como construgáo e movéis. 
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Anatomía, índices físicos e hidráulicos de la 
madera de Gliricidia sepium (Jacq.) Steud. 

Anatomy, hydraulic and physical indexes 
of Gliricidia sepium (Jacq.) Steud. wood 

Alejandra Quintanar-lsaías 1 - 3 , Guillermo Ángeles 2 
y José Alejandro Zavala-Hurtado 1 

RESUMEN 


Se presentan los resultados de un estudio de la variación anatómica y de los índices físicos e 
hidráulicos de la madera de albura y duramen de seis árboles de Gliricidia sepium (Jacq.) Steud. 
colectados en épocas de secas y lluvias de una selva baja de Veracruz, México. De rodajas cortadas 
a 1 ,30 m del suelo, se hicieron probetas de diversos tamaños para estudiar el contenido de humedad, 
la densidad relativa y la anatomía macroscópica y microscópica. Se usaron cubos con los que se 
elaboraron preparaciones fijas de cortes y material disociado, siguiendo los protocolos clásicos de 
laboratorio, así como las clasificaciones y nomenclaturas internacionales para los elementos celu- 
lares. Los cálculos de índices de vulnerabilidad, de resistencia al colapso y K^eórica se hicieron 
usando el diámetro, la densidad y el área efectiva de vasos. Los análisis de varianza y de compo- 
nentes principales mostraron que las diferencias entre árboles, zonas y estaciones están dadas por 
algunas características de las fibras, los vasos, los radios y las punteaduras. Aunque la variación en 
estas características es baja, explicaría la estabilidad dimensional de la albura y el duramen a los 
cambios de régimen hídrico ocasionados por la sequía y la lluvia. Los índices de resistencia al colapso 
y de vulnerabilidad, así como la K/, muestran un xílema mesomórfico muy vulnerable a embolismos, 
pero con mucho éxito para desarrollarse en ambientes perturbados, lo que podría estar explicado por 
la presencia de tilides, abundante parénquima axial y valores altos de la densidad relativa que 
aseguran resistencia de las paredes de los vasos y fibras al colapso. 

PALABRAS CLAVE: 

Anatomía de la madera, contenido de humedad, densidad relativa, índices de vulnerabilidad y colapso 

y ^h-teórica- 


ABSTRACT 

An anatomical, physical and physiological study of the heartwood and sapwood of Gliricidia sepium 
(Jacq.) Steud. was carried out using six trees collected during the rainy and dry seasons in a dry forest from 
Veracruz, México. Stem segments were cut at 1,30 m above the ground for obtaining wood samples for 
macroscopic and microscopio observations, and for the water contení and specific gravity analyses. Follo- 
wing classical microtech ñique procedures and ¡nternational classifications and nomenclatures for descríbing 
microscopio features, wood samples and permanent slides of wood sections and dissociated material were 
made. Estimations of vulnerability and collapse indexes, and K^heoreticai were made from vessel diameter, 
vessel number and effective vessel area Analysis of variance and principal components analysis showed 
that some fiber, vessels, rays and pit features influenced differences among trees, sites and seasons. 
Although variability of those features is relatively low, it would explain the sapwood and heartwood stability 
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during changes in the hydric regime 
caused by dryness and rain. Collapse and 
vulnerability indexes, as well as K/, show a 
mesomorphic xylem quite vulnerable to 
embolisms, but with high success to 
develop in disturbed environments. This 
could be explained by the presence of 
tyloses, abundant axial parenchyma and 
high valúes of specific gravity that assure 
resistance of fiber and vessel cell walls to 
collapse. 

KEY WORDS: 

Wood anatomy, moisture content, 
specific gravity, vulnerability and 
collapse indices y K h . te0 reth¡cai- 


INTRODUCCIÓN 

Gliricidia sepium ha sido una 
especie muy usada para leña, postería y 
artesanía; se le ha designado como una 
especie multipropósito por su empleo en 
agrocultivos, como árbol de sombra para 
el cacao, en barbechos, callejones, 
cortinas rompevientos, pero, sobre todo 
como cercas y postes vivos (Quintanar et 
al., 1997; Vázquez-Yánes et al., 1999). Se 
ha visto que sus hojas pueden ser tóxicas 
o con características alelopáticas, pero se 
usa como forraje ya que tiene altos rendi- 
mientos de biomasa (Elevitch y Francis, 
2006). En algunos países se prepara una 
pasta de hojas y corteza que es usada 
como raticida o insecticida. Es una buena 
fijadora de nitrógeno y tiene buena rege- 
neración después de incendios, los que 
aparentemente favorecen su estableci- 
miento. Su hábitat natural es semihú- 
medo, pero está ampliamente distribuida 
en México y Centroamérica, y dada su 
facilidad para establecerse y el uso exten- 
sivo que se le ha dado, se puede encon- 
trar en dunas costeras, bancos ribereños, 
planicies inundables, faldas de montañas, 
barrancos, terrenos abiertos, terrenos 
inestables y áreas perturbadas. 

Asimismo, se considera una especie 
potencialmente restauradora y reforesta- 
dora de ambientes perturbados 

(Vázquez-Yánes et al., 1999). Reciente- 


mente se ha propuesto para sustituir 
especies comerciales duras que están 
escaseando en el mercado internacional 
(Oluwafemi y Adegbenga, 2007). 

Se conoce la anatomía de la madera 
de aproximadamente 25 especies mexi- 
canas de la subfamilia Papilionoideae, a 
la que pertenece G. sepium (Cárdenas, 
1971; Barajas et al., 1979; Barajas, 1980; 
De la Paz Pérez et al., 1980; De la Paz 
Pérez y Corral, 1980; Corral, 1985; Rebo- 
llar et al., 1987; Barajas y León, 1989; 
López y Ortega, 1989; Chehaibar y 
Grether, 1990; Ortega et al., 1991; 
Santiago y Ortega, 1992; De la Paz 
Pérez, 1993; Rebollar et al., 1996-1997; 
Rebollar y Quintanar, 1 998 y Razo, 2003). 
Estos estudios destacan patrones celu- 
lares generales presentes en Gliricidia 
sepium. 

Carlquist (1977) estudió el impacto 
del diámetro, la longitud y la morfología 
de los elementos de vaso y su relación 
con algunos factores ecológicos, princi- 
palmente la precipitación media anual en 
diferentes comunidades vegetales. Ese 
autor analiza los valores de los índices de 
vulnerabilidad (I. V.) y mesomorfia (M) 
para especies que crecen en distintos 
hábitats. El autor sostiene que las espe- 
cies cuyo xilema presenta valores de I. V. 
por arriba de 1 son vulnerables o poco 
resistentes al estrés hídrico, y cuando los 
valores son menores a 1 entonces están 
bajo estrés hídrico, es decir, son resis- 
tentes al embolismo. Agrega que el 
producto del índice de vulnerabilidad y la 
longitud de los vasos es el valor hídrico 
mesomórfico de respuesta del xilema; de 
manera que por arriba de 200 el xilema se 
considera mesófito y por debajo de 75 
xerófito. En este sentido el xilema de G. 
sepium se consideraría como mesófito en 
condiciones de estrés hídrico, pues de 
acuerdo con Rebollar y Quintanar (1998), 
individuos de esta especie en Quintana 
Roo presentaron pocos vasos/mm 2 
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(4/mm 2 ) y un diámetro tangencial 
mediano (130 [75-232] jim). Asociado a 
estas características anatómicas, la 
densidad relativa de la madera y las 
dimensiones de las fibras son caracterís- 
ticas que se han ligado directamente con 
el colapso o implosión de los vasos, por 
las presiones negativas causadas por la 
evapotranspiración (Hacke et al., 2001; 
Jacobsen et al., 2005). De manera que la 
prevención de la implosión de un vaso 
depende de la relación entre el grosor de 
la doble pared de vasos vecinos y el 
diámetro del vaso (t/b) 2 y con el área de 
la pared de la fibra (Baas et al., 2004). La 
densidad del xilema secundario es una 
propiedad relacionada con las capaci- 
dades hidráulicas del tallo y la resistencia 
mecánica de la planta, es decir, entre más 
seguro un xilema, mayor será su 
densidad (Swenson y Enquist, 2007). 

Quintanar et al. (1997, 2003) estu- 
diaron algunas características histoquí- 
micas de albura y duramen de G. sepium 
que muestran los diferentes tipos de 
tílides, su abundancia y su desarrollo in 
vitro , el tipo de extractivos, cristales, 
algunos tipos de pectinas y procesos de 
esclerifícación en células de ambas 
zonas. Oluwafemi y Adegbenga (2007) 
estudiaron las características tecnoló- 
gicas de la madera de G. sepium en 
árboles del Campus de la Universidad de 
Ibadan en Nigeria, y han señalado que los 
extractivos encontrados en esta especie 
están relacionados con la alta densidad 
de esta madera y su posible resistencia al 
deterioro, por lo que proponen su empleo 
en exteriores. 

Gliricidia sepium (cocoíte) puede ser 
un árbol con tallas desde 5 m de altura y 
15 cm de diámetro a la altura del pecho 
(d.a.p.), hasta 20-25 m de altura y diáme- 
tros mayores a 40 cm de d.a.p., depen- 
diendo del tipo de vegetación al que 
pertenece. Forma ramas gruesas ligera- 
mente horizontales, pero las de nuevo 


crecimiento son verticales (Pennington y 
Sarukhán, 2005 y observaciones de 
campo de los autores). Esta especie está 
ampliamente distribuida en Mesoamérica 
y ha sido introducida a África e Indonesia 
(Elevitch y Francis, 2006). Los árboles de 
esta especie florecen y pierden sus hojas 
en época de sequía, mientras que en 
época de lluvia producen abundante 
follaje. El tronco puede torcerse, pero en 
general es recto en la mayoría de los indi- 
viduos y la corteza del tronco principal es 
ligeramente escamosa, en las ramas 
jóvenes es verde, lisa, con lenticelas 
verticales y de forma oval. 

OBJETIVOS 

El propósito de este trabajo es 
contribuir al conocimiento de la anatomía 
de la madera de Gliricidia sepium, a 
través del estudio de la variación celular, 
el índice de vulnerabilidad (I.V.), el conte- 
nido de humedad (C.H.), la densidad rela- 
tiva (P. E.) y valores hidráulicos como: la 
conductividad hidráulica (K| >t eór¡ca de 
vasos y por área) y los índices de resis- 
tencia al colapso o implosión de fibras y 
vasos (I. C. fibras e C. vasos ) entre albura 
y duramen en las épocas de secas y 
lluvias, con el fin de entender algunas 
estrategias del xilema de esta especie. 

METODOLOGÍA 

Sitio de muestreo. El sitio de colecta 
se encuentra aproximadamente a 2,6 km 
de la playa y a 2.8 km del río La Antigua, 
ubicada al Suroeste de la Laguna La Cata- 
lana, carretera Cardel-Veracruz, a una 
latitud de 19° 16’ 51.6” al Norte y longitud 
de 96° 17’ 20.28” al Oeste; con una eleva- 
ción de 76 msnm. En este sitio domina una 
selva baja caducifolia en condición pertur- 
bada por el pastoreo extensivo, la extrac- 
ción de madera para leña, postería y la loti- 
ficación. Además se encuentran legumi- 
nosas como Acacia y Mimosa, entre otros 
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géneros (Rzedowski, 1988). El suelo es de 
tipo regosol calcárico, con fase textural 
gruesa. La unidad geológica está formada 
por rocas sedimentarias y volcano sedi- 
mentarias debido al acarreo y retrabajo de 
arenas litorales por acción eólica (antiguas 
dunas fijas). Las arenas son de grano 
grueso, medio a fino, compuestas por 
feldespatos, micas, fragmentos de roca, 
cuarzo y conchas (inegi. Carta geológica y 
edafológica, s/a). Los promedios de preci- 
pitación y temperatura registrados para los 
meses de enero y febrero son: 30,5 mm y 
21,7° C; y 383 mm y 27,8° C respectiva- 
mente (Figura 1) (Mehltreter, 2007). 

Procesamiento del material 

En el sitio de muestreo se colectaron 
los fustes rectos de tres árboles maduros 
de 7 m de altura promedio, en el mes de 
enero (época de secas) y tres en el mes de 
agosto (época de lluvias), con las mismas 
características. Los árboles tuvieron de 20 
a 25 cm de ancho a la altura del pecho. Se 
cortaron rodajas transversales del tronco 
principal de 2 cm de grosor a 1,30 m del 
suelo para los siguientes estudios: 


1. Estudio anatómico 

Para el estudio anatómico microscó- 
pico se obtuvieron cubos de 1 cm 3 , de 
albura y de duramen. Cada cubo se 
ablandó usando una solución de alcohol- 
glicerol-agua (a.g.a.) infiltrado al vacío 
durante una semana. El a.g.a. se removió 
infiltrando una solución de peg 4000 al 
5% con golpes de vacío, durante varios 
días y se seccionaron a -20°C en un 
criostato marca Leica. Se elaboraron 
cortes transversales, tangenciales y 
radiales; las muestras se pegaron en 
portaobjetos, se tiñeron con azul de tolui- 
dina, se lavaron y se dejaron secar. 
Después se usó xilol y se montaron en 
resina Entellan. El material disociado se 
obtuvo empleando una solución de 
Jeffrey y se montó en gelatina glicerinada 
de acuerdo con Johansen (1940). 

En la sección transversal se 
midieron el diámetro de los vasos y su 
número/mm 2 ; en la tangencial la altura, la 
anchura de radios así como el número de 
radios/mm y la longitud de vasos. En el 
material disociado se midieron la longitud, 
grosor y diámetro de fibras. En la sección 
radial se midieron la altura y la achura de 



Figura 1. Diagrama ombrotérmico en el que se representan los meses más calurosos y 

los más húmedos para la zona de estudio. 
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las punteaduras intervasculares. Se 
hicieron al menos 30 mediciones por 
cada carácter. La clasificación de los 
tamaños se hizo de acuerdo a Chattaway 
(1932). Los caracteres cualitativos se 
describieron de acuerdo con la lista de 
características microscópicas del Comité 
de la iawa (1989). 

Para el estudio macroscópico se 
obtuvieron muestras de madera longitudi- 
nales radiales y tangenciales de 15 x 7 x 
1 cm; el hilo, la textura y el veteado se 
determinaron de acuerdo con Tortorelli 
(1956) y el color se determinó usando las 
Tablas Munsell de color de suelo (Munsell 
Color Co., 1954). 

2. Densidad relativa y contenido de 
humedad 

De las rodajas de los árboles 
cortados se obtuvieron cubos de albura y 
duramen de 2 cm3, con los que se calculó 
la densidad relativa y el contenido de 
humedad. Se obtuvieron los siguientes 
valores: peso del volumen desplazado, 
peso en verde y peso seco final. 

Para el cálculo de los valores se 
emplearon las siguientes ecuaciones: 

Densidad relativa (P.E.) = peso 
anhidro / peso del volumen verde 

Contenido de humedad (C.H.) = 
(peso verde-peso anhidro) / peso anhidro. 

3. Conductividad hidráulica teórica 

Kh -teórica 

Se obtuvo usando los datos del 
diámetro y número de vasos/mm 2 , con 
los que se calculó el área efectiva teórica 
de conducción, que se definió como la 
porción de la superficie de la sección 
transversal que está ocupada por vasos y 
que permite el paso del fluido. La repre- 
sentación de esta característica se asoció 
a una superficie de sección circular, que 
se llamó área efectiva, esto es, 
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Aefectiva (Ae)=(# poros/mm 2 * )*(jtd va so 2 )/4. 
Usando este valor se estimó el diámetro 
equivalente a través del cual ocurre el 
flujo y que está asociado a la sección 
transversal, quedando expresado como 
Deq sec= (4*A e f e ctiva^" 1 ) 0 - 5 - Con este 
diámetro, que equivale a la sección 
conductora, se obtuvo la ecuación de flujo 
que considera el principio de conserva- 
ción de la masa o continuidad y la ecua- 
ción de Darcy, esto es: m f- 1 = 
pjiD4(128»-ic/P d/-i=pjrH(8»-idP dH, 
donde m es la masa del fluido, t el tiempo, 
dP el cambio de la presión, di la longitud 
del segmento de tallo, D el diámetro del 
vaso o el equivalente (D eq S ec)> r el radio 
de esa misma sección, ¡.i la viscosidad 
dinámica del fluido y dh 1 la densidad del 
fluido (Streeter y Wyllie, 1985). 

4. índices de vulnerabilidad y de colapso 

El índice de vulnerabilidad se calculó 
de acuerdo con la ecuación: I.V. = 
(diámetro de vaso) (número de vasos)- 1 y 
el índice de resistencia al colapso o 
implosión: I.C. = (grosor de dos paredes 
contiguas / diámetro celular) 2 (Hacke et 
al., 2001). 

Los valores de K^. vaso , K h . área , los 
IC vaso , ^ fibra Y el l-V., solamente se 
calcularon para los vasos de la albura, 
pues el duramen ha dejado de conducir 
agua. Para comparar los valores de las 
conductividades de vasos y áreas 
(K/ 7 -vaso, K/várea) de G. sepium, se obtu- 
vieron los de otras especies de legumi- 
nosas procedentes de una selva baja 
caducifolia del estado de Oaxaca (datos 
inéditos). La relación K/ 7 _á rea /K/ >vas0) se 
usó para tener un parámetro comparativo 
entre las especies. 

5. Análisis estadístico 

Se estimaron medias, desviación 
estándar y error estándar, con los que se 
obtuvieron las categorías numéricas y la 
clasificación de las características anató- 
micas mensurables, así como los valores 
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medios de los pesos, contenidos de 
humedad, Kh_ {eórica , índices de vulnerabi- 
lidad e implosión, con los que se hicieron 
las gráficas considerando los siguientes 
niveles: entre árbol, entre zonas y entre 
estaciones. Para evaluar las diferencias 
entre características anatómicas micros- 
cópicas, la densidad relativa, el contenido 
de humedad, la conductividad hidráulica 
(K/,), el I. V. y los índices de resistencia al 
colapso (I. C.) de vaso y fibras, se reali- 
zaron análisis de varianza (anova) consi- 
derando la época de colecta, albura y 
duramen y los árboles como variables 
categóricas. Para ambos análisis se usó 
el programa Stata 7 (Stata Corp., 2001). 

Se realizó también un análisis multi- 
variado de las características anatómicas 
microscópicas con el propósito de 
detectar diferencias de los caracteres 
microscópicos mensurables en los 
niveles: 1) entre árboles, 2) entre albura y 
duramen y 3) entre la temporada seca y la 
de lluvias, e identificar a las variables 
causantes de estas diferencias. Se reali- 
zaron análisis discriminantes múltiples 
(adm) sobre una matriz de 396 muestras y 
11 variables anatómicas. El criterio utili- 
zado para evaluar las diferencias entre 
los grupos fue la X de Wilks, que tiende a 
tomar valores cercanos a cero cuando los 
grupos están bien definidos en función de 
las variables consideradas y tiende a uno 
cuando los grupos no están bien defi- 
nidos (Tatsuoka, 1970). 

Se realizaron tres análisis discrimi- 
nantes: adm (i) considerando a los seis 
árboles muestreados como grupos dife- 
rentes, (ii) separando a las muestras por 
zonas de albura y duramen, y (iii) 
evaluando las diferencias entre las mues- 
tras provenientes de la temporada de 
secas y la de lluvias. Se extrajeron las 
funciones discriminantes correspon- 
dientes (cinco en el primer caso y una en 
los otros dos), se evaluaron las diferen- 
cias entre los grupos analizados y se 


identificaron aquellas variables que más 
contribuyen a la separación entre grupos. 
Los adm se realizaron con el paquete 
estadístico spss 12.0 (spss Inc., 2003). 

Para el caso del análisis con solo 
dos grupos (entre albura y duramen y 
entre secas y lluvias), considerando que 
el adm solo extrae una función discrimi- 
nante, se generaron estimadores de 
densidad por kernel (edk) para poder 
observar las diferencias entre los grupos 
a lo largo de la función discriminante 
extraída. Se obtuvieron los edk a partir de 
la función ponderada de Epanechnikov, 
con una amplitud de ventana óptima de 
Haerdle, utilizando el programa edk 2000 
(Salgado-Ugarte, 2000). En el caso del 
análisis con seis grupos esta observación 
se realizó al graficar las dos primeras 
funciones discriminantes. 

Comparación de densidad relativa y 
contenido de humedad 

Usando un análisis de la t de 
Student, se comparó la densidad relativa 
y el contenido de humedad (con transfor- 
mación arco seno) de 112 especímenes 
entre zonas de albura y duramen (n=59 y 
n=53, respectivamente) y entre muestras 
de árboles en la época seca y la tempo- 
rada de lluvias (n=58 y n=54, respectiva- 
mente). El análisis se realizó utilizando el 
programa xlStat-Pro 3.0 (Fahny, 1997). 

RESULTADOS 

I. Anatomía de la madera de Gliricidia 
sepium 

Características macroscópicas 

La madera no presenta zonas de 
crecimiento definidas. Presenta diferen- 
cias marcadas entre albura y duramen. 
La albura presenta vetas de color castaño 
muy pálidas (10 YR 8/4) a amarillas (10 
YR 8/6) y el duramen vetas de color 
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castaño oscuro amarillentas (10 YR 4/4) y 
castaño amarillentas (10 YR 5/4). La 
textura es mediana, el hilo es entrecru- 
zado, el veteado es pronunciado, el olor y 
el sabor son característicos. No existen 
diferencias marcadas entre zonas de 
crecimiento. Los vasos y el parénquima 
axial son visibles a simple vista. El 
duramen presenta autofluorescencia. 

Características microscópicas (tablas 1 y 
2; figuras 2A-H). 

La porosidad es difusa. Los 
elementos de vaso son múltiples radiales 
de dos a tres y solitarios, muy pocos, 2 (1-5), 
extremadamente cortos, 180 (80-280) 
pm, y de diámetro mediano, 179 (68-272) 
pm. Las punteaduras intervasculares son 


ornamentadas y las de campo de cruza- 
miento vaso-radio tienen rebordes 
internos bien distinguibles, similares a las 
intervasculares. Las punteaduras miden 
de alto 4,3 (2-10) pm y de ancho 10,8 
(6-20) pm, por lo que se consideran de 
tamaño grande; se presentan tanto en 
caras radiales como tangenciales. La 
pared del vaso posee punteaduras inter- 
vasculares que presentan expansiones 
tubulares dirigidas hacia una membrana 
central con forma de disco. La pared del 
vaso posee punteaduras ¡ntervasculares 
que presentan expansiones tubulares diri- 
gidas hacia una membrana central con 
forma de disco (figuras 2G, 21); las de 
vaso-parénquima terminan en una sola 
membrana dilatada (Figura 2F). Las cavi- 


Tabla 1. Medias de caracteres microscópicos, densidad relativa, contenidos de 
humedad, índices de resistencia al colapso y vulnerabilidad, K^teórica para vaso y área 

de albura. 


ALBURA 

ÁRBOL 1 

ÁRBOL 2 

ÁRBOL 3 

ÁRBOL 4 

ÁRBOL 5 

ÁRBOL 6 

Vasos*mnr 2 

2,2 

2,1* 

2,1 

2,5* 

1,8 

2,3 

1 Diámetro de vaso 

169* 

195,5 

171,15 

192 

176,24* 

158* 

Longitud de vaso 

155 

166 

193 

174* 

211,6* 

183 

1 Altura de punteadura 

4,4 

3,51 

5,85 

7* 

6* 

4,07 

"Anchura punteadura 

11,9 

9,09* 

10 

13,8* 

11,5 

9,48 

Anchura de radio 

41* 

24* 

26* 

38 

33* 

37* 

Altura de radio 

149* 

187 

180 

178 

170* 

160* 

# de rad¡os*mm-i 

10* 

11* 

10 

8* 

10* 

11 

i Diámetro de fibra 

20,4 

18,9 

18 

19,3 

20,44 

20,29 

1 Grosor de la pared de la fibra 7,5 

8,6* 

4* 

8* 

5* 

6* 

1 Longitud de la fibra 

1283* 

1228 

1130* 

1115* 

1060* 

1163* 

í, v, 

93 

98 

82 

87 

107 

76 

í, C vasos 

0,028 

0,0254 

0,01213 

0,0228 

0,0088 

0,0,15 

1, C fibras 

0,40 

0,38 

0,25 

0,43 

0,0,27 

0,25 

Densidad relativa 

0,81 

0,65 

0,70 

0,762 

0,742 

0,713 

Contenido de humedad (%), 

27 

38 

34 

40 

38 

40 

K h -teórica (vaso) 

2,44 x 10-H 

4,02 x 10-H 

2,56 x 10-H 

4,4 x 10 -11 

2,72 x 10-H 

1,79 x 10-H 

Kft -teórica (área) 

1,30 x 10-10 

2,01 x 10-10 

1,38 x 10-10 

3,38 xlO-10 

1,05 x 10-10 

1,01 x 10-10 


*significativamente diferente con respecto al duramen. 1 pm. 
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Figura 2. Gliricidia sepium. A y C. Cortes transversales de vasos cercanos al cambium 
vascular. D. Corte tangencial albura media. F. Duramen. E. Punteaduras 
vaso-parénquima. G y H. Detalles punteaduras vaso-parénquima. T. Tílides. 

P. Pared de vaso. PA. Parénquima. 






Madera y Bosques 15(2), 2009:71-91 


79 


Tabla 2. Medias de caracteres microscópicos, densidad relativa, contenido de humedad, 
índices de resistencia al colapso y vulnerabilidad, K Me ó r¡ca para vaso 

y área de duramen. 


DURAMEN 

A1-D 

A2-D 

A3-D 

A4-D 

A5-D 

A6-D 

Vasos*mm- 2 

2,3 

2,6* 

2,5 

1,9* 

1,6 

2,2 

1 Diámetro de vaso 

197* 

191,3 

178 

199 

172* 

139* 

1 Longitud de vaso 

172 

162 

209 

189* 

237* 

179 

1 Altura de punteadura 

4,46 

3,57 

4,66 

2,33* 

4,84* 

3,39 

1 Anchura punteadura 

11,38 

12,54* 

10,66 

9,81* 

10,5 

9,82 

1 Anchura de radio 

47* 

29* 

30* 

36 

27* 

28* 

1 Altura de radio 

174* 

175 

177 

169 

212* 

126* 

# de radios*mm- 1 

9* 

10* 

10 

10* 

9* 

11 

1 Diámetro de fibra 

20,7 

16,36 

18,5 

20,24 

19,89 

19,82 

1 Grosor de la pared 







de la fibra 

8 

6* 

7* 

4* 

7* 

4* 

1 Longitud de la fibra 

1291* 

1223 

1114* 

1125* 

1298* 

1026* 

Densidad relativa 

0,87 

0,77 

0,77 

0,73 

0,749 

0,75 

Contenido de humedad 







Í%L 

16,4 

23,9 

21,7 

34,1 

32 

28,2 


*sign¡ficativamente diferente con respecto a la albura. V m - 


dades de los vasos presentan tílides 
comunes en la albura y de tipo esclero- 
sadas en el duramen. Las tílides 
provienen del parénquima axial y del 
radial. 

Las paredes de las tílides en la 
albura son muy delgadas y se engruesan 
hacia el duramen. Las tílides presentes 
en el duramen son de paredes muy 
gruesas esclerosadas con abundantes 
cristales romboidales y gomas que 
ocluyen completamente sus cavidades. 

El parénquima axial es aliforme y 
aliforme confluente de varias células de 
ancho. En caras tangenciales se 
presentan en cordones de más de 8 
células de longitud. Presentan cristales 
romboidales de oxalato de calcio. Los 
radios de 3 (2-4) series, moderadamente 


finos 33 (16-66) pm, extremadamente 
bajos 171 (96-360) pm y numerosos 
10 (7-13) / mm. Presentan estratificación; 
conteniendo almidón en albura y gomas 
en duramen. Las fibras son de tipo libri- 
forme, de longitud mediana 1178 (860- 
2000) pm, de diámetro fino 19(8-30) pm y 
paredes gruesas 7(2-14) pm. Presentan 
gomas en el lumen hacia el duramen. 

Densidad relativa: Albura 0,73 
(mediana) y duramen 0,75 (alta) (iawa, 
1989). 

2. Diferencias en las características ana- 
tómicas microscópicas (tablas 1 y 2) 

Se encontraron tres grupos de 
características significativamente dife- 
rentes. 
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a) Entre las zonas, entre árboles y en 
la interacción árboles-zona para la altura de 
las punteaduras (R 2 =0,40, p<0,01). 

b) Entre árboles y entre la interac- 
ción árboles*zona para: el número de 
vasos/mm 2 (R 2 0,12, p<0,01), el diámetro 
de vaso (R 2 = 0,18, p<0,01), la anchura 
de las punteaduras (R 2 = 0,18, p<0,01), el 
número de radios / mm (R 2 = 0,42, 
p<0,01), la altura de radios (R 2 = 0,21, 
p<0,01), la anchura de radios (R 2 = 0,48, 
p<0,01), la longitud de fibras (R 2 = 0,35, 
p<0,01) y el grosor de las fibras (R 2 = 
0,58, p<0,01), 

c) Entre árboles para: la longitud de 
vasos (R 2 = 0,26, p<0,01 ) y el diámetro de 
fibras (R 2 = 0,12, p<0,01). En general se 
observa que aunque existen diferencias 
significativas de las características anató- 
micas en los diferentes niveles, éstas son 
muy bajas. 


3. Análisis multivariado de las caracte- 
rísticas anatómicas 

Comparación entre árboles 

En la figura 3 se presenta la disper- 
sión de los 396 especímenes en un 
espacio definido por las dos primeras 
funciones del adm. Las dos primeras 
funciones discriminantes explican el 
72,6% de la variabilidad observada y 
ambas son significativas (p<0,0001). 
Aunque la k de Wilks para las dos 
primeras funciones es relativamente baja 
(0,303) y significativa (p<0,0001), se 
puede observar que existe una muy alta 
heterogeneidad dentro de los seis grupos 
(árboles) en las variables consideradas, 
de tal manera que los grupos no se 
separan claramente en el espacio discri- 
minante (Figura 3). 

En la tabla 3 se presenta la matriz de 
estructura con los valores de correlación 


4-r 
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Figura 3. Diagrama de dispersión de 396 muestras provenientes de seis árboles en un 
espacio discriminante definido por las dos primeras funciones del ADM que explican el 

73% de la variabilidad de la matriz. 
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entre las variables originales y las dos 
primeras funciones del adm. Las variables 
anatómicas que más contribuyen a la 
separación entre los grupos fueron, para 
la primera función en sentido positivo: el 
grosor de las fibras, la anchura de radios 
y la longitud de las fibras, y en sentido 
negativo la longitud de los vasos. Para la 
segunda función en sentido positivo: la 
anchura de los radios y la altura de las 
punteaduras y en sentido negativo el 
número de radios/mm. 


Comparación entre épocas y entre 
zonas 

El adm realizado para comparar 
registros de la estación seca con la esta- 
ción lluviosa muestra que, aunque la 
separación entre grupos es significativa 
(p<0,01), es muy poco clara, lo que se 
refleja en una k de Wilks de 0,774. Esto 
se observa claramente en la figura 4A, 
donde se aprecia la pobre, aunque signi- 
ficativa separación entre los grupos. 
Considerando las correlaciones entre las 


Tabla 3. Matriz de estructura del ADM realizado con 396 muestras y 11 variables para 
seis individuos (árboles), dos épocas (secas y lluvias) y dos zonas (albura y duramen). 
Los valores representan la correlación entre las variables y las funciones discriminantes 
extraídas. Se muestran con negritas las variables (r > 0.300) que más contribuyen a la 
separación entre los grupos para cada función discriminante. 




ÁRBOL 


ÉPOCA 

ZONA 

Características de elementos de vaso 

Variable 
Numero / mm 2 
Diámetro 
Longitud 

Altura punteadura 
Anchura punteadura 

ADM1 

0,091 

0,219 

- 0,349 

-0,113 

0,185 

ADM2 

-0,251 

-0,053 

0,399 

0,412 

0,084 

ADM1 

0,285 

0,260 

- 0,416 

-0,097 

0,060 

ADM1 

-0,021 

-0,073 

-0,289 

0,628 

0,132 

</j 

o 

Anchura 

0,557 

0,445 

0,005 

0,035 

'-o 

05 

o : 

Altura 

-0,101 

0,071 

0,098 

-0,051 


Número / mm 

-0,294 

-0,440 

-0,167 

0,251 

</5 

05 

Diámetro 

0,011 

0,242 

- 0,347 

-0,018 

s_ 

-Q 

Ll 

Grosor de pared 

0,569 

-0,199 

0,282 

0,314 


Longitud 

0,311 

0,089 

0,655 

-0,138 
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variables y la función discriminante (Tabla mayor longitud de los vasos (r = -0,416) y 

3), se obtiene que los tejidos de los diámetro de las fibras (r = -0,347), mien- 

árboles muestreados en la temporada de tras que en la temporada seca aumenta la 

lluvias estarían caracterizados por una longitud de las fibras (r = 0,655). 



Figura 4. Estimadores de Densidad por Kernel (EDK) a lo largo de la función 
discriminante extraída para (a) época (secas y lluvias) y (b) zona (albura y duramen). 
EDK generados a partir de la función ponderada de Epanechnikov con amplitud de 
ventana óptima de Haerdle de 0.837 (época) y 0.869 (zona). 
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En el caso de la comparación entre 
zonas de albura y duramen se presenta 
una situación similar, pues aunque la X, de 
Wilks es significativa (p<0,001), su valor 
es muy cercano a la unidad (0,838). Esto 
también revela una muy pobre separación 
entre grupos (Figura 4B). En este caso, la 
separación entre ambas zonas estaría 
determinada por la altura de la puntea- 
dura y el grosor de las paredes de las 
fibras, sin embargo, se sugiere un mayor 
tamaño de muestra para sustentar la 
tendencia de estas variables. 

índices de vulnerabilidad, de colapso 
de fibras y vasos 

Con los anova se encontró que para 
el índice de vulnerabilidad y de resis- 
tencia al colapso en vasos sólo hay dife- 
rencias significativas entre árboles (R 2 = 
0,72, p<0,01, R 2 = 0,18 p<0,01). Para el 
índice de colapso de fibras las diferencias 
estuvieron entre árboles, entre zonas y en 


la interacción árbol*zona (R 2 = 0,44 
p<0,01) (Figura 5 y tablas 1 y 2). 

Conductividad hidráulica teórica 

En la conductividad hidráulica de 
vaso (K h vaso) (R 2 = 0,16 p<0,01) y de 
área (K/ 7 ¿rea) (R 2 = 0.06 p<0,01) solo 
hubo diferencias entre árboles y en la 
interacción árbol*zona (tablas 1 y 2). 

Densidad relativa y contenido de 
humedad 

Los resultados de los análisis de 
varianza para la densidad relativa 
(R 2 = 0,58 p<0,01) y para el contenido de 
humedad: (R 2 = 0,85 p<0,01) sugieren 
que hay diferencias significativas entre 
árboles, entre zonas y entre su interac- 
ción (tablas 1 y 2). En la tabla 4 se 
presenta la comparación entre las medias 
de la densidad relativa por épocas (secas 
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Figura 5. Representación de los errores estándares de los índices de resistencia al 
colapso de vasos y fibras para los seis árboles. 
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Tabla 4. Comparación entre la época de secas y lluvias y entre albura y duramen en la 
densidad relativa y el contenido de humedad. Se presentan los valores promedio 

(error estándar) de las variables. 



SECAS 

LLUVIAS 

P 

ALBURA 

DURAMEN 

P 

Densidad relativa 
s/u 

0,759 

(0,0008) 

0,743 

(0,0437) 

0,7576 

0,728 

(0,0005) 

0,776 

(0,0006) 

0,0005 

Contenido de 
Humedad 

í%) 

27,174 

(7,410) 

35,714 

(3,085) 

0,00001 

36,174 

(4,069) 

25,856 

(6,568) 

0,0005 


y lluvias) y zonas (albura y duramen). Se 
puede apreciar que es significativamente 
mayor el contenido de humedad en las 
muestras provenientes de la temporada 
de lluvias que en las de secas, aunque no 
se observaron diferencias significativas 
en la densidad. Por otro lado, obser- 
vamos claras diferencias entre la albura y 
el duramen para ambas variables. 


DISCUSIÓN 

Características anatómicas 

Las características anatómicas de 
especies de la subfamilia Papilionoideae 
estudiadas en México presentan, en 
general, albura y duramen bien diferen- 
ciados por su color, con porosidad difusa, 
arreglos de poros solitarios a radiales; 
también se reportan tangenciales, agru- 
pados y diagonales. El tamaño de los 
vasos varía dependiendo del tipo de vege- 
tación del que procedan, pero de acuerdo 
con los datos extraídos de la literatura, se 
clasifican en moderadamente pequeños, 
medianos y algunos muy grandes, como 
es el caso de Erythrina (Barajas, 1989); la 
longitud de los vasos va de muy cortos a 
moderadamente cortos. El parénquima 
axial muestra una gran diversidad: en 
bandas confluentes, aliforme, aliforme 


confluente, difuso marginal y vasicéntrico. 
Los radios son en general homogéneos y 
la mayoría pueden tener de 2 a 4 series, 
son extremadamente finos a moderada- 
mente finos, extremadamente bajos y en 
algunos casos bajos y en general nume- 
rosos y en ocasiones estratificados. Las 
fibras son la mayoría libriformes y se 
detecta presencia de fibrotraqueidas. G. 
sepium se ajusta a esta estructura general, 
aunque resalta que sus vasos son extre- 
madamente cortos y en ocasiones puede 
presentar diámetros moderadamente 
grandes (>200 ¡im). 

En relación a la variación de la 
madera de los árboles estudiados de G. 
sepium, se observó que para el diámetro 
de los vasos no hay ninguna tendencia 
entre épocas, aunque el anova muestra 
una diferencia entre árboles y su interac- 
ción con las zonas, lo que podría estar 
explicado por el efecto de los individuos 1 
y 6 (tablas 1 y 2). Para la longitud de los 
vasos se observaron diferencias entre 
árboles, lo que se explica por el compor- 
tamiento del individuo 5 (tablas 1 y 2). 
Con el adm se observa la contribución de 
estos individuos para la separación entre 
árboles y estaciones (Tabla 3 y Figura 4). 
De lo anterior se observa que en las dos 
zonas se presentan diámetros y longi- 
tudes de vaso semejantes en secas y 
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lluvias, lo que sugiere que el volumen de 
agua que se transporta en los vasos 
funcionales es similar en ambas épocas, 
aunque el área conductora total dependa 
del número de vasos funcionales en 
secas o lluvias. 

El adm y el anova para la anchura de 
las punteaduras muestran que no hay 
diferencias claras entre árboles, zonas y 
épocas, sin embargo, en las tablas 1 y 2 
resalta que los árboles 1, 2, 4 y 5 
presentan diferencias. Para la altura de 
las punteaduras intervasculares se 
observa que las diferencias no son muy 
claras entre árboles y zonas, aunque el 
análisis numérico lo sugiere. 

Por otro lado, para estudiar el efecto 
que tienen las dos épocas estudiadas 
sobre este carácter, los autores sugieren 
un tamaño de muestra mayor. Las pun- 
teaduras intervasculares y de vaso- 
parénquima se presentan en toda la 
superficie interna del vaso (Figura 2B-H). 
Debido al tamaño y a la forma de las 
punteaduras se sugiere que la estructura 
cilindrica interna de cada una y su pared 
sólida pueden incrementar el valor del 
índice de prevención de colapso del vaso 
y las presiones internas del agua para 
reducir la formación de embolismos. 
Asociado a esto, las paredes de las tílides 
ejercerían una presión adicional sobre las 
paredes internas de los vasos, lo que 
incrementaría su resistencia al colapso 
(Figura 2E). Lo anterior sugiere que 
durante el desarrollo de las tílides que se 
forman cerca del cambium vascular, las 
paredes de estas células ejercen presión 
sobre la pared interna del vaso, que es 
adicional y en magnitud contraria a la de 
la transpiración en cavidades sin tílides. 

En cuanto al parénquima de radio, 
los análisis adm y anova muestran que la 
anchura del parénquima de radio es la 
que contribuye a que haya diferencias 
entre árboles y épocas, aunque estos 
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resultados pueden estar influidos por los 
valores del individuo 1 (tablas 1 y 2). 
Debido a esto no se observa una 
tendencia clara hacia el incremento de 
sus dimensiones o en número entre 
zonas o entre épocas, es decir, no se 
vislumbra un efecto de condiciones del 
crecimiento o de épocas (tablas 1 y 2). Es 
posible considerar que los radios son 
estructuras poco variables y estables y 
que alcanzan sus dimensiones finales 
relativamente rápido durante el 
desarrollo. 

En cuanto a las fibras, los análisis 
sugieren que los rasgos estudiados de las 
fibras son los que más contribuyen a 
explicar las diferencias entre zonas y 
épocas en la madera de Gliricidia sepium. 
De acuerdo con los datos, las fibras 
aumentan de longitud en la albura de la 
época de secas, pero en general son un 
poco más largas hacia el duramen. Las 
tablas 1 y 2 muestran que estas tenden- 
cias realmente dependen de la variabi- 
lidad de cada árbol y por tanto se sugiere 
que es un carácter que debe estudiarse 
con un tamaño de muestra mayor. De 
acuerdo con Chow (1971), las fibras 
pueden tener patrones de variación poco 
regulares en condiciones de tensión o 
compresión, lo que observó en la madera 
de especies de Liquidámbar, que 
presentó un incremento en la longitud de 
sus fibras hacia la corteza en los primeros 
10 cm y después decreció hacia la peri- 
feria. Se sugiere realizar estudios con un 
mayor tamaño de muestra para definir las 
tendencias en especies tropicales de 
selvas bajas caducifolias. 

En general se observa que solo 
algunos caracteres celulares logran 
separar la albura del duramen y resalta 
que las dimensiones celulares no son tan 
importantes para separar ambas zonas y 
épocas (tablas 1 y 2), lo que significa que 
tanto albura como duramen son seme- 
jantes en los periodos de secas y lluvias. 
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Cabe resaltar que las diferencias en los 
contenidos de humedad y la densidad 
relativa son variables más importantes 
para diferenciar zonas y estaciones 
(Tabla 5). 

Oluwafemi y Adegbenga (2007) 
encontraron que el contenido de humedad 
(C.H.) para G. sepium fue muy bajo (8,44- 
8,80%), con mayor humedad hacia la 
corteza, lo que contrasta con los valores 
obtenidos para G. sepium en este trabajo, 
que presentó una media de C.H. para la 
albura de 33% en secas y 39% en lluvias 
y para el duramen 20% en secas y 31% 
en lluvias. Asimismo, la diferencia de 
contenidos de humedad entre G. sepium 
de las plantaciones de Nigeria y la vera- 
cruzana puede estar explicada por el agua 
de los escurrimientos montañosos de la 
Sierra Madre Oriental. Elevitch y Francis 
(2006) y Quintanar etal. (1997) han seña- 
lado y discutido que la madera es pesada, 
lo que se debe a su estructura celular, 
como lo indican ampliamente los resul- 
tados de este trabajo: presencia de 
paredes gruesas en fibras y vasos, con 
presencia de tílides que se esclerosan 
hacia el duramen; extractivos y cristales 
presentes en cantidades abundantes 
dentro de los lúmenes de las fibras, de los 
vasos, de las tílides esclerosadas y en el 
parénquima axial y radial. A este respecto 
Oluwafemi y Adegbenga (2007) discuten 
que estos aspectos pueden estar expli- 
cando las propiedades mecánicas y 
físicas que encontraron para árboles de 
Gliricidia sepium en Nigeria. Estos autores 
mencionaron que no hay diferencias signi- 
ficativas en las contracciones radiales de 
albura y duramen. Los mismos autores 
destacan que estos individuos mostraron 
contracciones tangenciales más altas que 
las radiales (4,3%), las que disminuyeron 
hacia la corteza. Este valor fue conside- 
rado como moderado para una madera 
tan densa, lo que la hace muy estable 
para el servicio en exteriores. Asimismo, 
encontraron que los valores de mor y moe 


en condición verde eran más altos hacia la 
médula que hacia la periferia, 119,8 MNnr 2 
y 6629,2 MNnr 2 respectivamente; valores 
similares a los de especies comerciales 
como Diospyros y Lophira que se usan 
para construcción, aunque están por 
debajo de Teca (mor =136,1 MNm* 2 ). 
Debido a esto, los autores sugirieron que 
G. sepium puede sustituir la madera de 
especies de importancia económica actual 
y que ya son escasas. También agregaron 
que los altos valores de mor y moe hacia 
la médula podrían estar determinados por 
la densidad y que los valores de resis- 
tencia al impacto están muy correlacio- 
nados con la contracción radial, lo que 
asocian con la densidad. Estos autores 
destacan la relación mor/moe (0,83), para 
esta especie, lo cual es un valor de R 2 
muy alto. En este sentido, particularmente 
resalta lo que Swenson y Enquist (2007) 
sugieren acerca de la densidad como 
propiedad indicadora no solo de caracte- 
rísticas mecánicas, sino como una 
propiedad relacionada con la eficiencia 
hidráulica, lo que está intimamente rela- 
cionado con la flexión de vasos conduc- 
tores y de aquellos cavitados en árboles 
en pie. 

Como se ha resaltado, y de acuerdo 
con Oluwafemi y Adegbenga (2007), la 
densidad incrementa los valores de los 
mor y moe en la G. sepium de Nigeria, lo 
que sugiere que la especie veracruzana 
pueda tener comportamientos mecánicos 
similares, por lo que se recomienda 
comparar estas propiedades para México 
y Nigeria. Aunque estos autores no 
presentan valores de densidad relativa, la 
de los individuos veracruzanos es 
mediana para albura y alta para duramen 
(iawa, 1989). Debido a la baja variación 
de su estructura celular se sugiere que G. 
sepium es una madera muy estable y 
posiblemente los datos de Oluwafemi y 
Adegbenga (2007) dan luz sobre la rela- 
ción entre la estructura anatómica y las 
propiedades físicas y mecánicas. 



Madera y Bosques 15(2), 2009:71-91 

índices hidráulicos y de resistencia al 
colapso 

Dadas las características geomorfo- 
lógicas de las zonas costeras que 
incluyen la localidad de estudio, los escu- 
rrimientos provenientes de la Sierra 
Madre Occidental hacen un buen aporte 
de agua dulce que garantiza que la zona 
esté bien irrigada y bien drenada, debido 
a que los suelos de esta zona se caracte- 
rizan por ser poco desarrollados, 
presentar pocas arcillas, y por tanto 
muchos macroporos, lo que hace que 
sean de fácil infiltración. Por otro lado, de 
los datos de precipitación y temperatura, 
se observan diferencias marcadas entre 
los meses de enero y agosto, lo que 
podría sugerir un efecto en la estructura 
del xilema secundario. De las condiciones 
anteriores se percibe que el sustrato de la 
localidad recibe los escurrimientos de las 
montañas, y aunque los suelos son de 
tipo arenoso, el aporte por este proceso 
podría permitir que estas plantas no 
experimentaran estrés híd rico, lo que en 
principio podría explicar un efecto no 
apreciable de estos factores físicos en el 
xilema, dada la ausencia de porosidad 
anular y otros rasgos en los elementos de 
vaso y punteaduras, entre otras caracte- 
rísticas estructurales ya revisadas. 

Carlquist (1977) sostiene que entre 
más anchos y más cortos son los vasos, 
son más resistentes al colapso y entonces 
a la deformación, adicionando algunas 
otras características morfogenéticas entre 
las que resalta la naturaleza de las 
paredes de las punteaduras intervascu- 
lares, lo que hace más eficiente la conduc- 
ción y el ascenso del agua, por lo que a 
continuación se argumentan las razones 
por las cuales los autores consideran a G. 
sepium como una especie con xilema muy 
eficiente en la conducción hidráulica: El 
diámetro de sus vasos (179 ¿im en 
promedio) sugiere que el agua que 
asciende con presiones negativas, podría 
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ocasionar colapso, cavitación y embolismo 
en estas células, según lo muestran los 
valores de I. C. vasos y de vulnerabilidad 
(Tabla 1 y Figura 5). Los índices de resis- 
tencia al colapso expresan la capacidad de 
implosión de una estructura cilindrica 
cuando depende del grosor de sus 
paredes, lo que está asociado a la 
densidad (Hacke et al., 2001), por lo que 
resalta que el índice de resistencia 
ICyasos presenta un valor muy bajo con 
respecto al de las fibras (ICf¡b ras ). La defor- 
mación que sufre la pared de los vasos 
puede ser suficiente para colapsarlo 
cuando hay ascenso de agua; sin 
embargo, los valores altos de resistencia a 
la implosión y el grosor de las paredes de 
las fibras impiden el colapso de los vasos 
(Sperry, 2003; Larcher, 2003 y Jacobsen et 
al., 2005). Esto último también explica la 
relación entre el IC vasos , ICf¡bras y la alta 
densidad relativa de la madera de G. 
sepium, argumentado por Hacke (2001) y 
Swenson y Enquist. (2007). Por otro lado, 
el índice de vulnerabilidad está represen- 
tando pocos vasos en relación a su 
número/mm2, de manera que si un vaso se 
cavita, debido a su diámetro, entonces el 
área conductora del xilema se reduce. El 
desarrollo de tílides en vasos de albura 
cercanos al cambium vascular puede 
ocurrir después de que haya cavitación y 
embolismo. Este proceso se ve favorecido 
por la disposición y cantidad de parén- 
quima axial y radial (que reducen el 
número de vasos/mm 2 ), así como por el 
tamaño y abundancia de las punteaduras 
vaso-parénquima en la superficie del 
elemento de vaso. Como ya se mencionó, 
las tílides se expanden hacia el interior de 
la cavidad de los vasos hasta ocluirlos 
completamente, lo que de acuerdo con 
Braun (1984), puede incrementar el 
volumen de tejido accesorio. Las paredes 
de las tílides pronto se engruesan y, en el 
duramen, las cavidades de éstas se llenan 
de cristales y gomas. Esto último permitiría 
un incremento en el valor del IC en vasos 
disfuncionales de albura (Hacke et at., 
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2001) que ayudaría a resistir la presión de 
la deformación de las paredes de los 
vasos funcionales contiguos. Adicional- 
mente, el área de las paredes de las fibras 
contiguas a los vasos funcionales incor- 
pora resistencia al colapso, expresado 
macroscópicamente en la densidad 
(Sperry, 2003; Baas, 2004; Jacobsen etal., 
2005; Swenson y Enquist, 2007). 

La tabla 5 muestra valores inéditos 
de K/j.teórica de vaso y área respecto de 
algunas leguminosas de Oaxaca que se 
compararon con los de 6. sepium. Se 
observa que el valor de la relación Kfr_¿ rea 
/ kh- V a S0 de G. sepium es muy bajo con 
respecto al de las otras especies. G. 
sepium presenta en general un 
vaso/mm 2 , por lo que su Kft. te ór¡ca de vaso 
es muy similar a K^teórica P or área, situa- 
ción que contrasta con los valores de 
conductividad de las otras especies en 
las que la Kh-\ e ór\ca del vaso es muy 
similar a la de G. sepium pero la K^eórica 
del área se incrementa notablemente. 
Esto merece atención, puesto que las 
especies de Oaxaca proceden de selvas 
bajas caducifolias, donde se registra una 
precipitación pluvial mínima, en el mes 
más seco, de 1 ,3 mm (enero) y máxima 
de 93 mm (junio). Los valores de K^eórica 
por área están demostrando la impor- 
tancia del incremento en el número de 
vasos/mm 2 cuando existe estrés hídrico. 


De acuerdo con los resultados de los 
anova se observa que los valores de 
Kft-vaso y Kfi-á rea para G. sepium muestran 
conductividades similar en la albura de 
ambas épocas en todos los individuos, lo 
que sugiere que la fórmula xilemática de 
G. sepium no corresponde a condiciones 
de estrés hídrico o xéricas. Sin embargo, 
será importante que en futuras investiga- 
ciones se estudien los valores de K^. vaso , 
Kh-á rea. relacionándolos con el número de 
vasos funcionales por época. 

CONCLUSIONES 

La madera de G. sepium presenta 
vasos hidráulicamente eficientes para la 
conducción (K^teórica. diámetros, longi- 
tudes y valores bajos de resistencia al 
colapso en vasos), lo que sugiere un 
xilema muy conductor, muy vulnerable a 
la cavitación y al embolismo, conside- 
rando su hábitat; sin embargo, la albura y 
el duramen presentan valores medianos y 
altos de densidad relativa, que elevan la 
resistencia al colapso y explican su 
fórmula xilemática mesomórfica, con 
tílides que ocluyen las cavidades de los 
vasos cercanos al cambium vascular y 
abundancia del parénquima axial. Esta 
situación puede explicar el uso tan amplio 
que se le ha dado y su distribución 
geográfica actual, por lo que resulta una 
especie interesante para sustituir 


Tabla 5. Comparación de valores de K^. te órica (m*kg*MPa- 1 *s- 1 ) calculada para vasos 

y para el área de diferentes leguminosas. 


Valores de 
K h 

G. sepium 
Veracruz 

Piptadenia 

flava 

Oaxaca 

Prosopis 

laevigata 

Oaxaca 

Acacia 

farnesiana 

Oaxaca 

Cercidium 

plurifoliolatum 

Oaxaca 

de vaso 

2,98 x 10-u 

4,36 x 10-12 

1,16 x 10-n 

2,52 x 10-12 

5,05 x 10-12 

de área 

1,68 x 10-10 

8,42 x 10-10 

7,48 x 10-10 

2,03 x 10-10 

1,07 10-09 

K h área ! K h vaso 

5,4 

183 

67 

152 

199 
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maderas comerciales escasas. Esto 
último derivado de la similitud de las 
características anatómicas y la densidad 
de la albura y el duramen que la sugieren 
como una madera dimensionalmente 
estable. 
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del muérdago enano (Arceuthobium sp.) 
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RESUMEN 


El Nevado de Toluca fue decretado como Parque Nacional, reconociendo su importancia como 
sitio de recarga de los acuíferos, así como el valor de sus recursos naturales. Hoy es, sin embargo, 
escenario de una acelerada depredación tanto humana como natural. El muérdago es una de las 
plantas parásitas más importantes del grupo de las heterófitas debido a que causa serios perjuicios 
sobre otros vegetales y que se ha convertido un problema reciente debido a los niveles de infestación 
del parásito. En este trabajo se obtuvieron datos que permitieron conocer la distribución espacial de 
este parásito y mapas sobre su densidad en la zona de estudio; para esto se utilizaron plataformas 
de software ya existentes que han sido probadas ampliamente. El resultado de este estudio permitirá 
realizar programas de manejo de este parásito que sea sustentada por información confiable, que 
presumiblemente lleve a lograr un control adecuado de este problema. 

PALABRAS CLAVE: 

Coniferas, heterófitas, muérdago enano, Nevado de Toluca y SADIE. 


ABSTRACT 

The Nevado de Toluca was declared a National Park, recognizing its importance as a site for 
aquifer recharge, as well as the valué of its natural resources. Today, however, it is the scene of a rapid 
depredation both human caused and natural. The mistletoe is one of the most important plant para- 
sites group of heterophytes because they cause serious harm on other plants, which has become a 
problem recently because of the infestation levels of the parasite. In this study data were obtamed that 
permit to know the spatial distribution of this parasite and maps its density iri the study area; they can 
be used for software platforms already in place that have been tested extensively. The results of this 
study will permit management programs of this parasite that are broked up by reliable information, 
presumably leading to pretend adequate control of this problem. 
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Coniferous, heterophyte, mistletoe, Nevado de Toluca and SADIE. 
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INTRODUCCIÓN 

El Parque Nacional Nevado de 
Toluca, es una región montañosa que se 
ubica en la parte central del Estado de 
México, al suroeste del Valle de Toluca. 
Según protinbos, tiene un total de 
48,571 has de las cuales 77.42% tiene 
uso forestal y de esto más de 50% corres- 
ponde a bosques con poca densidad, 
combinados con usos agrícolas y gana- 
deros. El Nevado de Toluca o Xinantécatl, 
según cepanaf (1972), es un volcán 
extinto que alcanza una altura de 4558 
msnm. El Parque tiene una extensión de 
54,000 ha, pobladas por abundante vege- 
tación de pinos, cedros y abetos en 
donde habitan diferentes especies 
animales. 

Los pinos son un género ( Pinus ) de 
árboles o, raramente, arbustos, de la 
familia Pinaceae, subfamilia Pinoideae 
que presentan una ramificación frecuen- 
temente verticilada y más o menos 
regular. En México se tiene un gran 
número de especies diferentes de pinos 
(género Pinus), mayor que en cualquier 
otra área en el mundo, de tamaño seme- 
jante (Sánchez-González, 2008). Por 
consecuencia, es un buen lugar para 
llevar a cabo la observación e interpreta- 
ción de la variación y evolución de estas 
especies. 

En el Parque Nacional Nevado de 
Toluca (pnnt) se encuentran dos espe- 
cies de pino ( Pinus montezumae, Pinus 
hartwegii), las cuales son atacadas por el 
muérdago enano. El muérdago es una de 
las plantas parasitas más importantes del 
grupo de las heterófitas, debido a que 
causa serios perjuicios sobre otros vege- 
tales a pesar de que se sabe que existe 
una relación parcial de alimentación entre 
huésped y hospedero, la que nunca 
puede ser benéfica y sí a menudo muy 
dañina para este último (González, 1989). 
Los muérdagos contienen generalmente 


clorofila, la cual les permite elaborar 
algunos de sus propios alimentos a partir 
de materiales alimenticios inorgánicos 
tomados de su hospedero (conifera); 
viven sobre las ramas de árboles y 
arbustos, siendo tan extraordinaria la 
aptitud de esta fanerógama para el para- 
sitismo, que se les encuentra sobre hojas 
y coniferas de todas las alturas y latitudes 
y en algunos casos produciendo un auto- 
parasitismo en especies de su propia 
familia (Verduzco, 1976). 

En México existen tres especies 
principales de muérdago enano, Arceu- 
thobium abietis religiosae, el cual es pará- 
sito especifico de los árboles de oyamel, 
A. globosum subsp. grandicaule, que es 
la especie de mayor distribución en la 
parte central del país y A. vaginatum 
subsp. vaginatum, que es la especie de 
muérdago más ampliamente distribuida 
en México (Rodríguez, 1985) Aparente- 
mente no es más que una modesta 
especie de parásito de las coniferas, pero 
el muérdago enano posee un sistema 
único de reproducción que tiene maravi- 
llados a los científicos. A diferencia de los 
ejemplares de otras especies hermanas, 
cuyas semillas llegan a los árboles 
después de que un pájaro se las coma y 
las deposite en el lugar apropiado, este 
tipo de muérdago lanza su semilla a una 
distancia que puede llegar a los 20 
metros, gracias a una suerte de bomba 
de agua que se pone en marcha tras un 
periodo de gestación de 18 meses. 
Durante todo ese tiempo, la planta 
acumula el líquido elemento, hasta 
conseguir la presión necesaria para 
dispersar sus semillas. Los tallos o brotes 
aéreos son la parte del parásito que sale 
fuera del tejido de la hospedera, donde se 
desarrollan los órganos reproductores, se 
originan en yemas del sistema endofítico 
y una vez que éste se ha establecido, 
están cubiertos por una capa de tejido 
epidérmico; la formación epidémica 
original es sustituida por un tejido (epidér- 
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mico auricular secundario, pero los tallos 
no llegan a tener nunca una corteza 
normal). El muérdago enano puede ser 
parasitito y vivir indefinidamente sin la 
emisión de tallos. La presencia de 
almidón en el sistema endofítico de 
Arceuthobium sugiere la producción de 
carbohidratos por estas plantas (Gutié- 
rrez, 1970) 

Antes se creía que los muérdagos y 
hospederas vivían en simbiosis, hoy en 
día, definitivamente se reconocen las 
relaciones entre parásito y hospedera; 
hay sin embargo algo de duda respecto a 
la transmisión del muérdago hacia su 
hospedera (Gilí y Hawksworth, 1961) 
citado por Valdivia (1963). La velocidad 
de intercambio entre hospedera y pará- 
sito mediante sustancias coloreadas ha 
revelado que al paso de éstas es mucho 
más fácil de la hospedera al parásito que 
en dirección opuesta, pero no hay infor- 
mación cuantitativa sobre la velocidad del 
movimiento en cualquiera de las direc- 
ciones (Childs y Edgron, 1976). 

Las especies del género Arceutho- 
bium requieren para su desarrollo de 
abundante luz solar, estos requerimientos 
varían de acuerdo con el desarrollo 
mismo de la planta. La abundancia de luz 
solar, condición que se presenta en 
montes aclarados, es sumamente favo- 
rable para el desarrollo y diseminación 
del parásito, los brotes relacionados con 
vigor y abundancia en su presencia, 
mayormente cuando se encuentran infes- 
tando árboles residuales en áreas que 
han sido explotadas y donde el estrato 
inferior de la vegetación está formado por 
renuevo sobre el cual cae la semilla del 
parásito (Korstian, 1922) citado por 
Valdivia (1963) y Gutiérrez (1970). 

La edad del árbol tiene una impor- 
tancia muy grande respecto a la suscepti- 
bilidad a ser infestado, árboles viejos son 
menos susceptibles por la resistencia que 
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presenta su corteza engrosada a la pene- 
tración de la radícula de la semilla en 
germinación y por la carencia de tejido 
clorofiláceo expuesto. La susceptibilidad 
de una planta a ser parasitada disminuye 
proporcionalmente con su crecimiento, 
alcanzando un máximo aquéllas que 
tienen un diámetro de 10 a 15 cm (Gilí y 
Hawksworth, 1961). La variación y altitud 
dentro de la cual crece este parásito en 
México, según Gilí y Hawksworth, (1961) 
es de los 2,850 hasta 3,150 msnm. No 
hay una selectibilidad marcada con 
respecto a la exposición, encontrándose 
el muérdago en sitios de cualquier expo- 
sición. Las infestaciones del muérdago se 
desarrollan en pendientes moderadas 
como pronunciadas y los cambio bruscos 
de temperatura no afectan al muérdago 
enano (Gilí y Hawksworth, 1961) 

Varias pináceas son fácilmente 
infestadas por muérdago enano, ocasio- 
nándoles debilitamientos muy considera- 
bles que repercuten en su crecimiento e 
incremento. El primer signo visible de la 
infestación del muérdago es un ligero 
abultamiento causado por el estímulo 
local de la corteza, floema y xilema de la 
hospedera (Vega, 1976). A medida que 
transcurre el tiempo y el sistema endofí- 
tico crece dentro de la hospedera, la rama 
puede distorsionarse formando a veces 
“escobas de brujas”. Por regla general 
estas ramas infestadas sobreviven varios 
años a las que están libres del parásito, 
absorbiendo gran cantidad de nutrientes 
que normalmente irían hacia la parte 
superior de la copa de los árboles; estas 
ramas atrofiadas con hábitos de creci- 
miento anormal son de follaje muy denso 
y alcanzan un diámetro superior a los 
demás (Gutiérrez, 1970). El muérdago 
enano reduce el vigor del árbol, lo que 
trae como consecuencia una reducción 
en el crecimiento longitudinal que puede 
llegar a un 30-40% y el crecimiento en 
diámetro puede reducirse en un 40% en 
árboles altamente infestados (Hawks- 
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worth, 1961). Otra característica del 
ataque del muérdago es que predispone 
a los árboles infestados a un ataque 
intenso de descortezadores y barrena- 
dores, pudrición de la madera por 
presencia de algunos hongos, enferme- 
dades de la raíz, así como debilitamiento 
en el tronco y ramas que llegan a sufrir 
fácil rompimiento y caídas por el viento, y 
como consecuencia al caer árboles y 
ramas muertas son material propicio para 
incendios (Vega, 1976). 

Los principales problemas de las 
coniferas en ladera norte del Parque 
Nacional Nevado de Toluca (cauce del 
Lerma) son la deforestación, los incen- 
dios, las plagas y enfermedades. Dentro 
de estas últimas las crecientes pobla- 
ciones de muérdago enano se han 
convertido en un serio problema sanitario 
para los bosques de coniferas de la 
región. 

Los daños que causa el muérdago 
enano a los árboles son los siguientes: a) 
hipertrofia del fuste y de las ramas en el 
punto infestado; b) atrofia y posterior 
rotura del fuste y de las ramas de las 
partes infestadas hacia la cúspide debido 
al peso del parásito, de manera que los 
puntos del sujeto parasitado fácilmente 
terminan en penachos de muérdagos; c) 
produce la deformación de las ramas de 
los árboles afectados al culminar el 
ataque; d) mata de modo rápido principal- 
mente el arbolado joven, retarda grande- 
mente el crecimiento tanto en altura como 
en diámetro de los árboles de mayor 
edad; e) interfiere seriamente en la vida 
normal del árbol, provocándole a la postre 
la muerte al adquirir frondosidad el pará- 
sito y por consiguiente aumentar sus 
necesidades de sostén y nutrientes; f) al 
inicar la primavera, época de mayor creci- 
miento de los árboles, se retardan las 
funciones nulificando el desarrollo, altura 


y diámetro del sujeto, dándole a la postre 
una apariencia de achaparramiento; g) al 
debilitar grandemente a los árboles, éstos 
quedan predispuestos al ataque de 
insectos y microorganismos; h) en aque- 
llos lugares donde se han aplicado cortas 
como medida de control sobre los árboles 
intensamente atacados por el muérdago 
enano, pero que se han dejado algunos 
sujetos con parasitismo incipiente, viene 
una propagación extraordinaria del muér- 
dago enano como consecuencia del 
aclareo que sufrió la masa y la mayor 
cantidad de luz que reciben estas plantas; 
i) reduce considerablemente la produc- 
ción de semillas fértiles de los sujetos 
atacados y en muchas ocasiones llega a 
nulificarla totalmente. 

Como se aprecia, los daños 
causados por las poblaciones cada vez 
más alarmantes de muérdago enano en 
la región, han provocado un impacto 
negativo creciente sobre las poblaciones 
de coniferas en la zona del cauce del 
Balsas. Dicho impacto negativo ha sido 
observado pero no cuantificado hasta la 
fecha, hecho que ahora resulta de vital 
importancia para poder tomar medidas de 
control. 


OBJETIVOS 

Establecer el tipo de distribución 
espacial de las poblaciones de muérdago 
enano, parásitos de coniferas, en la 
ladera norte (cuenca del Lerma) del 
Parque Nacional Nevado de Toluca. 

Modelizar la distribución espacial de 
las poblaciones de muérdago enano 
mediante la utilización de las estadísticas 
no espacial (índices de dispersión y distri- 
buciones estadísticas) y espacial, con el 
Análisis Espacial por índices de Distancia 
(sadie). 
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Establecer la existencia de una esta- 
bilidad espacial y temporal de las pobla- 
ciones del muérdago enano. 

MATERIAL Y MÉTODOS 

El presente trabajo se realizó en el 
Parque Nacional Nevado de Toluca que 
está ubicado en el Estado de México, en 
los municipios de Zinacantepec, Toluca, 
Calimaya, Coatepec de Harinas, Tenango 
del Valle, Villa Guerrero, Almoloya de 
Juárez, Temascaltepec y Amanalco de 
Becerra; y situado en las coordenadas 
utm con 2 102 398 y los 2 134 700 de 
metros Norte y 399 833 y los 432 253 de 
metros Este. El Parque Nacional Nevado 
de Toluca es uno de los principales 
parques, ya que ocupa el primer lugar en 
su extensión, siendo su superficie de 
53 913 hectáreas, ocupando 52% de los 
diez parques del Estado de México. 

Para el estudio del comportamiento 
espacial del muérdago enano se estable- 
cieron en la ladera norte (cuenca del 
Lerma) del pnnt 10 parcelas experimen- 
tales de coniferas, de una hectárea cada 
una. Cabe señalar que cada una de las 
parcelas estaba situada a diferentes alti- 
tudes sobre el nivel del mar. Este estudio 
se realizó durante el año 2006. Se esta- 
bleció una malla rectangular de una 
hectárea de diámetro y se muestreó a 
cada 10 metros, en ambas direcciones, 
alrededor de 100 árboles por malla, cada 
árbol se georeferenció utilizando un dgps 
para obtener sus coordenadas espa- 
ciales, además se realizó un conteo de 
números de muérdagos en cada árbol 
muestreado y su identificación correspon- 
diente. Se realizó una exploración esta- 
dística de los datos originales de las 
poblaciones de muérdago enano para 
cada muestreó. Con base en lo anterior 
fue necesario la realización de una trans- 
formación logarítmica de los datos 
[logio(n+1)j para normalizarlos. 
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Estadística clásica 

Existen dos métodos generales para 
establecer la distribución espacial de los 
organismos (Binomial Negativa y 
Poisson) y los índices de dispersión (de 
dispersión y de Green); ambos se utili- 
zaron en este estudio para realizar la 
comparación pertinente entre lo hallado 
con la estadística no espacial y la estadís- 
tica espacial. Para todas las distribu- 
ciones estadísticas se utilizó el programa 
mlp de máxima verosimilitud (Ross, 1987) 
para ajustar los modelos a los datos obte- 
nidos. La bondad del ajuste fue exami- 
nada con un test c 2 (Sokal y Rohlf, 1995). 

Estadística espacial 

Análisis Espacial por índices de Distancia 

(SADlEj 

El objetivo de sadie es establecer el 
modelo espacial de una población mues- 
treada midiendo la distancia a la cual 
pueden desplazarse los individuos de la 
muestra observada. Perry y Hewitt (1991) 
desarrollaron un índice que utilizaba los 
datos colectados dentro de una cuadrí- 
cula ya preestablecida antes del mues- 
treo, determinando que era un índice 
biológicamente más descriptivo e infor- 
mativo de la distribución espacial de las 
poblaciones de insectos que el índice de 
dispersión, el cual depende directamente 
de la abundancia de la población. Perry 
(1995a) indicó que para datos recolec- 
tados en ubicaciones específicas, el uso 
de la distancia para la regularidad es muy 
adecuado. Demostró cómo distinguir no 
aleatoriedad en la forma de heteroge- 
neidad estadística, de la no aleatoriedad 
espacial, Perry (1995b) desarrolló y 
extendió el uso del índice de la distancia 
para la regularidad (l a ) para el estableci- 
miento de la estructura espacial de las 
poblaciones de insectos. Además, intro- 
dujo dos diagramas de diagnóstico como 
ayuda a la interpretación y un índice 
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nuevo para estimar el número de focos de 
agrupamiento de una población, el índice 
J a . Además, señaló que la escala de 
muestreo es un determinante importante 
del patrón espacial total. Alston (1996) 
corroboró que la distancia para la regula- 
ridad (D) proporciona bases más 
adecuadas para elaborar un índice, ya 
que el índice de agrupamiento desarro- 
llado por Perry y Hewitt (1991) tiene 
serias limitaciones para detectar múltiples 
grupos, amontonamientos o parches. En 
el presente trabajo se utilizaron el índice 
basado en la distancia para la regularidad 
l a y el índice J a , basado en la distancia del 
agrupamiento Perry (1995a,b) para esta- 
blecer el modelo de distribución de las 
poblaciones de muérdago enano. 

Estimación de los índices l a y J a 

Los datos recolectados en una 
cuadrícula prediseñada (conformada por 
unidades de muestreo), se consideran 
como un sistema de conteo de individuos, 
donde i = 1,„.,n unidades de muestreo. 
Se conoce además, la posición bidimen- 
sional (x¡, y¡) de cada unidad de muestreo 
y su conteo asociado, N¡. La distancia 
para la regularidad, D, es el valor mínimo 
de la distancia total que los individuos en 
la muestra pueden haberse movido, de 
una unidad de muestreo a otra, de modo 
que todas las unidades de muestreo 
contuvieran un número idéntico de indivi- 
duos. La solución se refiere a la manera 
óptima en la cual los individuos se move- 
rían desde cada unidad de muestreo con 
un conteo inicial más grande que la 
media, a otras unidades de muestreo con 
un conteo inicial más pequeño que la 
media. Si entonces los conteos obser- 
vados se permutan aleatoriamente entre 
las unidades de muestreo, de modo que 
la muestra que resulta sea un cambio o 
reajuste simple de la original, entonces P a 
(probabilidad de agregación) representa 
la proporción de muestras seleccionadas 
al azar con distancia para la regularidad 
tan grande como, o más grande que, el 


valor observado, D. Intuitivamente, con 
un valor grande de D podría esperarse 
que estuviera implicado un agregamiento 
o un agrupamiento, por ejemplo, un 
patrón espacial heterogéneo, e inversa- 
mente, un valor pequeño de D implicaría 
una regularidad, por ejemplo, un patrón 
espacial uniforme. Un valor de P a deri- 
vado de un número suficientemente 
grande de aleatorizaciones proporciona 
una prueba formal de aleatoriedad; la 
hipótesis nula de aleatoriedad espacial se 
puede rechazar, si P a < 0,025 (en favor 
de una hipótesis alternativa de agrega- 
ción), o si P a > 0,975 (en favor de la alter- 
nativa de regularidad) dando el usual 5% 
de probabilidad de rechazar la hipótesis 
nula cuando es verdad. Si la distancia 
media aritmética para la regularidad de 
las muestras aleatorias se denota como 
E a , entonces el índice de agregación, 
denominado / a , está definido como 
l a = D/E a . Usualmente se dice que una 
muestra es agregada si l a > 1 , la muestra 
es espacialmente aleatoria si l a = 1, y la 
muestra es regular si l a < 1. Un número 
total de 2000 aleatorizaciones son sufi- 
cientes para derivar los valores de los 
índices correspondientes. 

El término C denota la distancia para 
el agrupamiento, que es el valor mínimo 
de la distancia total que los individuos de 
la muestra deben moverse para congre- 
garse en una unidad. Este valor se 
encuentra más pronto que D, usando una 
búsqueda directa simple sobre todas las 
unidades de muestreo; a la unidad de 
muestreo con el valor mínimo se le deno- 
mina como el "foco” del agrupamiento. 
Permutaciones aleatorias de los conteos 
observados conducen a una proporción 
denominada Q a (probabilidad de agrupa- 
miento), con una distancia para el agru- 
pamiento tan pequeña, o más pequeña 
que el valor observado, C. Intuitivamente, 
para datos que comprenden un agrupa- 
miento individual, un valor pequeño de C 
denotaría un patrón espacial agregado; 
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inversamente un valor alto de C implicaría 
un patrón espacial regular. Análoga- 
mente, la hipótesis nula de aleatoriedad 
puede ser rechazada si Q a < 0,025 (en 
favor de la alternativa de agregación) o si 
Q a > 0,975 (en favor de la alternativa de 
regularidad), y si la distancia media para 
el agrupamiento para muestras aleatorias 
se denota como F a , entonces el índice de 
agregación J a , se define como J a = F a /C. 
Como en el caso del índice / a , valores de 
J a > 1 usualmente indican una muestra 
agregada, J a = 1 representan datos espa- 
cialmente aleatorios y J a < 1 muestras 
regulares. De esta manera, los valores 
del índice J a sirven para corroborar los 
resultados obtenidos con el índice l a . 
Además, este índice se utiliza para discri- 
minar entre patrones espaciales donde 
hay un único agrupamiento importante 
para el cual sus valores son significativa- 
mente mayores que la unidad, y en donde 
hay dos o más agrupamientos para los 
cuales su valor no es significativamente 
diferente de la unidad o incluso menor 
que ella. Para determinar la significación 
con respecto a la unidad se utiliza su 
probabilidad respectiva (Q a ) (Perry, 
1 998). Los valores de l a y J a para conteos 
aleatorios no están correlacionados, por 
lo que se pueden utilizar también 2000 
aleatorizaciones en el software empleado 
para obtener sus respectivos valores. El 
programa utilizado en este trabajo para 
determinar los valores y las probabili- 
dades de ambos índices fue el sadie 1.22 
(programa cedido por el doctor Perry). Se 
utilizaron mapas de densidad realizados 
con antelación mediante un análisis 
geoestadístico (Samper y Carrera, 1990) 
en el cual se empleó el Krigeado Ordi- 
nario*, para corroborar los resultados 
obtenidos con los índices del sadie. 


Método de interpolación desarrollado por el 
geólogo sudafricano D. G. Krige. 
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Cálculo de la superficie infestada con 
base en mapas de densidad 

Con base en mapas realizados con 
antelación utilizando el método del 
Krigeado Ordinario, se determinó el 
porcentaje de infestación establecido por 
el parásito en cada una de las parcelas 
experimentales. Para realizar lo anterior, 
se utilizó el programa Surfer 9.0. 

Estabilidad espacial y temporal 

Se compararon los mapas elabo- 
rados con anterioridad con el método del 
Krigeado Ordinario, con la finalidad de 
establecer si existía o no una estabilidad 
espacio-temporal de las poblaciones del 
muérdago enano. La comparación para 
establecer la estabilidad espacio- 
temporal de los mapas citados se llevó a 
cabo mediante el índice de Asociación del 
sadie llamado l m (Perry y Klukowsky, 
1997 y Korie etal., 2000). Si l m > 0, indica 
que existe una asociación o una estabi- 
lidad espacio-temporal entre los mapas sí 
existe un nivel de significación de 
Pm < 0.025. La determinación de la esta- 
bilidad se realizó con el sadie 1.22. 

RESULTADOS 

Con los datos obtenidos fue posible 
generar la modelización y mapeo del 
comportamiento espacial del muérdago 
enano. Se estableció el porcentaje de 
infestación en cada una de las parcelas 
de estudio y se logró determinar si existía 
estabilidad espacio-temporal de este 
parásito en las parcelas estudiadas, y con 
ello se consiguió determinar realmente el 
patrón de comportamiento espacial 
presente por parte del parásito. Los resul- 
tados obtenidos en cada una de las 10 
parcelas de estudio de una hectárea en la 
ladera norte (cuenca del Lerma) del 
Parque Nacional Nevado de Toluca, se 
detallan enseguida: 
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El objetivo general planteado en 
este estudio era modelizar y generar 
mapas de la distribución espacial del 
muérdago enano, parásito de coniferas, 
del Parque Nacional Nevado de Toluca. 
Para tal hecho se utilizaron dos métodos 
que permitieron establecer dicho compor- 
tamiento espacial. Además, se realizó 
una comparación entre los resultados de 
las parcelas, de tal manera que dicha 
actividad permitiera conocer con mayor 
precisión la distribución espacial de este 
parásito y al mismo tiempo se estable- 
cieran las ventajas y desventajas de cada 
método empleado. 

Se realizaron apartados específicos 
de los resultados alcanzados con cada 
uno de los métodos utilizados y de esta 
manera dejar más claro qué se logró 
alcanzar con cada uno de ellos. 

Distribución espacial del muérdago 
enano 

Estadística No Espacial 

En el Tabla 1 se observa el resumen 
estadístico del número de muérdagos 


enanos presentes en los árboles de coni- 
feras muestreados. Se logró determinar 
gracias al Coeficiente de Curtosis que no 
existía normalidad en los datos recolec- 
tados, por lo tanto fue necesario realizar 
una transformación de los mismos para 
poder utilizar los datos en la estadística 
no espacial y en el análisis. Para llevar a 
cabo esto se utilizó una transformación 
logarítmica de los mismos [(logio( n + 1)]. 
una vez realizado este proceso se cons- 
tató que los datos ya se distribuían de 
formal normal aplicando nuevamente el 
Coeficiente de Curtosis. Dicho coeficiente 
indica normalidad cuando registra valores 
entre -3 y 3. 

El intervalo de densidad de muér- 
dagos enanos por conifera varió entre 
4,282 y 1,472. Las parcelas con mayor 
densidad de muérdagos (3 y 10) se 
encuentran ubicadas a altitudes más 
bajas que las otras tres, habría que deter- 
minar si la altitud en realidad juega algún 
papel importante para determinar la 
densidad de este parásito, esto no se 
realizó en este estudio, pero sería muy 
interesante corroborarlo en estudios 
posteriores. En este ambiente hetero- 


Tabla 1. Resumen estadístico del parásito muérdago enano. 


PARCELAS 

X 

S 

S2 

COEFICIENTE COEFICIENTE 

DE CURTOSIS (1 a )DE CURTOSIS (2 a ) 

P1 

3,727 

3,981 

15,851 

4,471 

1,085 

P2 

3,381 

3,440 

11,834 

6,983 

-1,902 

P3 

4,282 

2,934 

8,608 

10,627 

0,441 

P4 

3,300 

3,333 

11,111 

19,552 

0,833 

P5 

2,300 

2,395 

5,735 

16,561 

-1,953 

P6 

1,472 

2,285 

5,224 

6,381 

1,267 

P7 

1,818 

3,059 

9,361 

10,772 

0,821 

P8 

2,300 

3,520 

12,395 

5,222 

-0,661 

P9 

3,281 

6,175 

6,175 

16,192 

1,993 

PIO 

4,036 

7,332 

53,760 

19,356 

0,240 
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géneo donde se ubicaron las distintas 
parcelas y en donde se realizaron por 
tanto los diferentes muéstreos, era de 
esperarse que existiera una considerable 
variabilidad de la repartición espacial de 
este parásito. La varianza de los datos se 
halló entre 53,760 y 5,224. 

En cuanto al tipo de distribución 
espacial de los muérdagos enanos 
encontrada en cada una de las parcelas 
experimentales, los resultados se hallan 
registrados en la tabla 2. El índice de 
dispersión señala que en todas las 
parcelas se encontró una distribución 
agregada del muérdago enano, con 
excepción de las parcelas 4, 6 y 9, donde 
indica una distribución aleatoria. Por su 
parte, el índice de Green confirmó dicha 
agregación, pero en este caso para todas 
las parcelas, indicando, en función de su 
bajo valor (Tabla 2), que la agregación era 
leve. Estos datos un tanto contradictorios 
señalan algunas de las limitaciones de 
estos métodos para determinar con preci- 
sión la distribución espacial del parásito. 
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Por otro lado, con base en las distri- 
buciones estadísticas, los resultados 
indican que en las parcelas 1 , 2, 4, 5, 6, 7, 
9 y 10 se ajustó una distribución binomial 
negativa (agregación) a los datos. En la 
parcelas 2, 5 y 6 también se efectuó un 
ajuste a la distribución de Poisson (alea- 
toria), por lo que en este caso no se 
puede discernir la real distribución espa- 
cial del muérdago enano en estas 
parcelas, ya que no es posible que el 
parásito tenga al mismo tiempo dos tipos 
de distribución espacial. Quizás la expli- 
cación a dicha dualidad puede estar en 
que para valores altos del parámetro k, 
como es relativamente el caso, la Bino- 
mial Negativa tiende a la distribución de 
Poisson, por lo que al realizar el ajuste 
estadístico por máxima verosimilitud se 
ajusten a los datos ambos modelos indis- 
tintamente. Los resultados obtenidos por 
el índice de dispersión y el índice de 
Green tampoco dejan clara la situación, 
ya que en los dos primeros casos indican 
una agregación y el caso de la parcela 6 
el índice de dispersión indica aleatoriedad 


Tabla 2. índices de dispersión y distribuciones estadísticas del parásito muérdago 

enano. 


PARCELAS ÍNDICE DE 

DISPERSIÓN 

ÍNDICE DE 
GREEN 

POISSON 

BINOMIAL 

NEGATIVA 

K 

P1 

2,32 s 

0,002 

NS 

S 

1,62 

P2 

1 ,56 s 

0,003 

S 

S 

6,92 

P3 

3,22 s 

0,004 

NA 

NA 

- 

P4 

0,77 ns 

0,001 

NS 

S 

4,73 

P5 

4,89 s 

0,005 

S 

S 

9,73 

P6 

0 , 86 ns 

0,001 

S 

S 

11,35 

P7 

2,71 s 

0,003 

NS 

S 

5,22 

P8 

3,39 s 

0,006 

NA 

NA 

- 

P9 

0,65 ns 

0,004 

NS 

S 

2,80 

PIO 

3,02 s 

0,007 

NS 

S 

5,65 


S=significativa; NS= no significativa; nivel de significación al 5%; NA= No Ajustada 






102 


Estudio de la distribución espacial del muérdago enano (Arceuthobium sp.) 


y el de Green, agregación. Esto 
demuestra nuevamente las contradic- 
ciones en las que pueden caer estos 
métodos y, por lo tanto, señalan sus limi- 
taciones. 

Los datos obtenidos en las parcelas 
3 y 8 no se pudieron ajustar a ninguna 
distribución estadística, ya que no fue 
posible obtener la convergencia de los 
algoritmos de ajuste por máxima verosi- 
militud. Esto es una muestra más de las 
limitaciones de los métodos de la estadís- 
tica no espacial. 

Estadística Espacial 

Análisis Espacial por índices de Distancia 
(sadieJ 

Los datos obtenidos con la aplica- 
ción de los índices del sadie se concen- 
tran en la tabla 1, en ella se observa que 
en el caso del índice l a el valor más alto 
se registró en la parcela 2, que fue de 


1,76 y el más bajo en la parcela 10, igual 
a 1 ,27. En todos los casos el índice l a fue 
significativamente superior a 1 (Tabla 3) 
lo cual señala una distribución agregada. 

En lo referente al índice J a el valor 
mayor se dio en la parcela 2, que fue de 
1 ,25 y el más bajo en la parcela 7, con un 
valor de 1,07. El índice J a en todas las 
parcelas experimentales también fue 
superior a la unidad, este dato confirma la 
agregación detallada en el índice anterior. 

Con los datos obtenidos con el 
índice J a , fue posible establecer la 
cantidad de centros de agregación que se 
hallaron en cada una de las parcelas 
estudiadas. En cuanto a este dato, el 
valor del índice J a en todas las parcelas 
no fue significativamente superior a 1 , lo 
que permite determinar que la distribución 
espacial de las poblaciones del parásito 
se concentraba en diferentes centros de 
agregación. 


Tabla 3. Valor de los índices l a y J a y sus respectivas probabilidades P a y Q a 

en la población de muérdago enano. 


Parcelas 

la 

Pa 

^a 

Qa 

P1 

1,66 

0,007 s 

1,21 

0,241 ns 

P2 

1,76 

0,015 s 

1,25 

0,274 ns 

P3 

1,32 

0,003 s 

1,15 

0,162 ns 

P4 

1,47 

0,010 s 

1,17 

0,205 ns 

P5 

1,70 

0,002 s 

1,23 

0,196 ns 

P6 

1,51 

0,008 s 

1,11 

0,254 ns 

P7 

1,47 

0,01 2 S 

1,07 

0,286 ns 

P8 

1,30 

0,004 s 

1,09 

0,128 ns 

P9 

1,55 

0,009 s 

1,16 

0,216 ns 

PIO 

1,27 

0,017 s 

1,19 

0,249 ns 


ns: no significativo al 5% s: significativo al 5% 
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Cálculo de la superficie infestada con 
base en mapas de densidad 

Lo detallado con los índices del 
sadie se observa con precisión en los 
mapas de densidad elaborados con este 
programa. 

Los mapas elaborados con el 
método del Krigeado Ordinario se deta- 
llan en la figura 1. La distribución agre- 
gada de las poblaciones del parásito 
queda de manifiesto en cada uno de 
ellos, lo que corrobora lo establecido por 
los índices l a y J a - 

En los mapas realizados se alcanzan 
a apreciar los diferentes focos de agrega- 
ción o de agrupamiento de las pobla- 
ciones del parásito en todas las parcelas 
estudiadas. Se distingue de forma 
general una cierta relación entre el 
número de centros de agregación exis- 
tentes y el promedio de muérdagos 
detectados en cada caso. Mayor cantidad 
de centros de agregación se relacionan 
con mayor densidad del parásito. Que 
fuera posible observar en cada uno de los 


103 

mapas realizados la presencia de los 
centros de agresión bien definidos, 
permite afianzar lo hallado con los índices 
del sadie, como se mencionó anterior- 
mente. 

Los numerosos focos de agregación 
detectados en cada casó, corroboran 
particularmente lo detectado por el índice 
J a . En todos los casos los centros de 
agregación del parásito se hallaron repar- 
tidos más o menos por toda la parcela 
experimental. Habría que determinar en 
estudios posteriores si la altitud y la 
ubicación de las parcelas tienen algún 
efecto sobre el comportamiento espacial 
detectado en las poblaciones del muér- 
dago enano. 

De manera general en los mapas 
elaborados con el método del Krigeado 
Ordinario se observó que la superficie 
libre de infestación disminuyó conforme la 
densidad media era mayor. De esta 
manera se detectó que en las parcelas 3 
y 5 fue donde existió la menor superficie 
sin infestación, por el contrario, en las 
parcelas 9 y 10 la superficie no infestada 


Tabla 4. Porcentaje de superficie infestada y no infestada por el muérdago enano 


MAPAS ELABORADOS CON EL KRIGEADO ORDINARIO 


PARCELAS 

% INFESTADO 

% NO INFESTADO 

P1 

82,6 

17,4 

P2 

75,8 

24,2 

P3 

94,1 

5,9 

P4 

87,1 

12,9 

P5 

92,9 

7,1 

P6 

75,1 

24,9 

P7 

69,5 

30,5 

P8 

70,1 

29,9 

P9 

65,4 

34,6 

PIO 

64,9 

35,1 
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Figura 1. Mapas de densidad de las poblaciones del muérdago enano elaborados con 
antelación con el método del Krigeado Ordinario, para cada una de las parcelas 

experimentales. 
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fue mayor. En el tabla 4 se detallan los 
resultados obtenidos del porcentaje de 
superficie infestada. En el caso de los 
mapas el porcentaje de superficie sin 
infestar se halló en el rango entre 35,1% 
y 5,9% del área total con un valor medio 
de 22,25%. De manera general se 
encontró que a mayor densidad menor 
superficie sin infestar. Lo descrito resulta 
interesante, ya que en ningún caso la 
infestación provocada por el parásito, 
alcanzo 100%. Lo anterior resulta de gran 
relevancia en el manejo integrado del 
muérdago enano, ya que permitiría 
realizar tácticas de control sobre las 
áreas específicas de infestación. 


Estabilidad espacial y temporal 

En la tabla 5 se concentraron los 
resultados de la comparación entre los 
mapas de la distribución espacial de las 
poblaciones del parásito, en cada parcela 
estudiada. Para llevar a cabo la compara- 
ción entre los diferentes mapas se utilizó 
el índice l m del sadie; si los valores de 
este índice eran mayores que 0, esto indi- 
caba que no existía diferencia significa- 
tiva entre los mapas comparados. 

Una vez realizadas las 45 compara- 
ciones posibles entre los mapas elabo- 
rados de las respectivas parcelas de 


Tabla 5. índices de asociación de distribución espacial comparando parcelas con 

muérdago enano. 


PARCELAS 


PARCELAS 


PARCELAS 


COMPARADAS 


COMPARADAS 


COMPARADAS 


1 vs. 2 

- 0,16 

2 vs.9 

- 0,32 

5 vs.6 

- 0,22 

1 vs.3 

- 0,35 

2 vs.10 

- 0,39 

5 vs.7 

- 0,88 

1 vs.4 

- 0,22 

3 vs.4 

- 0,67 

5 vs.8 

- 0,38 

1 vs. 5 

- 0,76 

3 vs. 5 

-0,53 

5 vs. 9 

- 0,11 

1 vs. 6 

- 0,31 

3 vs.6 

- 0,14 

5 vs.10 

- 0,40 

1 vs.7 

- 0,65 

3 vs.7 

- 0,69 

6 vs.7 

- 0,33 

1 vs. 8 

- 0,72 

3 vs.8 

- 0,12 

6 vs.8 

- 0,51 

1 vs.9 

- 0,95 

3 vs.9 

- 0,44 

6 vs.9 

- 0,92 

1 vs. 10 

- 0,27 

3 vs. 10 

- 0,95 

6 vs. 10 

- 0,80 

2 vs. 3 

- 0,54 

4 vs.5 

- 0,28 

7 vs.8 

- 0,26 

2 vs.4 

- 0,70 

4 vs.6 

-0,57 

7 vs.9 

- 0,19 

2 vs.5 

- 0,94 

4 vs.7 

- 0,93 

7 vs.10 

- 0,97 

2 vs.6 

- 0,53 

4 vs.8 

- 0,28 

8 vs.9 

- 0,37 

2 vs. 7 

- 0,74 

4 vs. 9 

- 0,76 

8 vs. 10 

- 0,20 

2 vs.8 

- 0,66 

4 vs.10 

- 0,48 

9 vs.10 

- 0,26 


Valores de lm > 0 indican asociación espacial 
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estudio, de acuerdo con el índice l m , no 
se halló asociación espacial alguna entre 
las comparaciones realizadas. Esto indica 
que no fue posible detectar una estabi- 
lidad espacio-temporal de las poblaciones 
de la plaga. 

DISCUSIÓN 
Estadística Clásica 

Los índices de dispersión utilizados 
en este trabajo mostraron un escenario 
general de agregación en el comporta- 
miento espacial del muérdago enano. 
Aunque el índice de dispersión (cociente 
media-varianza) y el índice de Green 
mostraron resultados similares en cuanto 
a la detección del agregamiento, la fuerza 
de dicha agregación la detectaron de dife- 
rente manera. En general el índice de 
Green encontró un mayor número de 
agregaciones leves. Si bien este indice es 
presentado como el más adecuado de los 
índices de dispersión, ya que es indepen- 
diente del número total de individuos 
capturados y del valor medio de su abun- 
dancia, no es independiente del tamaño 
de la muestra (De los Santos, 1982). Por 
su parte, el índice de dispersión (cociente 
varianza-media) al utilizar la varianza 
para evaluar el grado de agregación 
presenta graves inconvenientes, ya que 
depende directamente de la abundancia 
(Leveche, 1972). Lo que hará que pobla- 
ciones con densidades medias altas 
presenten valores altos del índice, tal 
hecho limita en gran medida la fiabilidad 
del resultado obtenido. 

El uso de distribuciones estadísticas 
presentó resultados similares a los de los 
índices de dispersión, con un ajuste casi 
general a la distribución Binomial Nega- 
tiva. Dicho modelo tiene una larga tradi- 
ción en la literatura ecológica para 
explicar el comportamiento espacial de 
los organismos (Pielov, 1977). 
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Asociado a la distribución Binomial 
Negativa tenemos el parámetro k, que ha 
sido muy utilizado como otro índice de 
agregación (Southwood, 1978 y Taylor, 
1984). Lo valores de k obtenidos en este 
trabajo se encuentran dentro del rango de 
valores generales para organismos vivos. 
Este índice, sin embargo, está influen- 
ciado por el tamaño de la muestra y de la 
unidad de muestreo (Colé, 1946; Morris, 
1954). Si se logran mantener constantes 
estos factores, este índice proporciona 
una media muy útil del nivel de agrega- 
ción de una población, siendo a la vez 
muy sensible al tipo de hábitat y al 
estadio de desarrollo de las poblaciones 
(Hairston, 1959; Waters, 1954). 

Los métodos estadísticos clásicos 
utilizados para establecer la distribución 
espacial de las poblaciones de orga- 
nismos vivos tienen deficiencias y limita- 
ciones. Tales deficiencias se han puesto 
de manifiesto en el presente trabajo. En 
muchas ocasiones los resultados obte- 
nidos con ellos parecían ser aparente- 
mente contradictorios, en otros casos no 
era posible ajustar los modelos de distri- 
buciones espaciales y en otras circuns- 
tancias se ajustaban los datos a los dos 
modelos de distribución. A ello hay que 
añadir otra gran limitación de estos 
métodos, como es que no tienen en 
cuenta la naturaleza bidimensional de las 
poblaciones, es decir, la localización 
exacta de los individuos muestreados. 
Por consiguiente, no utilizan toda la infor- 
mación espacial disponible, lo que limita 
su alcance. 

La distribución espacial agregada 
del muérdago enano fue estudiada por 
Gutiérrez (1970) y Andrade (1981); en 
ambos casos los autores utilizaron distri- 
buciones espaciales en sus estudios y 
registraron algunos problemas con el 
ajuste estadístico correspondiente, como 
sucedió en este estudio. Rodríguez 
(1985) y Hernández Montiel (1999) deter- 
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minaron el mismo comportamiento espa- 
cial del parásito utilizando diversos 
índices de dispersión y encontraron 
asimismo, problemas en la determinación 
de la distribución espacial del muérdago 
enano. 

Análisis Espacial por índices de 
Distancia (sadie) 

Este tipo de índices tienen la gran 
ventaja sobre los estudiados con anterio- 
ridad, que tiene en cuenta la localización 
espacial de la muestra, además de ser 
más intuitivos y biológicos que los índices 
no espaciales (Perry, 1995a,b). Los resul- 
tados obtenidos con estos índices (l a y J a ) 
mostraron en todos los casos una estruc- 
tura espacial en agregados de las pobla- 
ciones del muérdago enano, por lo que 
demuestran una mayor eficacia que los 
índices no espaciales. El índice J a 
también indica si la estructura espacial 
agregada se encuentra dispuesta en uno 
o varios centros de agregación. En este 
trabajo se encontró en la mayoría de los 
casos la existencia de varios centros de 
agregación. 

También en el caso de la utilización 
del sadie para establecer el tipo de distri- 
bución de este parásito sería interesante, 
en trabajos posteriores, analizar el patrón 
espacial con datos a diferentes escalas y 
a través de diversos periodos de tiempo, 
tal y como sugieren de forma general 
para el comportamiento espacial de orga- 
nismos vivos Perry (1998), Ferguson et 
al., (2000) y Thomas et al., (2001). De 
esta manera se podrían comparar los 
resultados obtenidos y conocer con 
mayor detalle, en este caso, la distribu- 
ción del muérdago enano y las bondades 
del método. Hay que destacar el hecho 
de que las diferentes altitudes en las 
cuales fueron establecidas las parcelas 
de estudio y el diámetro del tronco del 
pino, no influyeron en el establecimiento 


de la agregación por parte de las pobla- 
ciones de muérdago enano, tal y como se 
aprecia en los mapas. 

Aunque no se determinó estadística- 
mente, se pudo observar en campo que, 
tomando en cuenta que el diámetro del 
tronco de los pinos está en relación con 
su edad, el parásito no parece tener 
preferencias de infestación en cuanto a la 
edad y talla del pino. Lo anterior resulta 
interesante, ya que Gutiérrez (1970) 
señala que la infestación del muérdago 
enano influye sobre el desarrollo en 
grosor del fuste del pino, por lo tanto, el 
ataque a pinos jóvenes reduciría el 
desarrollo posterior de estos árboles. 
Tomando en cuenta que una parte consi- 
derable de pinos muestreados en este 
estudio eran relativamente jóvenes, el 
impacto sobre su desarrollo tendrá gran 
importancia dentro del contexto ambiental 
del pnnt. La anterior aseveración habría 
que corroborarla estableciendo un 
estudio en el que se demostrara estadís- 
ticamente que el parásito infesta de forma 
generalizada sin importar edad y talla del 
pino. 

Cabe señar que en la ladera sur del 
Parque Nacional Nevado de Toluca 
(cuenca del Balsas), se han registrado 
resultados similares en cuanto a la distri- 
bución espacial del muérdago enano; en 
el caso mencionado se utilizaron también 
los índices del sadie (Ramírez-Dávila y 
González-Andujar, 2006). Con lo anterior 
se puede afirmar que el comportamiento 
espacial de este parásito en la totalidad 
del Parque Nacional Nevado de Toluca es 
agregado y distribuido en varios centros 
de agregación. 

Los índices del sadie se han utili- 
zado con éxito en otros estudios para 
establecer la distribución espacial de 
organismos biológicos: Ferguson et al., 
2000; Thomas et al., 2001; Winder et al., 
2005; Alexander et al., 2005; Conrad, et 
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al., 2006 y Ramírez-Dávila y González- 
Andujar, 2006. 

Superficie infestada 

Pese a que el muérdago enano se 
encontraba en condiciones óptimas para 
su desarrollo, el porcentaje de infestación 
alcanzado en las diferentes parcelas de 
estudio nunca llegó a 100%. Por lo tanto, 
las acciones de control deberán ser diri- 
gidas sobre zonas específicas de infesta- 
ción y no sobre la suposición de una colo- 
nización de 100% por parte del parásito, 
tal y como proponía en su momento 
González (1989), con la utilización de 
herbicidas selectivos. Es importante 
resaltar que, contando con las áreas 
específicas de infestación de este pará- 
sito sobre la zona de estudio, se podrían 
establecer programas de manejo de este 
problema sanitario del pnnt en su ladera 
norte, utilizando técnicas de precisión, 
como por ejemplo la adecuación de las 
propuestas sugeridas por Fleischer et al., 
1999. Con lo anterior se podrían reducir 
los costos económicos y medioambien- 
tales, al reducir el uso de herbicidas 
selectivos en el control de esta plaga. 

Estabilidad espacial y temporal a corto 
plazo 

Se esperaba que existiera alguna 
estabilidad espacio-temporal en las 
poblaciones de muérdago enano, pero 
los resultados indicaron lo contrario. 
Habría que establecer si la altitud o la 
edad de los árboles influye en la manera 
de agregarse por parte de las pobla- 
ciones del parásito, de ser así, ésto expli- 
caría por qué la manera de establecerse 
los agregados dentro de las parcelas de 
estudio fue diferente. Otra posibilidad 
podría ser, que el tamaño de las parcelas 
experimentales no permite detectar simili- 
tudes entre la forma de agregarse de las 
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poblaciones de muérdago enano, por lo 
tanto, estudios posteriores podrían 
explorar la posibilidad de aumentar el 
tamaño de las parcelas para determinar si 
se puede detectar o no esa estabilidad. 
Como se mencionó anteriormente, la 
distribución agregada del parásito ya fue 
estudiada por Gutiérrez (1970), Andrade 
(1981), Rodríguez (1985) y Hernández 
(1994). Sin embargo, estos autores repor- 
taron problemas con los métodos de 
Estadística No Espacial. En este estudio, 
el sadie permitió sin lugar a dudas deter- 
minar el comportamiento agregado de las 
poblaciones del muérdago enano. 

CONCLUSIONES 

1 . Los métodos no espaciales 
presentaron serias deficiencias y limita- 
ciones al tratar de determinar el compor- 
tamiento de las poblaciones de muérdago 
enano. 

2. Los análisis con sadie presen- 
taron un comportamiento espacial en 
agregación, con las poblaciones distri- 
buidas en varios focos de concentración. 

3. Los índices l a , J a e l m demostraron 
ser eficientes en la determinación de la 
estructura espacial de las poblaciones de 
muérdago enano. 

4. No se logró detectar estabilidad 
espacial y temporal en las poblaciones de 
la plaga. 

5. Se identificaron áreas libres de 
infestación que permitían en su momento 
una aplicación dirigida de las medidas de 
control establecidas para este parásito. 
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